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1.  INTRODUCTION 


1.1  PREAMBLE 

There  are  many  natural  and  man-made  events  which  subject  soils 
to  large  dynamic  forces  and  the  survival  of  buildings  in  the  vicinity  is 
strongly  dependent  upon  one  manner  in  which  these  soils  respond.  To 
realistically  assess  vulnerability  it  is  therefore  essential  to  have  a 
good  quantitative  knowledge  of  soil  dynamic  response.  Of  particular  interest 
to  the  US  Air  Force  is  the  response  of  soil  to  the  loads  associated  with  a 
nuclear  explosion.  An  accurate  knowledge  of  this  response  is  required  for 
the  design  and  assessment  of  many  military  structures  such  as  command 
centers,  missile  silos  and  other  surface  or  under-ground  facilities. 

Modern  numerical  tools  allow  performing  sophisticated  dynamic 
analyses  which  model  the  interaction  of  soils  and  structures.  Indeed, 
the  US  Air  Force  is  presently  conducting  an  effort  to  improve  its 
numerical  capabilities  in  regard  of  that  specific  problem,  both  in  two 
and  three  space  dimensions.  In  fact,  the  numerical  problem  can  be 
considered  almost  solved.  Most  of  the  required  capabilities  already 
exist  and  the  effort  remaining  is  essentially  that  of  blending  them  into 
a  single  computer  program  of  the  required  generality. 

Of  greater  concern  is  the  difficult  problem  of  obtaining  an 
accurate  constitutive  description  of  the  soil.  To  the  writers’  knowledge, 
there  is  not  sufficient  confidence  in  any  of  the  existing  soil  constitutive 
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laws  to  allow  making  reliable  predictions  of  the  behaviour  of  soil  when 
subjected  to  complex  stress  paths  under  dynamic  condit.ons.  Even  at 
the  fundamental  level,  it  is  not  obvious  whether  concepts  such  as  the 
traditional  "effective  stress"  still  apply  under  the  high  transient 
stresses  generated  by  blast  loading. 

It  is  fair  to  state  that  in  the  last  two  decades  numerical 
capabilities  have  improved  at  a  much  faster  rate  than  soils  testing. 

At  present,  most  soils  testing  is  conducted  in  the  laboratory  on  relatively 
small  soil  samples.  These  samples  suffer  from  the  unknown  effects  of 
disturbance  and  stress  relief  which,  due  to  the  strong  history  dependence 
of  soils,  often  introduce  large  uncertainties  in  the  validity  of  the 
results  obtained  by  those  tests.  A  clear  example  of  this  is  the  large 
difference  in  shear  modulus  which  is  normally  observed  at  identical 
strain  levels  between  laboratory  resonant  column  tests  and  in-situ 
seismic  measurements. 

There  are  many  other  uncertainties  inherent  in  present  day 
soils  testing  methods.  Variables  can  only  be  measured  at  internal  or 
external  boundaries  in  the  material,  the  influence  of  which  is  not 
always  obvious.  The  effect  of  sample  size  and  soil  inhomogeneity  is 
difficult  to  assess,  and  still  more  important,  testing  has  normally  only 
one  or  two  (seldom  three)  degrees  of  freedom,  whereas  a  true  constitutive 
description  requires  at  least  six  components  of  the  stress  and  strain 
tensor.  In  the  absence  of  better  data  one  must  make  assumptions  to 
generalize  the  results  of  tests  such  as  triaxial  and  simple  shear  to  the 
full  three-dimensional  situation. 
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The  conclusion  is  that  too  much  is  often  expected  from  the 


tests.  It  is  equally  clear  that  it  is  not  a  simple  matter  to  overcome 
these  problems.  The  work  described  in  this  report  will  hopefully  constitute 
a  step  towards  a  better  material  description  by  providing  a  method  which 
minimizes  some  of  the  problems  of  the  present  testing  procedures  for 
derivation  of  the  dynamic  shear  modulus  of  soils. 


1.2  PURPOSE  OF  THE  WORK  PERFORMED 


Dynamic  calculations  in  soil  are  very  often  based  on  generalisations 
or  direct  use  of  a  one-dimensional  constitutive  law  which,  in  shear  is 
characterised  by  a  secant  shear  modulus,  expressed  as: 


_y_ 


V 


r 


(1) 


where  U  is  the  secant  shear  modulus  (p  =  x/e,  where  x  is  the 
shear  stress) 


o’  is  the  effective  confining  pressure 

e  is  the  shear  strain 

n  is  a  constant 

the  subscript  r  applies  to  values  of  the  variables  at  the  reference 

state  (’Characterized  by  small  strains  and  o’  =  o’  ,  is  a  reference 

rf 

effective  confining  pressure) . 


The  objective  of  the  apparatus  described  in  the  present  report 
is  to  produce  the  function  y(e)  along  the  soil  profile  for  the  complete 


range  of  strains  of  interest.  The  function  p (e)  is  intended  to  be 
representative  of  the  in-situ  soil  conditions;  hence,  an  attempt  has 
been  made  to  minimize  the  soil  disturbance  during  the  test. 

1.3  SCOPE  OF  WORK 

In  trying  to  provide  a  basis  towards  the  operation  of  a  device 
able  to  perform  the  aforementioned  task,  the  scope  of  the  present  work 
can  be  summarized  in  the  following  items: 

(a)  A  brief  review  of  the  state-of-the-art  as  regards  the 
determination  of  the  dynamic  shear  modulus  of  soils. 

(b)  A  feasibility  analysis  of  a  borehole  device  for  in-situ  determination 
of  the  dynamic  shear  modulus. 

(c)  Design,  construction  and  laboratory  testing  of  the 
mechanism  which  couples  the  device  to  the  borehole  wall. 

The  feasibility  analysis  includes  several  different  aspects: 

(I)  Conceptual  design  of  the  overall  device  and  its  operation. 

(II)  Analytical  and  numerical  studies  to  determine  the  stress, 
strain  and  displacement  fields  at  the  strain  levels  of 
interest;  also,  studies  to  predict  the  normal  operation  of 
the  device  as  well  as  the  influence  on  its  performance  of 
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factors  such  as  soil  layering  and  soil  disturbance  around  the 
borehole. 


(Ill)  The  construction  of  a  procedure  to  obtain  the  shear  modulus  - 

shear  strain  curve  in -situ,  based  on  the  values  of  the  variables 
during  operation  in  the  field. 


It  must  be  however  remembered  that  the  device  proposed  in  this 
report  ic  r.ot  intended  to  solve  the  problem  of  measuring  the  shear 
modulus  of  soil  completely,  but  simply  to  alleviate  some  of  the  most 
serious  inconsistencies  of  present  practice. 


2.  STATE-OF-THE-ART 


Dynamic  analyses  of  soil  response  are  frequently  required  for  a 


variety  of  problems,  including  blasting,  earthquakes,  etc.  Such  analyses 
must  always  be  based  on  the  dynamic  properties  of  the  soil.  The  present 
chapter  briefly  reviews  those  properties  and  the  current  methods  of  obtaining 
them. 


The  information  which  is  ideally  required  for  performing  a  dynamic 
soil  analysis  consists  of  the  full  constitutive  equations  of  the  soil.  This 
is  too  formidable  a  task  to  be  attempted  as  yet.  Idealised  forms  of  the 
constitutive  laws,  which  approach  the  soil  behaviour  over  the  range  of  conditions 
of  interest,  must  therefore  be  assumed  in  order  to  decide  which  properties  are 
to  be  determined. 


By  passing  the  completely  linear  approaches,  which  are  not  of  great 
interest  here,  the  first  non  linear  approximations  which  would  be  applied  with 
some  generality  were  based  on  the  assumption  that  a  cne-dimensional  propagation 
of  shear  waves  would  be  an  adequate  representation  of  soil  dynamics  problems 
(Schnabel  et  al,  1972''*’.  Under  these  circumstances,  the  soil  is  totally  defined 
by  the  relation: 


u(e  )  e  +  n(e  )  e 
e  e 


where  t  is  the  shear  stress  (*) 


e  is  the  shear  strain,  e is  the  rate  of  increase  of  shear  strain 
M  is  a  strain-dependent  secant  shear  modulus 


n  is  a  strain-dependent  viscosity 


ee  is  an  equivalent  strain 


(*)  All  symbols  are  defined  in  their  first  appearance  and  in  Appendix  A. 


-  . .  .v,  -v. 


Such  an  approach  has  constituted  the  basis  of  many  lines  of  \ 
analysis,  primarily  due  to  the  fact  that  it  can  be  formulated  in  an  \ 

"equivalent  linear"  manner  by  using  for  the  values  of  the  modulus  and  the 
vicosity  those  corresponding  to  a  certain  percentage  of  the  maximum  shear 
strain  (equivalent  strain)  developed  at  each  point  in  the  soil  profile.  The 
approach  has  been  generalized  to  more  dimensions  and  other  alternative  procedures 
have  been  devised;  but  the  determination  of  the  function  p(e)  still  remains 
a  basic  concern  in  practically  any  soil  dynamics  analysis. 

The  function  p(e)  is  known  to  be  more  or  less  constant  for  strain 

levels  smaller  than  10  ^  or  even  almost  10  At  this  point,  non-linearities 

appear  and  the  shear  modulus  starts  decreasing;  clearly,  it  should  tend  to 

zero  as  the  strain  level  increases  indefinitely  (see  for  example:  Seed  and 
2 

Idriss,  1970)  . 


The  standard  procedure  for  determining  p(e)  at  present  is  explained 
below.  Geophysical  techniques  are  used  to  measure  the  in-situ  velocity  of 
compressional  and  shear  waves.  Those  two  velocities  can  then  be  used  together 
with  the  well  known  elastic  relationships  to  derive  the  value  of  the  dynamic 
Poisson's  ratio  and  of  the  shear  modulus  for  the  lower  range  of  strains. 

(2) 


V  =  pv2 
s 


(v  /v  )  -  2 

p  s 


(3) 


v  = 


2[W2  -i] 


where 


p  is  the  density 

v  is  the  dynamic  Poisson's  ratio 

Vp'vs  are  the  velocities  of  propagation  of  compressional  and  shear 


waves  respectively. 


In  materials  which  are  suitable  for  sampling,  laboratory  tests 

are  then  conducted  in  conditions  which  attempt  to  reproduce  the  in-situ 

configuration-  Shear  tests  of  the  resonant  column  type  are  used  for  strain 
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levels  up  to  10  or  10  .  Those  are  complemented  either  with  simple  shear 

tests  or  with  strain-controlled  dynamic  triaxial  tests  for  the  upper  range 
of  strains,  including  failure  of  the  material  if  necessary.  In  case  of  using 
triaxial  tests,  the  property  obtained  is  the  Young's  modulus;  assuming  that 
the  dynamic  Poisson's  ratio  is  constant  and  hence  independent  of  the  strain 
level,  the  value  of  the  Poisson's  ratio  derived  in  the  field  is  then  used 
together  with  the  Young's  modulus,  E,  in  order  to  derive  the  value  of  the 
shear  modulus  sought. 


2 (1+v) 

The  above  procedure  presents  a  number  of  conceptual  inconsistencies 

due  to  the  three  diverse  experimental  philosophies  embodied,-  as  it  could  be 

expected,  such  inconsistencies  are  clearly  displayed  by  the  results  obtained. 

At  the  overlapping  range  of  strains,  field  and  laboratory  techniques  have 

been  found  to  yield  values  of  the  shear  modulus  which  differ  by  factors 

varying  from  1.2  to  4.0  (Stokoe  et  al,  1978)^.  Such  differences  are  believed 

to  be  due  to  the  disturbance  which  occurs  during  sampling  and  to  the  stress 

relief  undergone  by  the  sample  when  removed  from  its  native  environment;  in 

general,  they  show  the  difficulty  of  simulating  field  conditions  in  a 

4 

laboratory  test  (Stokoe  and  Richart,  1973)  .  A  number  of  procedures  have 
been  devised  to  correct  the  laboratory  curve  so  that  it  passes  through  the  field 
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value  of  the  modulus  :  both  parallel  and  proportional  corrections  have  been 
advocated  (Anderson  et  al,  1978)"’.  Theories  have  also  been  proposed  which 
postulate  that  a  thixotropic  effect  stiffens  the  soil  with  time  until  an  energy 
balance  within  the  soil  is  attained  (Anderson  et  al,  1976)  .  Such  procedures 
and  theories  are  valuable  but,  used  in  this  context,  they  simply  constitute 
attempts  to  patch  up  the  weakness  inherent  in  the  experimental  tricotomy  which 
the  method  suffers. 

Let  us  examine  in  more  detail  some  of  the  procedures  used  for 

in-situ  determination  of  the  shear  modulus  of  soils.  The  most  common  ones 

are  based  on  seismic  wave  propagation  phenomena.  The  instrumentation 

required  for  those  tests  is  basically  composed  of  a  seismic  wave  source  and 

one  or  more  pick-up  devices;  nature  and  location  of  these  components  requires 

a  careful  understanding  of  wave  generation  and  propagation  phenomena.  General 

theory  of  propagation  of  waves  in  elastic  media  can  be  found  in  Timoshenko 
7  8 

and  Goodier  (1951)  and  Ewing  et  al  (1957)  .  Since  the  material  behaviour 
remains  essentially  elastic  during  most  wave  propagation  tests,  these  general 
studies  provide  the  necessary  background. 

Seismic  tests  in  which  source  and  detectors  are  located  on  the 

surface  of  the  soil  were  developed  first  and  have  been  discussed  by  Duke 
9 

(1969)  .  Methods  based  in  reflection  and  refraction  techniques  were  described 
by  Richart  et  al  (1970)1-0.  In  any  case,  the  input  source  can  be  impulsive  or 
vibratory.  In  case  of  impulse  tests,  the  energy  source  is  provided  by 
detonation  of  explosives  or  by  mechanical  means.  The  arrival  times  of  the 
direct  and/or  reflect-!  and/or  refracted  body  or  surface  waves  are  measured 
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in  an  array  of  transducers  and  used  to  determine  average  values  of  wave 


velocities  alony  the  different  paths  assumed.  Refraction  surveys  give 
in  addition  the  depths  to  the  different  refractors,  which  can  be  used  to 
check  the  reliability  of  the  results  obtained. 

The  impulse  source  is  chosen  depending  on  the  type  of  waves 
desired.  Shear  waves  are  the  most  suitable  ones  for  determination  of  the 
shear  modulus  because  their  velocity  is  directly  related  to  it.  Also, 

r 

at  high  degrees  of  saturation,  the  compressional  wave  is  likely  to  represent 
more  the  velocity  of  propagation  through  the  water  and,  therefore,  produces 
little  information  on  the  behaviour  of  the  soil  skeleton. 

A  description  of  various  suitable  sources  of  shear-waves  is  given 
U 

by  Mooney  (1974)  .  Impulse  methods  often  produce  considerable  amounts  of 

P-  waves  and  consequently  S-wave  arrivals  are  difficult  to  identify;  this 

problem  can  be  somewhat  overcome  by  using  a  seismic  source  which  allows 

12 

reversal  of  the  energy  input.  Schwarz  et  al  (1974)  used  this  technique 
to  determine  the  shear  wave  velocity  in  marine  foundations.  The  procedure 
consists  on  subtracting  seismic  records  produced  with  opposite  energy  inputs; 
the  reason  is  that,  under  those  conditions,  P-wave  components  should  be 
equal  .vhile  S-wave  components  should  be  exactly  opposite.  Lande  (1974)^ 
describes  a  technique  in  which  the  reaction  force  on  a  horizontal  plate 
produced  by  a  water  cannon  is  used  as  a  source  of  horizontally  polarized  waves. 
Different  orientations  of  the  water  cannon  allow  production  of  P-  and/or  S- 
waves  as  veil  as  application  of  the  phase  reversal  method. 
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Surface  tests  with  sustained  seismic  sources  have  been  performed 
using  constant  or  systematically  varying  sinusoidal  inputs.  The  source  is 
normally  a  mechanical  or  electro-mechanical  vibrator  and  the  signals  are 
recorded  at  various  points  at  or  below  the  surface  of  the  soil.  A  procedure 
has  been  described  by  Kurzeme  (1971^in  which  horizontally  polarized  shear 
waves  are  produced  by  means  of  a  torsional  vibrator.  The  method  which  has 
been  primarily  used  in  pavement  investigation,  consists  in  analysing  the 
measured  velocity  dis.^rsicn  with  frequency  using  an  idealised  material 
and  structural  nxxL-1.  Shear  wave  velocity  and  thickness  of  the  various 
layers  can  be  determined;  the  quality  of  the  results  depends  on  how  closely 
the  model  simulates  the  real  soil  conditions. 

All  of  the  aforementioned  tests  utilise  body  wave  velocities 
to  determine  the  shear  modulus  of  soils.  A  method  using  Rayleigh  waves  was 
described  by  Jones  (1960^51962)16.  This  technique,  which  was  initiated  in 
Germany  during  the  thirties,  consists  in  measuring  the  wave  length  and  phase 
velocity  of  the  vibrations  propagated  along  the  surface  of  the  soil.  Shear 
wave  velocity  can  then  be  calculated  assuming  a  ratio  of  v  to  v  resulting 
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from  elastic  theory  (Richart  et  al  1970)  .  For  the  range  of  Poisson's  ratios 

commonly  found  in  soil  materials,  errors  due  to  the  above  assumption  are  less 

than  10%.  Due  to  the  fact  that  the  shear  modulus  of  soil  is  not  constant  with 

depth,  the  phase  wave  velocity  -  depends  on  the  frequency  cf  vibration.  Fowler 
17 

(1972)  has  found  empirically  that  the  velocity  determined  from  propagation  of 
a  surface  wave  is  applicable  to  a  depth  of  equal  to  half  of  the  wave  length. 


A  method  using  this  technique  is  described  by  Ballard  et  al,  (1975) 
in  which  an  impulsive  source  is  used  to  generate  Rayleigh  waves.  Instrumen¬ 
tation  and  field  procedures  are  similar  to  those  used  in  the  refraction  method; 
phase  velocity  can  be  determined  by  numbering  the  peaks  and  troughs  of  the 
oscillations  and  following  them  through  adjacent  traces-  This  method  is 
limited  in  application  by  the  inability  to  exercise  selective  control  over  the 
generated  frequencies  and,  hence  to  the  depths  associated  with  those  frequency 
values. 


Various  methods  have  been  developed  to  measure  the  direct  shear 

wave  velocity  between  surface  and  subsurface  points.  They  normally  utilize 

a  single  borehole;  the  seismic  source  is  located  on  the  surface  of  the 

soil  (downhole)  or  in  the  borehole  (uphole) .  In  the  first  case,  transducers 

are  placed  in  the  borehole;  in  the  second,  they  are  at  the  surface  and/or 

in  the  borehole.  Several  types  of  seismic  sources  have  been  used  and  references 

19 

are  given  by  Stokoe  et  al,  (1972) 


The  crosshole  seismic  technique  is  a  method  which  has  received 

a  great  deal  of  attention  in  the  last  few  years.  It  basically  consists  in 

measuring  the  tine  that  body  waves  take  to  travel  along  a  horizontal  path 

between  two  points  in  the  soil  mass.  The  source  and  one  or  more  transducers 
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are  located  in  separate  boreholes.  Swain  (1962)  and  Stokoe  et  al,  (1972) 

used  the  technique  with  impulsive  notions  generated  by  impact.  Tanimoto 
21 

et  al,  (1973)  discussed  the  problems  of  using  only  two  boreholes;  those  arise 
due  to  the  differences  between  the  characteristics  of  the  motions  of  the 
transducer  in  the  operati- „  source  and  those  of  the  transducer  in  the  receivins 
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hole  as  well  as  from  large  differences  in  the  strain  level  along  the  wave 

22 

path.  A  low  energy  impulse  source  has  been  developed  by  Auld  (1977);  it  is 
composed  of  a  stationary  part  and  a  hammer,  the  stationary  par  consists 
of  a  hydraulic  cylinder  block  with  eight  horizontal  pistons,  four  pistons 
for  expansion  to  grip  the  borehole  walls  and  the  others  for  contraction. 

The  hammer  can  strike  either  downwards  or  upwards  thus  allowing  application 
of  the  phase  reversal  technique. 


The  use  of  vibratory  sources  of  shear  waves  in  crosshole  configurations 

has  been  tried  by  Ballard  et  al  (1969 f\nd  it  has  teen  successfully  developed 
18 

by  Ballard  (1975) ,  It  consists  of  a  vibratory  source  in  conjunction  with  a 
spring  wedged  device  in  contact  with  the  borehole  wall.  Vertically  polarized 
shear  waves  are  propogated  and  the  results  obtained  are  processed  for  determining 
the  shear  wave  velocity. 


All  methods  mentioned  so  far  are  used  for  determining  the  in-situ 
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shear  modulus  for  small  strains  (up  to  10  ) .  Some  of  them  only  generate 

small  strains;  in  cases  where  large  strains  are  generated,  the  configuration 
of  the  test  is  chosen  so  that  the  large  strains  do  not  influence  the 
measurements.  The  measurement  of  the  shear  modulus  in-situ  for  large  strains 
has  proved  to  be  a  difficult  task  indeed.  Large  strain  amplitudes  can 
be  obtained  using  explosives  as  seismic  sources.  The  main  problems  associated 
with  explosives  is  that  they  generate  large  amounts  of  P-waves  and  that  the 
results  produced  are  difficult  to  analyze  and  non-repeatable. 


Some  investigators  have  tried  to  produce  large  strains  by  means  of 
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vibratory  sources.  Miller  et  al,  (1972)  utilized  a  vibratory  high-energy 
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driving  source;  velocity  transducers  were  located  in  two  adjacent  holes. 
Shear  waves  originated  from  the  sinusoidal  input,  are  propagated  to  the 
monitoring  boreholes.  Although  the  source  was  initially  designed  to  produce 
large  strains,  it  has  been  used  only  for  small  strains  due  to  difficulties 
of  interpretation  of  the  signals  related  to  interference  with  background 
noise. 


An  attempt  to  measure  the  shear  modulus  as  a  function  of  strain 
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level  using  an  in-situ  technique  has  been  presented  by  Miller  et  al,  (1975) 

26 

and  Troncoso  (1975)"  .  The  method  essentially  consists  in  a  crosshole  impulse 

test;  the  generating  source  is  composed  of  an  anchor  firmly  coupled  to  the 

surrounding  soil  and  a  striking  hammer  which  impacts  the  top  of  the  anchor. 

Coupling  of  the  anchor  to  the  soil  is  accomplished  by  means  of  an  radial  stress 

The  recording  stations  consist  in  velocity  transducers  firmly  coupled  to  the 

soil  by  means  of  high  pressure  rubber  balloons.  The  shear  wave  velocity 

between  two  stations  is  calculated  from  arrival  times  of  wave  peaks. 

The  strain  level  developed  at  each  station  is  calculated  using  the  elastic 

relationship  which  equals  it  to  the  ratio  of  the  measured  particle  velocity 

divided  by  the  shear  wave  velocity.  The  procedure  is  designed  to  work  up 

to  strain  levels  of  about  10  ^  although  it  does  not  appear  to  have  been 
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used  in  natural  soils  beyond  10 

Several  important  problems  are  associated  with  the  method  as  reported 
The  large  strains  occur  close  to  the  source.  At  this  point,  particle  velocity 
is  equated  to  source  velocity;  this  often  leads  to  errors  due  to  slip  which 
is  difficult  to  avoid  and  ’as  greater  tendency  to  appear  as  the  strain  and 


hence  the  input  energy  increase .  Fui  thermore ,  the  device  generates  P-waves 
which,  close  to  the  source,  have  not  yet  separated  from  the  S-waves,  thus 
introducing  uncertainties  in  interpretation.  The  formula  for  calculating 
local  strain  values  also  starts  introducing  errors  as  the  material  behaves 
more  non-linear ly,  which  is  the  behavior  of  greater  interest.  Finally, 
it  must  be  noted  that  wave  forms  are  preserved  only  in  linear  materials; 
non-linear  behavior  distorts  the  wave  forms.  Hence,  wave  parameters 
calculated  from  arrival  times  of  wave  peaks  are  likely  to  introduce  errors 
increasing  with  the  strain  level. 

Independently  of  the  type  of  borehole  test  for  high  strains,  it  can 
easily  be  seen  that,  even  under  linear  propagation,  strains  will  decrease 
as  Bessel  functions;  as  a  consequence  of  the  cylindrical  symmetry  this  decay 
will  be  progressively  faster  as  greater  non-linearities  (or  greater  strains) 
are  induced  in  the  soil.  If  it  is  also  realised  that  the  soil  has  only  a 
limited  shear  strength  around  the  borehole  which  can  be  activated  to  propagate 
the  input  to  the  progressively  larger  regions  of  surrounding  soil,  it  is 
quite  clear  that,  in  most  realistic  situations,  the  zone  of  large  strains 
is  going  to  be  confined  to  the  immediate  vicinity  of  the  source.  This  is 
confirmed  by  the  analyses  presented  in  this  report.  Any  type  of  crosshole 
methods,  even  with  very  short  distances  between  the  source  and  transducers 
(which  creates  many  other  experimental  problems)  are  likely  to  have  difficulties 
arising  from  vast  non-uniformities  in  wave  propagation  arising  from  the 
widely  varying  strain  levels.  Also,  any  method  relying  on  arrival  times  of 
any  kind,  is  likely  to  introduce  errors  associated  with  the  different 
velocities  of  propagation  of  different  strain  levels.  Because  of  the  above 
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3.  ANALYTICAL  SOLUTIONS 


It  seems  quite  clear  that  the  full  problem  of  a  finite  cylinder 
rotating  in  an  inhomogeneous  and  non-linear  domain  cannot  be  solved  by 
present  analytical  procedures.  However,  since  close-form  solutions  are 
very  useful  for  studying  the  influence  of  the  different  parameters,  it 
has  been  considered  appropriate  to  introduce  certain  simplifications  to 
the  problem  in  order  to  allow  the  derivation  of  a  number  of  analytical 


solutions. 


The  first  simplification  included  is  that  the  length  to  diameter 
ratio  of  the  shoe  (rotating  cylinder)  is  large.  This  assumption  effectively 
decreases  the  order  of  the  problem  from  two  spatial  dimensions  to  one;  the 
time  dimension  is  independent  of  that  simplification. 


Second,  as  the  maximum  frequencies  foreseen  for  the  device  are  on 
the  range  of  lO  to  20  Hz  and  typical  shear  wave  velocities  will  span  from  one 
to  a  few  hundred  metres  per  second,  it  appears  that  associated  wave  lengths 
will  be  substantially  larger  than  the  domain  under  investigation.  In  that 
situation,  the  static  strain  distribution  cannot  be  expected  to  differ 
greatly  from  the  real  part  of  the  dynamic  one. 


Finally,  the  soil  will  be  considered  as  a  homogeneous  half-space; 
although  linearity  and  full  isotropy  are  not  required  at  this  stage. 


axisymmetric  property  distributions  are  assumed  throughout  this  report. 
Not  all  of  the  above  .-.sumptions  are  needed  for  allowing  derivation 
of  dost  -form  solutions,  as  will  be  seen  later.  However,  those 
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obtained  will  be  presented  here  following  an  approximate  order  of  increasing 
complexity  or,  equivalently,  of  decreasing  number  of  assumptions. 


3.1  LONG  CYLINDER  -  STATIC 


Let  us  consider  an  infinitely  long  cylinder  of  radius  r^  embedded 


in  a  homogeneous  continuum.  Let  be  the  static  shear  stress  transmitted 


at  the  interface  between  the  cylinder  and  the  continuum.  The  equilibrium 


of  any  annulus  boundred  by  r^  and  an  arbitrary  radius  r  requires  (Figure  1) 


r(r)  = 


It  is  then  clear  that  the  only 


non-vanishing  strain  e  .  t called  e 


for  simplicity)  fulfills,  independ¬ 


ently  of  the  constitutive  relation: 


=  Tkt) 


Figure  1.  Shear  stress  distribution 


where  i  =  tit)  is  the  constitutive  relation  of  the  materi-1!  simple-shear . 


Equation  (2)  allows  the  derivation  of  the  strain 


distribution  under  fairly  general  constitutive  laws.  Five  such  laws 


have  been  used  as  examples.  The  laws  are  summarized  in  Table  1, 


together  with  their  formulation  and  main  characteristics: 


Ufa 


Table  l.  Constitutive  laws  used  in  the  derivation  of  strain  fields 


ELASTICITY 


HYPERBOLIC  LAW 


POWER  LAW 


HYPERBOLIC  TANGENT 


RAMBERG-OSGOCD 


FORMULATION 


t  =  u  e 
o 


1  +  y  e/x 
o  m 


t  =  Ae  (o<n<l) 


y  e 

t  =  t  tanh  - 

m  i 

m 


e  =  At  +  t/y  (n>i) 
o 


INITIAL 

MODULUS 


LIMITING 

STRESS 


For  each  one  of  the  above  laws,  the  corresponding  strain  distributions 


are  presented  below: 


a)  Elasticity 


The  constitutive  law  is: 


t  =  y  e 
o 


The  strain  distribution  becomes: 


x  r2 
o  o 

£  =  - 2~ 

V 


b)  Hyperbolic  Law 


Substituting  the  constitutive  law  into  equation  (2) : 
*-2 


1  +  —  c 
T 

m 
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Hence,  the  strain  distribution  is: 


e  = 


IV  To1 

“o  -2”  -  r 

r  m 

L  o  -I 


(4) 


or,  for  easier  comparison  with  equation  (3) : 

2 


e  = 


T  r  i 

o  o  1 


t  r 


r2  -  °  ° 


which  reverts  to  the  elastic  distribution  for  large  or  at  large  distances 


for  the  cylinder. 


<_■ )  Power  Law 


Operating  as  above: 
2 


=  Aen  (n<l) 


e  = 


2 

t  r 
o  o 


1  1/n 


Ar 


(5) 


which  approaches  the  elastic  distribution  as  n  approaches  unity. 


d)  Hyperbolic  Tangent 
2 


Vo 


r  tanh 


V 


m 

T 

m 

I* 

a  tanh 

r  2 
Vo 

Vo 

2 

T  r 
m 

(6) 


Again  the  elastic  distribution  is  recovered  at  large  distances  or. 


in  any  case,  for  large  *  - 

m 
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The  equation  proposed  is  a  typical  skeleton  curve  used  in  studies 
of  cyclic  behaviour. 


The  strain  distribution  becomes: 


Elastic  behavior  is  obtained  as  n  approaches  1  or  A  approaches  zero 
or  the  distance  becomes  large. 


In  order  to  study  the  effect  of  the  material  description  in  the 
strain  distribution,  a  typical  stress-strain  curve  was  fitted  by  the  proposed 
models  (Figure  2) .  The  resulting  strain  distributions  are  plotted  together 
in  Figure  3.  As  can  be  seen,  elastic  strains  decrease  with  the  inverse 
of  the  square  of  the  distance  from  the  axis  of  the  borehole;  the  effect  of  the 
non-linearity  of  the  stress-strain  law  is  essentially  felt  as  a  faster  than 
elastic  rate  of  decrease  of  the  shear  strain.  When  the  initial  modulus 
coincides,  all  stress- strain  laws  predict  identical  strains  at  large  distances. 


Once  the  strain  distribution  is  known,  the  torque  required  and  the 
corresponding  displacement  field  can  be  easily  obtained. 


\  \\\ 


BOREHOLE  SHEAR  STRESS  ^=0.65  kg/cm 
BOREHOLE  RADIUS  rD  =O.IOm 


DISTANCE  FROM  THE  AXIS(m) 


-  LINEAR  ELASTICONITIAL  MODULUS) 


HYPERBOLIC 
POWER  LAW 


HYPERBOLIC  TANGENT 
RAMBERG- OSGOOD  TYPE 


LINEAR  ELASTIC  (SECANT  MOOULUS) 


Fiquro  Comparison  of  the  corrosron-ii !«  strain 
distributions 


The  torque  per  unit  length  ir»: 


Te 


2 


2 

ir  r  t 


o  o 


(8) 


The  displacement  field  is  derived  from  the  relation: 


-1H.  1L 

■'•r  r 


(9) 


where  u  is  the  ‘  -  component  of  the  displacement  (all  other  components  vanish). 


(lO) 


For  example,  for  the  elastic  case: 

2 
o 
2 

<• 

This  equation  can  bo  readily  integrated-  The  partial  derivative 
becomes  an  ordinary  one,  as  u  only  depends  on  r.  Assuming  that  the  displacement 


uu 

3r 


u 

r 


t  r 
o  o 


vanishes  at  infinity,  the  solution  is  (see  Appendix  IV  for  the  derivations) 

2 


u  =  - 


t  .r 
o  o 

2yor 


(11) 


The  above  result  can  also  be  expressed  as  a  function  of  the  displace¬ 
ment  at  the  interface  of  »  or  the  angle  rotated  Gq.  Table  2  provides 
those  forms  of  the  solution. 


Table  2.  Summary  of  elastic  results 


AS  A 

FUNCTION  OF 


VARIABLE 


2p  u  r 
o  o  o 


2u  6  r2 

COO 


2u  r 
o  o 
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-4xu  u  r 
o  o  o 


.  2 
-4rue  r 
o  o 


h  H 

£  B 


The  above  formulae,  in  spite  of  their  simplicity  end  the  number  of 
assumptions  involved,  are  of  notable  use  in  later  considerations. 


Derivations  similar  to  the  elastic  ones  can  be  carried  out  for 
other  constitutive  relations.  The  hyperbolic  model  will  be  studied  here 
because  of  its  later  utilization  for  interpretation  of  field  data.  The 
displacement  distributions  fulfills  now  the  equation: 


_  t  r 
2  o  o 


again  subject  to  the  boundary  condition: 


iim  u  =  o 


xyj*  V>*, .-■«>«  A*~ 


a2u 


3t 


2 

3  u. 


3t 


■=  (A+2g) 


3A 

3r 


_  3w 
2g  _ z 

r  30 


3  w 

A+2p  3A  „  z 

~  re  +  2y  aT 


(15) 


where : 


3ur  Ur  1  3u0 

A  3r  +  r  +  r  36 


1  ( 

w  =  — —  ] 

z  Zr  l 


3u, 


r  3r  +  U8  30 
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and  A,p  are  complex  Lame  constants 


In  the  present  case  u  =Oand  •—  =  0;  then,  the  first  of  equations 

IT  du 

(15)  becomes  an  identity  while  the  second  one  takes  the  form: 


3V 


3t 


=  v 


3?'u, 


3r 


i  !!§ 

r  3r 


(16) 


where  v  ~  =  —  is  the  square  of  the  complex  velocity  of  shear-wave  propagation, 
s  p 


The  boundary  conditions  for  equation  (16)  are  vanishing  displace¬ 
ments  at  infinity  and  harmonic  input  at  the  cylinder-soil  interface,  whether 
in  the  form  of  controlled  displacements  or  controlled  stresses.  Again  represent¬ 
ing  by  u  the  only  non-vanishing  component  of  the  displacement,  those  boundary 
conditions  can  be  expressed  mathematically  as: 
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lim  u  =  o 
r  -*■  « 


and 


=  u  e 
o 


iwt 


for  r  =  r 

o 


(controlled  displacements) 


or 


T 


.ii.it 

t  e  for  r  =  r 

o  o 


(controlled  stress) 


Since  only  the  steady-state  solutions  are  sought,  initial  conditions 
are  not  required. 

Equation  (16)  can  be  integrated  the  standard  method  of  separation 

of  variables,  that  is,  finding  the  solutions  of  the  form: 

u(r,t)  =  R(r)  T(t) 

where  R  depends  only  or  r  and  T  only  on  t  (the  bar  is  used  simply  to 
differentiate  the  symbol  from  that  used  for  torque) 

Substituting  into  equation  (16) : 


ipll 

T 


R" 


-1- 

r  R 


(17) 


Since  the  left  hand  side  of  equation  (17)  is  only  a  function  of 

time  and  the  right  hand  side  depends  only  on  the  distance,  both  sides  are 

2 

equal  to  a  constant,  henceforth  named  -»  •  Equation  (17)  then  generates: 
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(IB) 


'|E!1  +  I  51 

R  +  r  R 


All  solutions  of  equation  (18)  are  of  the  form: 


T  =  T  e 
o 


where  T  is  a  constant, 
o 


Equation  (19)  can  be  regrouped;  by  changing  variables  to  the 

dimensionless  distance—  the  equation  takes  the  form  of  a  Bessel  differential 

s 

equation : 


d2R 

/  \2  + 


© 


f)  -f\ 

. s/  dj'^-j  [Vs/ 


1  =  0 


the  general  solution  of  which  is: 


(?)  ■  ■-,,(») . 


where  J,  is  a  Bessel  function  of  the  first  kind 
Y.  is  a  Bessel  function  of  the  second  kind 


A1  and  B’  are  arbitrary  complex  constants 
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( I ) 

It  is  then  readily  apparent  that  H  corresponds  to  waves 

(2' 

converging  towards  the  axis  r=o  while  Hj  '  represents  waves  diverging 
from  the  axis.  It  can  further  be  noticed  that  real  wave  numbers 
correspond  to  travelling  waves  while  imaginary  wave  numbers  represent 
standing  waves.  In  this  situation,  for  displacement  controlled  boundary 
conditions,  the  general  solution  for  the  displacement  becomes: 


u 


r  e 
o  o 


(24) 


where  <»  is  the  input  frequency. 


Simple  differentiation  of  equation  (24)  provides  the  strain  field: 


e  = 


and  similarly: 


T 


T 


9. 


(25) 


(26) 


(27) 
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Although  it  is  not  shown  here,  it  is  immediate  to  notice  that 
the  static  solutions  are  easily  recovered  from  the  above  formulae  by 
taking  limits  when  rQ  ->  “  (substitute  the  Bessel  functions  by  their 
assymptotic  approximations) . 


Since  the  material  is  linear,  the  above  equations  (24)  and  (27) 
can  be  used  to  obtain  the  results  for  the  stress  controlled  boundary. 

For  r  =  r  ,  equation  (22)  becomes: 


h(2)/^\ 

Mro0ow  2  \Vs  /  i 


which  substituted  into  equations  (24)  to  (27)  yields: 


t  =  T 


,v  <2,(~) 

_p  _s  _ \_s_/  i 

.  ■;”(?) 

£ _ \  s  /  iw'; 


H22>  “ 

\s  /  irnt 


■!"(?) 


T  =.  2*r  2j  ciu,t 
Z  o  o 


hhM' 


3.3  FINITE  CYLINDER  -  STATIC  AND  DYNAMIC 


For  an  elastic  half-space  bounded  by  a  free  or  fixed  surface, 

the  problem  of  a  finite  cylinder  rotating  in  a  borehole  normal  to  the 
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surface  has  been  partially  solved  by  Boyer  and  Oien  (1972)  -The  full 
fields  of  each  variable  were  not  derived  but  torque-angle  ratios  were 
provided  as  a  function  of  the  parameters  of  interest.  The  problem  will 
be  presented  here  following  the  reference  mentioned.  The  approach  used 
there  is  somewhat  similar  to  that  used  in  the  previous  section,  but  the 
mathematical  treatment  becomes  rather  involved  and  the  reader  is  referred 
to  the  original  paper  for  more  specific  details. 

The  equation  to  be  solved  is: 


1  1_  /  „  3u  \  u  A  a2u  j  ,  2 

r  Or  y  j  2  '  2  'C  u  = 


.i  -  —  ,  u-  being  the  input  frequency,  n  the  height  of  the  cylinder 
s 

vg  the  shear  wave  velocity, 
u  is  the  0  component  of  the  displacement 

r,  G,  z  are  the  cylindrical  coordinates  and  all  distances  are 
scaled  by  division  by  h/2. 


The  boundary  conditions,  placing  the  origin  of  coordinates  at  the 
center  of  gravity  of  the  cylinder,  take  the  form: 


u  =RQ  at  r  =  R  !z(  "  h 


'  -  =  0 
rG 


at  r  =  R  |z|  >  ’> 


u  =  0 


at  r  *  “ 
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INFLUENCE  FUNCTION 


Fiqure  7 


t'iquru 


Ue.i  I  part.  of  t  hi'  influence  function  vs.  frequency 
for  a  cylinder  with  qoomotry  h/r  -=  2.0  at  various 
distances  from  either  a  free  or  °  fixed  surface. 
(After  Boyer  and  Oi.en,  1072) 


0  12  3  4  5 
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Norma  1 1 sod  imaqinery  part  of  the  influence  function 
vS ;  frequency  for  a  cylinder  with  qeometry  h/r  -  2 
at  various  distances  from  either  a  free  or 
fixed  surface.  (After  Boyer  and  Oirn,  1972) 


4.  numerical  analysis 


It  was  pointed  out  at  the  onset  of  Chapter  3  that  close-form 
solutions,  in  spite  of  their  usefulness,  required  for  their  derivation 
the  incorporation  of  a  number  of  assumptions.  Numerical  analysis  has  the 
disadvantage  of  providing  answers  only  at  discrete  points  of  the  space 
of  input  parameters;  however,  it  allows  solving  the  problem  without 
unnecessary  assumptions  for  each  specific  set  of  input  data. 

The  purpose  of  the  numerical  analysis  reported  in  the  prasent 
chapter  is  to  assess  the  magnitude  of  the  error^  introduced  in  the  analytical 
treatment  and  to  correct  them  when  deemed  necessary.  Those  activities 
are  orientated  towards  the  final  ooal  of  providing  a  realistic  method  of 
obtaining  the  value  of  the  secant  shear  modulus  of  ;'ne  soil  as  a  function 
of  strain  level,  based  on  field  measurements  of  torque  and  rotation  during 
the  operation  of  the  BSD. 

4.1  MEl  IDOLOGY 

The  field  operation  of  the  BSD  has  two  characteristics  which  are 
of  decisive  importance  in  the  selection  of  the  numerical  approach:  those  are 
the  facts  that  the  actuator  is  cycled  at  a  single  frequency  and  under  steady- 
state  conditions.  In  those  circumstances,  a  finite  element  analysis  working 
in  the  frequency  domain  shows  its  best  advantages.  Non-linearities  can 
then  be  handled  via  an  equivalent-linear  method,  where  the  solution  procedure 
is  iterated  until  the  defr  nation  levels  achieved  are  consistent  with  the 
values  of  the  modulus  used. 
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Based  on  t'-.e  above  considerations,  the  program  FILEL,  developed 
by  Mr.  D.  Hitchings  at  Imperial  College,  London,  was  selected  for  carrying 


out  the  computations.  General  information  on  FINEL,  including  its  methodology 
an<I  capabilities,  can  be  found  in  Appendix  i>. 

The  version  of  FINEL  used  in  the  present  project  assumes  axial 
symmetry  of  geometry  and  input  conditions.  In  order  to  adapt  it  best  to 
the  particular  configuration  of  the  problem  under  consideration,  several 
modifications  were  introduced  in  the  computer  program: 

The  equivalent- linear  procedure  was  implemented.  Aj  mentioned, 
this  allows  the  use  of  a  shear  modulus  "consistent"  with  the 
strain  level  produced  at  each  point.  This  consistency  will 
be  discussed  in  further  sections. 

Special  non-reflecting  boundaries  were  introduce’  to  allow  the 
use  of  a  small  finite-element  mesh  while  appropriately  modelling 
the  radiation  of  energy  out  of  the  finite-element  region.  The 
cylindrical  surface  of  the  region  was  provided  with  a  boundary 
similar  to  that  proposed  by  Lysaer  and  Waas  (1072)  but  modified 
for  axisymmetric  as  opposed  to  plane-strain  conditions;  such 
boundary  is  theoretically  perfect  as  long  as  the  soil  outside 
the  boundary  behaves  linearly  and  is  horizontally  layered. 

The  top  boundary  was  equipped  with  an  optional  non-reflecting 
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boundary  of  the  type  described  by  Lysmer  and  Kuhlemeyer  (1969)  ; 

energy  is  absorbed  at  the  boundary  by  means  of  suitably  selected 


viscous  dasnpots. 
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Finally,  the  nodes  at  the  borehole  wall  in  contact  with  the 

BSD  are  subjected  to  the  boundary  condition: 

iwt 

u  =  u  e 
o 

while  all  other  nodes  in  the  borehole  walls  are  considered  free  to  move. 

It  is  interesting  to  note  here  that,  although  the  strains  considered 

are  large  enough  to  introduce  constitutive  non-linearities,  geometrical 

non-linearities  are  not  taken  into  account  in  the  analysis  due  to  the 

relatively  small  rotations  involved.  Consequently,  the  only  non -vanishing 

displacement  is  u.  and  the  only  non-zero  strains  are  e  .  and  e  these 
o  rt)  z  0 

strains  are  related  to  the  displacement  field  by  the  small  displacement 
formula cf on. 

Finally,  because  of  the  small  size  of  the  domain  investigated  in 
comparison  with  the  wave  lengths  induced  (small  dimensionless  input  frequency) , 
the  effect  of  material  damping  becomes  negligible  with  respect  to  that  of 
radiation  damping.  It  has  therefore  been  neglected  in  the  numex'ical  calculations. 

4.2  DUPLICATION  OF  ELASTIC  RESULTS 

The  purpose  of  duplicating  e  results  already  known  from  close- 
form  solutions  is  essentially  that  of  verifying  the  adequacy  of  the  methodology 
selected  for  the  numerical  analysis.  Of  particular  interest  in  this  respect 
is  the  validation  of  the  non-ref lectir.g  boundaries  w.iich  were  added  to  FINEL 
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for  the  present  study.  Apart  from  routine  checks  of  little  interest  here, 
two  basic  tests  will  be  included  in  the  present  section  because  of  the 
insight  that  they  provide  into  the  general  treatment  of  the  problem  and 
because  they  can  be  used  as  proof  of  the  validity  of  *-he  answers  generated 
by  the  computer  jr  roc  ram. 

The  two  basic  tests  mentioned  are  described  in  detail  in  the  next 
subsections.  They  include: 

a)  Infinitely  long  cylinder  rotating  harmonically  in  an  elastic 
medium.  All  variables  are  known  from  analytical  solutions. 
This  test  is  designed  for  general  verification  of  the  computer 
program  and,  in  particular,  of  the  axisymmetric  non-reflecting 
boundary.  Although  the  modelling  is  two-dimensional,  the 
problem  is  only  one-dimensional. 

b)  Finite  cylinder  rotating  harmonically  in  an  elastic  medium. 

The  torque-angle  relationships  are  known  from  analytical 
solutions.  This  test  is  designed  for  validation  of  the 
program  under  a  truly  two-dimensional  situation  and,  in 
particular,  for  verification  of  the  plane  non-reflecting 
boundary.  It  will  be  seen  later  that,  based  on  this 
example,  upper  bounds  can  be  found  to  the  effects  of 
layering  in  distorting  the  field  measurements. 
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4.2.1  Infinitely  Long  Cylinder 

The  following  input  parameters  were  used  for  the  case  of  an  infinitely 
long  cylinder  rotating  harmonically  in  an  elastic  medium: 


Soil  density 

P 

=  2000  Kg/m3 

Shear  modulus 

P 

=  1.8  x  108  N/m2 

Borehole  radius 

r 

o 

=  0.1  m 

Input  frequency 

0) 

=  100  rad/sec 

Angular  displacement 

u 

-4 

=  10  m  (half-amplitude) 

The  above  conditions  are  representative  of  a  typical  field  situation, 
assuming  that  the  soil  behaves  elastically.  Shear-wave  propagation  velocity  is 
300  m/sec,  based  on  the  values  of  u  and  p  selected. 

Analytical  solutions  were  provided  for  those  conditions  in  Chapter  3. 
Figures  lo  and  11  present  the  comparison  of  analytical  and  numerical  solutions 
for  the  real  and  imaginary  parts  of  the  displacement  field,  respectively.  Figures 
12  and  13  depict  similar  results  for  the  strain  field.  Note  that,  under  the 
conditions  examined,  the  imaginary  part  of  the  response  is  very  small  while  the 
real  part  is  practically  equal  to  the  static  solution;  this  is  true  for  both 
displacements  and  strain  fields. 

As  can  be  seen  in  the  referenced  figures,  the  agreement  of  the  finite- 
element  calculations  and  the  close-form  solutions  is  extremely  good.  In  order  to 
give  an  idea  of  the  importance  of  the  cylindrical  non-reflecting  boundary,  it 
should  be  mentioned  here  that  errors  in  excess  of  100%  were  not  uncommon  in  the 
neighborhood  of  the  boundary,  when  fixed  or  free  conditions  were  used.  A  notable 
improvement  over  free  of  fixed  boundary  conditions  could  be  obtained  by  using  the 


w  =  100  rod /tec 
=  300  m/sec. 

©  =  F.E.  SOLUTION 


"10 -I 
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Figure-  11.  Comparison  of  analytical  and  minor i cal 

solutions  -  imaginary  pari  of  t ho  disrd  acoment 


a)  -  100  rod/tec 
:  300  m/sec 
®  =  F.E.  SOLUTION 

-  :  EXACT  SOLUTION 

- =  STATIC  CASE 


Comparison  of  analytical  and  numerical 
solutions  -  real  part  of  the  strain 


plane-strain  Lysmer  -  Waas  boundary;  this  condition  is  perfect  as  the  radius 
grows  to  infinity  but  still  generates  substantial  errors  at  the  relatively 
small  radii  of  the  geometries  of  interest.  It  is  due  to  those  circumstances 
that  a  decision  was  made  to  develop  and  implement  the  axisymmetric  boundary. 

4.2.2  Finite  Cylinder 

The  case  of  a  harmonically  rotating  cylinder  of  finite  length  was 
studied  under  the  same  conditions  and  parameters  used  in  the  analysis  of  the 
infinite  cylinder.  A  number  of  tests  were  carried  out  in  order  to  ascertain 
the  validity  of  the  results  of  the  program  and  the  adequacy  of  the  plane  non¬ 
reflecting  boundary;  this  was  done  by  comparison  with  the  close  form  solutions 
obtained  in  the  previous  chapter. 


Of  special  interest  here  is  a  series  of  tests  performed  on  the  mesh 
shown  in  Figure  14  .  The  cylinder  is  0.5  m  long  and  the  distance  from  its 
centre  to  the  top  boundary  is  again  0.5  m.  Three  alternative  conditions  were 
used  at  the  top  boundary:  free,  simulating  the  free  surface;  fixed,  modelling 
an  infinitely  rigid  layer  and  non- ref lecting,  to  represent  infinite  embedment. 
Table  3  presents  the  comparison  of  the  torque  -  rotation  ratios  obtained  by  the 
program  and  from  the  close-form  solutions.  Values  of  the  influence  function 
F  defined  in  Figure  5  are  used  for  this  com,  arisc.n. 

Table  3  Comparison  of  the  real  part  of  the  influence  function  F  for 
d/h  =  l 


!  F.E.  SOLUTION 


5 

r 

o 


ANALYTICAL  SOLUTION 
(after  Boyer  &  Oien,  1972) 


5 

r 

o 


*=  2 
r 

o 


FIXED  SURFACE 
FREE  SURFACE 
INFINITE  EMBEDMENT 


2.35 

2.31 


2.32 


2. 36±0.05 


2.90+0.05 
2. 78±0.05 
2.85+0.05 


v^-~ ***y * +r*-~ ■- ■¥■»»«-».  / 


Notice  that  the  analytical  solution  for  fixed  and  free  boundaries 
cannot  be  expected  to  coincide  with  the  numerical  results  since  both 
correspond  to  different  cylinder  geometries;  besides,  the  numerical  results 
include  two  rather  than  one  fixed  or  free  planes  due  to  the  symmetry 
conditions  imposed  (unlike  the  analytical  values) .  The  results  in  the 
third  column  correspond  to  the  geometry  which  is  closest  to  the  present 
BSD  design  within  the  set  of  analytical  results  available.  The  tolerances 

shown  in  the  analytical  values  represent  the  possible  errors  in  reading 
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from  the  figures  of  Boyer  and  Oien  (1972)  . 


The  strain  fields  developed  in  the  above  cases  are  shown  in  Figures 
lb  to  20,  It  can  be  noticed  that  the  higher  strains  are  barely  affected 
by  the  presence  of  a  fixed  or  free  boundary  located  half  a  cylinder  length 
away  from  the  end  of  the  cylinder.  This  conclusion  can  also  be  drawn  from 
the  torque-rotation  ratios.  As  expected,  the  results  for  infinite  embedment 
lie  in  between  the  cases  of  fixed  and  free  boundary  conditions  on  the  top 
plane.  It  cam  therefore  be  stated  that  although  the  plane  non- reflecting 
boundairy  is  working  adequately,  the  importance  of  its  implementation  on  a 
satisfactory  m  deiling  of  this  problem  is  considerably  less  than  in  the  case 
of  the  cylindrical  boundary. 


One  more  point  worth  noticing  is  that  e0z  becomes  infinity  at  the 
extreme  of  the  cylinder.  However,  for  the  geometries  of  interest,  its 
contribution  to  the  required  torque  is  small.  Its  importance  would  increase 
as  the  cylinder  becomes  very  short  compared  to  its  diameter. 


Figure  17.  6r©  strain  distribution 

t  -.*p  sus  face  non-ref iecting 


Figure  IS.  GOZ  strain  distribution 
top  surface  non-reflecting 


53 


V 


**V 


4.3  NON-LINEAR  ANAL x SIS 


The  tests  reported  under  the  present  section  are  orientated  to  two 

purposes : 

To  verfiy  the  implementation  of  the  "equivalent-linear"  method 
and  its  ability  to  represent  the  constitutive  non-linearities 
of  the  problem. 

To  study  the  influence  of  the.  non-linearity  of  the  problem  on 
the  ratio  of  torques  required  in  the  cases  of  infinite  and 
finite  cylinder. 


The  object  of  the  two  goals  mentioned,  particularly  the  second  one, 
will  become  more  clear  in  Chapter  6.  There,  a  procedure  will  be  proposed 
to  construct  an  arbitrary  u-e  (shear  modulus  -  strain  level)  curve  for  the 
soil,  based  on  the  r -c>  (torque  -  rotation)  measurements  obtained  in  the 
field.  The  present  section  is,  in  a  way,  the  justification  of  the  assumptions 
embodied  in  that  procedure. 

Two  points  must,  however,  be  made  here  concerning  the  "equivalent 
linear"  method.  The  first  one  is  that  the  method  originates  from  a  one¬ 
dimensional  philosophy  and  its  generalization  to  two  dimensions  is  by  no 
means  unique.  The  modulus  for  each  shear  strain  may  be  related  to  the 
corresponding  strain,  to  a  specific  one  of  the  two  strains,  to  the  modulus 
of  both  or  to  any  other  shear  strain  measure.  Each  hypothesis  has  advantages 
and  disadvantages  but  the  present  state  of  the  art  in  dynamic  soil  modelling 
does  not  provide  a  definite  argument  in  favour  of  any  of  them.  Under  these 


conditions,  and  given  the  fact  that  e  is  larger  than  almost  everywhere 


and  the  most  important  one  to  reproduce  correctly,  we  have  made  the  somewhat 
arbitrary  decision  of  assuming  that,  at  each  point  within  the  soil  and 
each  moment  in  time,  there  is  a  single  value  of  the  secant  modulus  relating 
shear  strains  and  the  corresponding  shear  stresses,  furthermore,  the 

deterioration  of  such  modulus  is  totally  controlled  by  alone. 

/ 

s' 

The  second  point  is  that  the  "equivalent  linear"  procedure  was 
developed  and  has  been  primarily  applied  in  the  field  of  earthquake 
engineering?  there,  the  dynamic  and  inertial  aspects  of  true  wave  propagation 
prevail  over  the  static  or  quasi-static  distributions  of  stress  which  are 
only  controlled  through  stiffness  characteristics.  Furthermore,  the 
behavior  of  the  soil  consists  mainly  of  transient  responses  to  trains  of 
irregular  cycles.  ”  consequence,  in  earthquake  engineering,  the  "consistency" 
of  the  strain  level  and  the  secant  modulus  is  customarily  satisfied  by  taking 
for  the  latter  the  modulus  corresponding  to  65%  (approximately)  of  the  maximum 
strain  developed  at  each  point  of  the  soil.  The  above  conditions  do  not 
apply  exactly  to  the  problem  being  examined  here.  Hence,  some  research 
was  done  in  optimizing  the  strain  level  at  which  the  secant  modulus  is  taken 
as  "consistent".  It  is  clear  that,  under  static  conditions,  100%  of  the 
strain  would  be  the  optimum  value.  The  conclusions  of  our  investigation 
were  that,  although  convergence  was  somewhat  hindered,  values  close  to  100% 
of  the  strain  provided  the  best  answers.  As  a  consequence,  a  value  of  100% 
was  used  in  all  the  runs,  including  those  in  the  present  section. 

One  more  aspect  should  be  mentioned  before  going  into  the  tests. 
Since  a  certain  constitutive  law  has  to  be  assumed  in  order  to  be  able  to 
obtain  some  quantitative  results,  a  hyperbolic  stress-strain  law  has  been 
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postulated  on  the  soil  for  the  purpose  of  the  numerical  tests  included 
under  the  present  section.  It  is  reasonable  to  assume  that  the  verifications 
performed  for  a  hyperbolic  law  would  in  fact  hold  if  other  non-linear 
stress-strain  laws  had  been  used.  Although,  strictly  speaking,  since  no 
single  law  is  universally  applicable,  one  would  need  to  repeat  the  process 
with  all  other  possible  laws,  this  is  not  actually  required  since  the  final 
recommended  procedure  is  based  on  the  hyperbolic  law.  However,  it  is  important 
to  stress  that,  in  spite  of  the  fact  that  the  procedure  will  be  based  on 
the  hyperbolic  law,  it  will  not  be  assumed  that  the  soil  behaves  according 
to  such  law;  rather,  as  will  be  seen  in  Chapter  6,  the  constitutive  law  will 
be  reproduced  by  a  series  of  hyperbolic  approximations. 

4.3.1  Infinitely  Long  Cylinder 

The  primary  purpose  of  these  tests  is  to  check  the  "equivalent- 
linear"  procedure.  We  have  neither  been  able  to  locate  nor  deemed  sensible 
to  try  and  produce  any  non-linear  close-form  solutions  for  the  present 
boundary  value  problem  under  general  dynamic  conditions.  However,  in  view 
of  the  fact  that,  at  the  range  of  frequencies  of  interest,  the  static  solution 
closely  approximates  the  real  part  of  the  dynamic  solution,  comparisons  have 
been  established  between  analytical  solutions  of  the  former  and  the  finite-element 
predictions  of  the  latter. 

In  the  next  two  tests,  the  soil  constitutive  law  is  hence  hyperbolic 
and  can  be  described  by  means  of  the  }i-e  curve  shown  in  Figure  21  .  The 


input  parameters  common  to  both  tests  are: 


Soil  density 


P 


Initial  shear  modulus  y  = 

o 

Limiting  shear  strength  x  = 

m 

Borehole  radius  r  = 

o 

Two  different  angular  displacements  were  given  at  the  borehole  wall 

-4  -4 

in  the  first  case,  u Q  =  10  m  and  in  the  second,  uq  =  3.5  x  10  m.  Two 

different  input  frequencies  were  also  used:  a>  =  10  rad/s  and  w  =  100  rad/s 

The  comparisons  between  the  analytical  static  solution  and  the  real  part  of 

the  dynamic  numerical  solution  are  presented  for  strains  (Figrre  22  and  23) 

and  displacements  (Figure  24  and  25) .  The  required  torques  are  shown  in 

Table  4. 


2  x  10  Kg/m3 
1.8  xlO8  N/m2 
5  x  105  N/m2 


0.1  m 


Table  4.  Comparison  of  torques  in  hyperbolic  soil 


INPUT  ANGULAR 
DISPLACEMENT, U 
(m)  ° 

REQUIRED  TORQUE 

,  T(Nm/m) 

CLOSE-FORM 

STATIC 

F.E.  REAL  PART  OF  DYNAMIC 

w=lOrad/s 

w=100rad/s 

-i  ~  -4 

4 

4 

4 

1.0  x  10 

1.61  x  10 

1.64  x  10 

1.64  x  10 

-4 

4 

4 

4 

3.5  x  10 

2.89  x  10 

3.15  x  10 

3.15  x  lO 

4.3.2  Finite  Cylinder 

The  above  tests  were  repeated  under  identical  conditions  with  a 
cylinder  of  finite  length  (0.5  m)  located  at  large  depth  in  a  borehole.  The 
finite-element  mesh  "ced  in  the  calculations  is  shown  in  Figure  26.  Again, 


STRAIN  (- 


Fiqure  23.  Comparison  between  analytical  ana  numerical 
strain  values 
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HALF  LENGTH 


0,75m 


the  real  component  of  both  strain  distributions  (e^g  and  e^)  has  been 
plotted  by  the  computer  (Figures  27  and  30) .  As  it  was  mentioned 
earlier  in  this  report,  it  can  be  easily  observed  that  the  greater  the 
non-linearity  of  the  part  of  the  stress-strain  law  involved,  the  faster 
the  decay  of  the  strains  with  distance  from  the  cylinder. 


The  values  of  torque  required  to  induce  the  imposed  boundary 
displacements  cannot  be  verified  against  analytical  solutions  and  are 


given  in  Table  5. 


Table  5.  Torque  values  for  finite  cylinder  in  hyperbolic  soil 


INPUT  ANGULAR 


TORQUE  REQUIRED 


DISPLACEMENT,  u  TO  PRODUCE  IT,  T 

(m)  °  (N  m) 

«  » 

1.0  x  10~4  !  9.20  x  lO  3 

j  3.5  x  10  ‘i  1.79  x  10  1  , 

Quite  naturally,  the  torque  required  divided  by  the  displacement 

produced  decreases  as  the  soil  becomes  more  non-linear  due  to  the 

degradation  of  the  secant  stiffness. 


4.3.3  Torque  Ratios 


It  was  mentioned  at  the  onset  of  Section  4.3  that  one  of  the 


purposes  of  the  non-linear  analysis  was  to  determine  the  influence  of  the 
non-linearity  on  the  ratio  of  torques  for  the  case  of  long  and  short 
cylinders.  This  point  constitutes  the  object  of  the  present  section. 


Figure 


Gre  strain  distribution  for  a  non  1 inea: 
-  *itutivc  law  (U  =  3.5slO-4a5 


The  torque  T  required  to  produce  a  steady-state  rotation  0 

on  a  finite  cylinder  of  length  h  can  be  expressed  in  terms  of  the  torque 

T  required  to  induce  the  same  effect  in  an  infinitely  long  cylinder. 

L 

We  will  concentrate  on  the  real  part  of  the  torque  since  it  was  already 
shown  that  the  imaginary  part  is  very  small.  The  expression  will  be  of 
the  type: 


T  =  T  F 
F  L 


From  dimensional  analysis,  in  the  elastic  case  F  (the  subscript 

£ 

E  stands  for  elastic  is  only  a  function  of: 


A  quick  look  to  Figures  5  to  8  shows  that  the  influence  of 

d  can  be  neglected  when  <•*.  >1;  since  this  is  always  true,  F  can  always  be 
h  h 

taken  as  its  value  when  d  -**>.  Similarly,  due  to  the  facts  that  a)  Figures 

h  wr  urQ 

5  and  7  start  the  curves  horizontally  and  b)  - is  very  small,  - - 

s  s 

can  be  considered  zero  for  the  purpose  of  determing  F.  Hence: 


F  =  F 
E  E 


(0 


where  ^  =  5  in  the  present  design  of  the  BSD. 


In  the  case  of  non-linear  soil,  the  small  influence  of  frequency 

d 

was  already  seen  in  Section  4.3.1.  As  for  the  influence  of  h,  it  is  clear 
that,  although  it  cannot  be  directly  quantified,  it  must  be  even  smaller 


than  in  the  elastic  case  due  to  the  fact,  that  the  strains  decrease  faster 
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with  distance  from  the  cylinder.  One  more  parameter  has  to  be  considered 
in  the  case  of  hyperbolic  soil;  several  equivalent  ones  could  be  used,  for 
example,  e  (the  value  of  the  maximum  e  .  strain  developed) .  Then,  in  the 

0  To 

case  of  hyperbolic  soil: 


F 


H 


(4) 


The  value  of  h__  is  still  5. 
r 

rotation  level;  clearly. 


The  value  of  e  increases  with  the 
o 


Fy  =  Fg  when  eQ  <<  1  (5) 

since  the  soil  behaves  elastically  for  sufficiently  small  strains  (Section 
3.1b).  The  question  is  how  does  F  vary  with  e  as  the  latter  increases. 

n  0 

The  results  obtained  in  the  tests  are  shown  in  Table  6. 

Table  6.  FH  factor  for  hyperbolic  soil 


MAXIMUM  e  „  STRAIN 
r0 

DEVELOPED,  CQ 
(-) 

F  FACTOR 

H 

(-) 

SMALL 

2.31 

-3 

2.93  x  10 

2.25 

-2 

3.17  x  10 

2.27 

As  can  be  seen  in  the  above  Table,  the  factor  F  remains  essentially 

H 

constant  over  the  range  of  strain  levels  of  interest;  the  differences  are 
so  small  that  they  are  almost  on  the  order  of  the  accuracy  of  the  results 


presented. 


4.4  INFLUENCE  OF  LAYERING 


To  this  point,  all  the  analysis  has  been  made  on  the  assumption  that 
the  soil  is  a  homogenous  medium.  In  all  probability,  layers  will  exist  in 
the  real  soil  profiles  where  the  BSD  is  to  be  used.  The  present  section 
aims  to  define  the  influence  of  layering  on  the  torque  and  rotation  measure¬ 
ments  on  which  the  soil  modulus  calculations  will  be  based. 


It  has  been  shown  that  the  influence  of  a  fixed  or  free  boundary 
on  the  torque  required  to  produce  a  certain  rotation  of  the  cylinder  is 
greater  in  the  linear  than  in  the  non-linear  case.  In  fact,  it  is  easy 
to  see  that,  as  the  strains  increase,  t  around  the  cylinder  approaches 

TV 

the  limiting  strength  of  the  material,  irrespective  of  any  free  or  fixed 
surfaces,  no  matter  how  close  to  the  cylinder.  Therefore,  the  elastic 
case,  eind  hence  the  small  strain  operations,  are  the  most  sensitive  to  the 
proximity  of  very  hard  or  very  soft-  layers  above  or  below  the  BSD. 


Fortunately,  even  in  the  elastic  case,  there  does  not  seem  to 

be  reason  for  concern.  The  worst  condition  will  arise  if  an  infinitely 

rigid  or  infinitely  soft  layer  are  at  a  short  distance  from  the  BSD. 

But  even  for  this  worst  condition  of  the  worst  case.  Table  3  clearly 

shows  that,  as  long  as  the  extreme  of  the  BSD  is  more  than  25cm  away 

from  the  infinitely  rigid  or  soft  layer,  the  influence  in  the  torque 

is  less  than  1%  (half  of  what  is  shown  in  Table  3  which  includes  two 

symmetric  fixed  or  free  surfaces  rather  than  one) .  Even  when  the  distance 

is  on  the  order  of  lOci., ,  the  torque  can  only  be  seen  to  vary  on  the 

order  of  10%  in  the  far  more  disadvantageous  configuration  of  Figure  7 

(h/r  =  2,  rather  than  5). 
o 


msis 


!jsj| 
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vertical  direction  of  the  t :  strains,  which  are  responsible  for  transmitting 

the  BSD  project  to  levels  higher  or  lower  than  the  location  of  the  BSD. 


In  summary,  since  a)  the  geotechnical  profile  will  be  known, 
b)  the  operation  will  concentrate  on  relatively  large  strains  and  c)  there 
cure  no  infinitely  rigid  or  soft  layers  in  nature,  it  is  not  likely  that 
errors  as  high  as  2%  will  be  introduced  in  the  measurements  as  a  consequence 
of  soil  layering. 


4.5  EFFECT  OF  DISTURBANCE 


It  is  well  known  that  any  method  of  creating  an  opening  in  the 
ground,  no  matter  how  careful,  will  inevitably  produce  a  zone  of  disturbance 
around  the  opening.  Such  an  effect  ha-  been  considered  at  all  stages  of 
the  BSD  project,  although  not  always  with  identical  priority.  At  the 
early  stages  of  the  project,  when  the  device  was  expected  to  propagate 
waves  from  the  BSD  to  a  pick-up,  the  importance  of  a  narrow  zone  of 
disturbance  was  very  limited.  However,  as  preliminary  analyses  showed  the 
impossibility  of  that  operation  and  the  conceptual  design  shifted  to  a 
single-point  measurement  configuration,  the  minimization  of  borehole  disturbance 
became  paramount.  It  was  at  that  time  that  a  decision  was  made  to  incorporate 
a  self-boring  tool  at  the  front  of  the  BSD. 

Figure  31  presents  two  ideal  curves  as  typically  originated  by  a 
standard  pressuremeter  and  a  self-boring  pressuremeter.  Standard  pressuremeters 
completely  relieve  the  horizontal  stress;  the  initial  modulus  obtained 
is  low  compared  to  the  undisturbed  one,  due  to  both  the  stress  relief  and 
the  disturbance  caused.  It  is  usually  considered  that  the  undisturbed 
situation  has  been  reached  where  the  o  -  u^  curve  becomes  a  straight  line 
as  predicted  by  elastic  considerations. 

A  self-boring  pressuremeter  however  does  not  produce  a  total 

rel  ef  of  the  horizontal  normal  stress  and  causes  much  less  disturbance 
30 

(Hughes  et  al,  1977);  also,  the  deterioration  as  a  fixed  rather  than  a 
radially  smaller.  If  the  BSD  had  been  designed  as  a  fixed  rather  than  a 
radially  expanding  surface,  it  would  have  been  appropriate  to  perform  an 
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advised. 


5.0  SYSTEMS  DESIGN  STUDY 


5.1  INTRODUCTION 


The  purpose  of  this  part  of  the  project  was  to  examine  the  practical 


aspects  of  the  complete  field  measurement  system.  The  scope  of  the  study 


ranged  from  problems  such  as  cost  and  complexity  in  manufacturing  processes, 


to  those  related  to  handling  and  operation  of  the  equipment  in  the  field. 


Clearly,  for  such  a  broad  scope,  the  aim  was  not  to  consider  each  aspect  in 


detail  but  rather  to  endeavor  to  concentrate  on  fundamental  problems  effecting 


the  feasibility  of  th«-  proposed  test;  also,  to  identify  areas  where  further  work 


either  experimental,  analytical  or  in  systems  design  was  necessary  or  would 


contribute'  to  the  further  development  of  the  project. 


The  BSD  is  a  new  type  of  in- situ  testing  apparatus.  Little  was 


Known  about  the  detailed  behaviour  of  the  device  until  the  analytical  aspects  of 


the  work  had  been  under  way  for  some  time.  Consequently  there  have  been  some 


fundamental  discoveries  that  have  influenced  significantly  both  the  configuration 


of  the  originally  proposed  apparatus  and  the  testing  proc.dures. 


A  brief  history  of  the  developments  leading  to  the  current  concept  of 


the  test  is  given  in  the  next  section.  The  following  sections  discuss  various 


aspects  of  these  developments  in  detail.  Also  discussed  are  the  problems  involved 


in  incorporating  these  developments  into  the  design  of  a  practical  field  measure¬ 


ment  system. 


f 


5.2  EVOLUTION  OF  THE  BSD  TEST 

The  principal  factors  affecting  the  configuration  of  the  testing 
apparatus  and  the  fundamental  nature  of  the  test  itself  have  emerged  largely 
from  the  analytical  aspects  of  the  work.  They  can  be  summarized  as  follows: 

(a)  It  was  found  that,  for  most  realistic  geometries  of  the  originally 

proposed  apparatus,  the  dominant  strain  induced  in  the  surrounding 
soil  by  the  coupling  cylinders  is  the  component  and  not,  as 

at  first  proposed,  the  component. 

(b)  Radiated  power  occurs  only  at  relatively  high  frequencies. 

For  a  significant  portion  of  the  power  applied  to  the  coupling 
cylinder  to  be  converted  to  radiated  power  (i.e.  to  generate  shear 
waves) ,  a  relatively  high  frequency  of  excitation  is  required. 

At  the  torsional  stiffnesses  calculated  for  the  in-situ  coupling 
cylinder,  an  extremely  high  mechanical  power  would  be  required 
by  the  system.  Much  of  this  power  would  be  effectively  reactive, 
that  is,  exchanged  within  the  system  and  not  contributing  to 
generating  seismic  waves. 

Regarding  the  strain  distribution,  it  was  found  that,  for  the  original 

configuration  of  the  apparatus  where  two  coupling  cylinder:  were  driven  by  a 

downhole  motor,  the  distance  between  the  cylinders  would  have  to  be  of  the 

order  of  half  the  radius  of  the  borehole  for  the  strains  to  be  comoarable 

9z 

with  the  component.  For  su<~h  a  geometry,  a  number  of  fundamental 
problems  arise.  Firstly  the  measurement  region  is  small  and  concentrated 
in  the  disturbed  regi<— ■;  also,  the  strain  field  becomes  very  complex 
and  the  edge  effects  introduced  by  local  yielding  of  the  soil  at  the  ends 
of  the  cylinders  influence  the  measurements  very  strongly.  Furthermore, 
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there  is  the  practical  problem  of  generating  and  coupling  the  drive  torque 
required  in  such  a  confined  space. 

Hence  the  Cg^  strain  concept  was  abandoned  and  the  analytical  and 
design  efforts  were  concentrated  around  a  device  which  was  to  be  optimized  for 
irducing  the  e  component  of  strain.  The  following  were  the  major  repercussions  / 
on  the  systems  design  work  of  this  change  in  the  orientation  of  the  project: 

(a)  A  two  cylinder  configuration  was  no  longer  essential.  For 

the  e  strains,  the  zone  of  high  strain  surrounds  the 
r9 

coupling  cylinder  and  extends  radially  outward.  Hence  a 
single  coupling  cylinder  driven  either  inertially  or  from  a 
remote  reaction  system  would  induce  the  required  strain 
equally  well.  As  a  consequence  of  this,  the  systems  design 
study  was  expanded  to  include  the  practical  considerations 
involved  in  developing  a  single  coupling  cylinder  borehole 
shear  device. 

(b)  Control  over  the  e  „  strain  distribution  is  not  achieved  bv 

rO 

varying  the  geometry  of  the  in-situ  apparatus,  at  least  not  to 
the  extent  of  the  original  component.  The  region  of 
high  strain  intimately  surrounds  the  coupling  cylinder  and 
thus  the  characteristics  of  the  coupling  and  the  disturbance 
induced  by  placement  are  potentially  limiting  factors  to  the 
integrity  of  the  test. 

(o)  At  the  commencement  of  the  project,  the  choice  between 
up/down  1  ^ie,  crosshole  and  poin"  measurement  or  any 
combination  therefore  was  a  deliberately  open  issue. 


Following  the  change  to  the  e^g  strains,  however,  the  up- 


hole  technique  was  immediately  eliminated  and  within  a  short 

period  of  time,  the  crosshole  method  was  also  ruled  out. 

Obviously  an  uphole  detector  would  have  to  rely  on  the 

unwanted  e.  component  and  for  the  crosshole  method,  the  analytical 
02 

solutions  showed  that  for  the  detector  to  be  measuring  in  the 
zone  of  high  strain,  it  would  have  to  be  impractically  close 
to  the  coupling  cylinder. 

The  effect  of  power  limitations  on  the  frequency  of  operation  of 
the  test  meant  that  it  was  not  possible  to  modify  the  strain  distribution 
by  operating  at  high  levels  of  radiated  power.  Here  again  a  blow  was  dealt 
to  the  option  of  the  crosshole  method.  Had  it  been  possible  to  operate  at 
frequencies  where  the  crosshole  distance  was  less  than  a  wave-length,  the 
method  would  have  been  worth  pursuing. 

Thus  the  BSD  test  evolved  into  a  semi-static  point  measurement. 

Material  properties  are  derived  from  the  torque  and  angular  displacement 
functions  measured  directly  at  the  coupling  cylinder.  A  sinusoidal,  or 
more  correctly,  a  periodic  excitation  is  applied  to  the  coupling  cylinder/s 
in  the  form  of  a  feedback  controlled  torque.  The  input  to  the  system  is 
monotonically  increased  to  progressively  load  the  soil  throughout  the 
strain  range  of  interest.  Both  torque  and  the  resultant  angular  displacement 
functions  are  recorded  throughout  the  loading  history  and  used  subsequently 
to  derive  the  u-t  relationship  for  the  soil  under  test.  The  original 
concept  of  two  coupling  cylinders  driven  by  a  downhole  actuator  is  still 
discussed  in  this  rer  _c,  but  it  is  dealt  with  as  a  technique  which 
produces  two  simultaneous  point  measurements. 
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in  the  low  level,'  near-elastic  range.  In  the  high  range,  the  strains 
indicated  by  the  rotation  signal  using  an  elasvic  relation  would  always 
be  smaller  than  the  highest  strain  in  the  surrounding  soil.  Hence  the 
danger  would  exist  that  the  apparatus  could  advance  the  soil  too  rapidly 
through  the  loading  history.  Under  these  conditions,  the  ri *  '  of  increase 
of  strain  pfer  cycle  of  load  could  be  large  enough  to  limit  seriously  the 
amount  of  useful  data  obtainable  from  the  measurement. 

If,  on  the  other  hand,  torque  is  used  as  the  control  variable,  it 
appears  more  likely  that  the  appropriate  level  of  loading  can  be  estimated 
better  and  that  it  can  be  more  closely  controlled  throughout  the  loading 
history.  This  is  felt  to  be  so  because  there  will  always  be  a  direct  relation 
between  the  mean  shear  stress  at  the  coupling  mechanism/soil  interface  and  the 
applied  torque.  The  relation  of  this  shear  stress  to  shear  strain  will  depend 
heavily  on  the  soil  properties,  but  a  reasonable  upper  bound  can  be  easily 
estimated  based  on  competent  geotechnical  data  for  the  site.  This  upper 
bound  is  determined  from  the  shear  strength  of  the  soil,  since  the 
apparatus  cannot  couple  torques,  with  their  corresponding  shear  stresses, 
which  exceed  the  shear  strength  of  the  soil.  Thus,  in  order  to  ensure 
loading  of  the  soil  at  high  strains,  it  is  reasonable  to  assume  that  the 
loadings  should  extend  to  yield  of  soil.  As  mentioned  previously,  estimation 
of  this  bound  is  relatively  easy  given  good  site  data;  not  only  is  the 
procedure  routine,  but  much  practical  experience  coupled  with  field  test 
results  is  available  to  assist  in  the  calculation. 

Given  an  upper  bound  or  maximum  value  for  the  applied  torque,  the 
operation  procedure  would  be  to  start  the  loading  history  at  zero  and,  under 
controlled  torque,  increase  the  load  at  an  appropriate  rate  to  the  calculated 


yield  value.  The  zero  starting  value  would  of  course  in  practice  correspond 
to  the  resolution  of  the  system.  This  may  be  determined  by  backlash,  striction 
or  hysteresis  in  mechanical  components;  or,  in  a  higher  performance  system, 
by  random  noise  introduced  by  one  or  more  of  the  system  components. 

It  may  be  appropriate  to  modify  the  feedback  signal  in  the 
control  loop  according  to  the  rate  of  increase  of  strain  but,  fundamentally, 
the  system  would  be  described  as  a  single-input,  single-output  torque 
control  system. 


5.4  RECOVERY  OF  SOIL  PROPERTIES  FROM  FIELD  DATA 


The  basic  purpose  of  the  development  and  operation  of  the  BSD 
is  the  generation  of  a  curve  relating  the  secant  shear  modulus  to  the 
strain  level.  The  present  section  is  hence  of  the  outmost  importance 
since  it  is  the  one  which  provides  the  link  between  the  in-situ  operation 
and  measurements,  and  the  construction  of  the  curve  sought. 

It  has  already  been  shown  that,  at  the  design  frequencies  of 
the  apparatus,  the  associated  wave  lengths  are  long  and  the  real  part 
of  the  dynamic  strain  distribution  approaches  the  static  one.  These 
facts  permit  proposing  a  relatively  simple  method  of  deriving  the  p-e 
curve. 


-  4 .  j  Proposed  Procedure 

From  Table  4  the  torque  required  to  induce  a  rotation  0  on 
a  long  cylinder  in  elastic  soil  is; 


T  -  4nr  ji  h© 
EL  o  o 


(1) 


It  was  also  shown  (Section  3.3)  that,  if  the  cylinder  can  no 
longer  be  considered  long,  such  torque  becomes: 


T 


EF 


2 

2iir  p  h©F 
o  o  E 


(2) 


where  F  is  an  influence  function  given  in  Figures  5  to  8. 
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This  means  that,  in  elastic  soil,  the  solution  for  a  finite  cylinder 
can  be  obtained  from  that  for  an  infintely  long  cylinder  as: 


T  =  T  F 
EF  EL  E 


For  the  cases  of  interest,  it  was  justified  in  Section  4.3.3  the 

fact  that  F£  can  be  considered  to  depend  on  —  alone.  For  the  present  design, 

o 

=  5,  for  which  value  F  =  2.31.  In  the  case  of  a  long  cylinder  in 

£ 

o 

soil  with  a  hyperbolic  stress-strain  law,  the  torque  at  small  dimensionless 
frequencies  (Section  3.1  and  4.3.1)  can  ne  expressed  as: 


Tm  =  2*r  ht  1  -  e 
HL  o  ml 


It  was  also  shown  in  Section  4.3.3  that  the  equivalent  of  eq.  (4) 
for  a  short  cylinder  can  be  written: 


2  /  .3sL\ 

r  nT  (  1  -  e  t  I 
o  m  \  m  / 


that  is: 


T  =  T  .F  (< 

HF  HL  H 

where  F  has  been  deduced  numerically  (Table  6)  and  appears  to  remain 
H 

practically  constanty  over  the  range  of  strains  of  interest. 


Since  F  appears  to  be  approximately  the  same  for  elastic  and  hyperbolic 
soil  models,  even  at  high  strains,  let  us  assume  for  the  moment  that  it  will 
also  be  the  same  for  other  non-linear  soil  models.  This  assumption  will  be 
verified  in  the  next  subsection. 


■  %  mm  ‘SKissg 


As  mentioned  in  the  previous  section,  the  field  procedure  recommended 
consists  of  cycling  the  torque  at  gradually  increasing  levels  starting 
from  the  lowest  values  allowed  by  the  limit  of  resolution  of  the  apparatus. 


When  the  maximum  strains  are  small,  the  initial  modulus  p  is 
easily  determined  from  eq.  5: 


THF  =  Wr, 


fx  -  (i  -  % 

o  m[  \  Tir 


=  2irr  hp  F 
o  o  e 


as  could  be  expected  from  eq.  (1).  Hence: 


2irr  hF  0 
o  E 


For  larger  strains,  it  must  be  noticed  that,  in  order  to  obtain  the 
relationship  between  T  and  p  at  a  given  strain  level,  the  form  of  the  cons¬ 
titutive  law  must  be  known.  The  problem  is  hence  implicit.  The  following 
procedure  is  proposed. 


Suppose  the  stress-strain  law  is  truly  hyperbolic.  Then,  once 

p  is  known,  a  second  test  at  higher  0  is  sufficient  to  determine  x  from 
o  in 

eq.  (5)  and  hence  to  provide  p  vs  e  at  all  values  of  e.  Indeed: 


h 

1+  —  e 

T 

m 


where  p  is  already  known,  as  ror  t  't  can  be  found  solving  by  iteration 
o  m 

the  implicit  equation  (5)  ,  which  can  be  v/ritten  in  the  form: 
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T 

m 


_ T _ 

2^0 

irx2h  F  fl-e  T  1 

o  H  ]_  raj 


(9) 


A  first  approximation  for  is  easily  obtained  by  adding  one 
term  to  the  series  expansion  previously  used: 


which  yields: 


T  = 
m 


2  ii  0 

o 


1  + 


HF 


27rr  hu  F  0 
o  o  H 


(lO) 


Such  an  approximation  is  always  sufficient  for  the  first  iteration 
due  to  the  good  convergence  characteristics  of  eq.  (9) . 


Even  if  the  stress-strain  law  were  not  truly  hyperbolic,  the 
procedure  for  obtaining  the  ii-e  curve,  as  described  above,  would  probably 
sufficient,  particularly  for  the  lower  strains  (probably  up  to  about  10  . 


What  such  procedure  would  obtain  is  essentially  a  best  fit  of  the  real 
stress-strain  law  by  means  of  a  hyperbola.  In  order  to  avoid  imposing  a 


hyperbolic  stress-strain  law  onto  the  soil,  the  following  improvement  to  the 
procedure  described  is  proposed. 


Suppose  a  test  at  low  strains  has  already  been  conducted  and 

M  is  known.  The  next  test  is  carried  out  at  slightly  higher  strains  and 
o 

allows  the  computation  of  a  value  by  means  of  eq.  (9)  and  using  eq.  (10) 

as  a  first  approximation.  The  rotation  is  again  increased;  the  former  value 

of  t  is  now  used  as  a  first  approximation  to  solve  eq.  (9) .  This  is  continued 
ra 

until  the  soil  fails. 

Clearly,  at  each  rotation  level,  the  stress-strain  law  of  the  soil 

has  been  fitted  with  a  hyperbolic  law  over  the  range  of  strains  associated  with 

the  rotation  imposed.  Such  hyperbolic  law  allows  the  derivation  of  the  p-e 

relationship  for  all  values  of  e  smaller  than  the  maximum  strain  imposed. 

It  is  clear  that  when  very  high  strains  are  involved,  the  approximation  will 

be  less  satisfactory  than  when  all  strains  are  low;  this  is  so  because  in  the 

former  case  a  larger  portion  of  the  stress-strain  law  is  fitted  with  the 

hyperbola.  This  means  that  the  value  of  p  corresponding  to  a  strain  z  will 

be  optimum  if  derived  from  the  value  of  i  obtained  during  the  test  in  which 

m 

z  was  the  maximum  strain.  Obviously,  if  the  material  is  hyperbolic,  all 

the  t  values  will  be  identical, 
m 


In  summary  the  secant  shear  modulus  u  at  each  strain  level  z, 
based  on  the  torque  T  required  to  induce  a  i-otation  G,  is  taken  as: 


p 


JL  + 


Poe<e) 
r  (S) 


(11) 
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where  e  and  t^  are  obtained  from  6  via  the  hyperbolic  relationships. 


5.4.2  Verification  of  the  Proposed  Procedure 

In  order  to  verify  the  adequacy  of  the  procedure  described,  some 
"numerical  field  tests"  have  been  performed.  The  purpose  of  these  tests  is  to 
reproduce  numerically  the  field  operation.  The  soil  is  assumed  to  have  a  certain 
stress-strain  law;  in  the  soil,  the  BSD  is  operated  numerically,  thus  producing 
torque  and  rotation  values.  The  procedure  described  in  the  latter  section 
is  then  used  to  generate  a  y-t  curve  which  is  compared  to  the  original  one 
in  order  to  assess  the  errors  made  in  the  process.  Tests  with  hyperbolic 
and  non-hyperbolic  soils  have  been  performed  and  are  described  in  this  section. 

It  should  be  reminded  however  that  the  examples  included  in  the  section 
provide  an  upper  bound  rather  them  a  precise  estimate  of  the  errors  involved. 

This  is  so  because  the  soil  in  the  field  will  respond  exactly  according  to 
its  true  constitutive  law,  while  the  numerical  soil  includes  errors  in  the 
response  generated  by  the  finite-element  procedure  itself.  In  a  way,  the 
present  validation  compounds  the  errors  of  the  finite-element  method  with  those 
of  the  procedure  for  generation  of  the  y-e  curve  from  T-0  data. 

The  first  series  of  "numerical"  field  tests  were  run  with  a  hyperbolic 
soil.  Input  data  were  as  follows: 


h 

W 


m 


0.1  m 
0.5  m 

1.8x10^  N/m^ 

5x]^5  M/m2 


(o  =  100  rad/sec 
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From  the  elastic  run  at  strains  smaller  than  lO  ,  p  was  exactlv 

o 

recovered,  as  could  be  expected  from  the  results  shown  in  Section  4.2. 

Then  a  second  test  was  run  with  identical  input  data  except  that  the 
displacement  at  the  BSD-soil  interface  was  now: 

-4 

u  =  10  m 
o 

corresponding  to  a  rotation: 

0  =  10  ^  rad 

The  torque  "measured"  by  the  computer  program  was: 

3 

T  =  9.2x10  Nm 

Solving  then  eq.  (9)  by  iteration: 

T  =  4. 83x10 5  N/m^ 
m 

In  Figure  33  the  input  curve  is  shown  together  with  the  one  derived 
from  the  test.  As  could  be  expected,  the  agreement  is  very  good  for  a 
hyperbolic  stress-strain  law  since  the  approximation  is  also  hyperbolic. 

The  true  test  of  the  procedure  is  to  verify  that  it  recovers  a 
p-e  curve  even  if  the  latter  does  not  correspond  to  a  hyperbolic  stress- 
strain  law.  The  same  configuration  of  the  above  example  was  used  again; 
in  this  test  however,  the  stress-strain  law  was  now  defined  by  the  thick 
line  in  Figure  34.  It  can  be  noticed  that  such  curve  intersects  the 
hyperbolic  ones  and  includes  a  horizontal  segment  at  the  end.  Besides  the 
elastic  run,  tests  were  performed  at  the  following  displacements  and 


rotations: 


-  INPUT  CURVE  (HYPERBOLIC,  rm*8xlO®  N/m2) 

- F.E.  SOLUTION  (HYPERBOLIC,  rm  «  4.83x10®  N/m2) 


U  =  3.5x10  n 
o 


It  is  not  difficult  to  reproduce  the  noculus  values  at  strain 
levels  at  which  the  initial  value  has  not  deteriorated  excessively.  The  cosi 
severe  tests  axe  the  cases  shown  here  where  the  secant  modulus  is  50% 
or  less  of  the  initial  value .  Table  7  presents  a  sussary  of  the  results 
obtained  by  the  program  and  the  values  calculated  by  the  procedure 
explained  in  Section  5.3.1. 


Table  7.  Results  of  the  non-hvoerbolic  tests 
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The  u-c  curve  given  as  input,  the  four  hyperbolic  laws  used  in  the 
successive  approximations  and  the  four  points  in  which  the  hyperbolic  laws 
are  used  to  obtain  a  p-e  point  are  all  plotted  together  in  Figure  34. 


As  can  be  seen,  the  agreement  is  better  than  could  be  expected  from  an 
approximate  method,  particularly  since  the  errors  compound  those  in 
the  F.E.  secant  linerarisation  with  those  in  the  derivation  of  p-e 
values  via  a  hyperbolic  approximation.  As  limitation,  it  should  be 
pointed  out  that  it  is  probably  dangerous  to  use  this  scheme  when  the 
values  of  the  secant  modulus  are  below  10%  of  its  initial  value. 


5.5  BSD  TESTING  HARDWARE 


5.5.1  General 

Figure  35  is  a  schematic  diagram  comparing  the  in-situ  hardware 
for  the  two-cylinder  and  single  cylinder  BSD's.  Both  devices  employ 
exactly  similar  coupling  cylinders  since,  from  the  design  or  systems 
specifications  point  of  view,  it  is  immaterial  whether  they  are  driven  from 
the  surface  or  from  a  downhole  actuator.  For  all  realistic  geometries  of 
the  apparatus,  the  torque-displacement  relations  derived  from  each  coupling 
cylinder  can  be  considered  to  be  independent  measurements.  There  are  no 
special  considerations  such  as  modified  stress  distributions  or  design 
torques  that  influence  the  design  of  the  coupling  cylinder  for  either  the 
single-/  or  two-cylinder  measurement  systems. 

5.5.2  The  Self  Boring  Attachment 

The  self  boring  attachment  is  shown  fitted  to  both  systems.  The 
purpose  of  this  equipment  is  to  remove  and  convey  through  the  apparatus 
sufficient  material  to  enable  it  to  enter  the  virgin  soil  below  the 
borehole  casing  with  an  absolute  minimum  of  disturbance.  The  helical  auger 
which  is  used  to  transport  the  soil  cuttings  through  the  apparatus  is 
concentrically  located,  passing  through  the  coupling  mechanism  and,  in  the 
case  of  the  two-cylinder  arrangement,  passing  also  through  the  actuator. 
Provision  of  this  free  and  uninterrupted  central  bore  through  the  apparatus 
has  been  a  major  design  constraint,  particularly  for  the  two-cylinder 
arrangement.  The  access  bore  for  the  auger  must  be  of  sufficient  diameter 
to  allow  a  reasonable  s5  of  auger  and  yet  not  to  limit  too  severely 
the  radial  space  available  for  the  torque  motor.  The  less  radial 
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space  available  for  the  actuator,  the  longer  it  becomes  in  order  to  develop 


a  given  torque.  The  auger,  on  the  other  hand,  must  transport  a  volume  of 
soil  cuttings  of  the  order  of  that  of  the  in-situ  system.  For  the  downhole 
acutator  system,  the  problem  is  accentuated  since  the  apparatus  is 
considerably  greater  in  volume  than  the.  single-cylinder  system,  with 
correspondingly  greater  distances  over  which  the  auger  must  convey  the  soil 
cuttings.  Clearly,  in  the  design  of  a  two-cylinder  BSD,  a  careful  trade¬ 
off  between  the  conflicting  requirements  of  the  components  of  these  two  system 
would  be  required.  Both  systems  are  shown  equipped  with  a  catch  tank  to 
store  the  soil  cuttings.  In  the  case  of  the  two-cylinder  system,  it  may  be 
more  convenient  to  use  a  bailer  tank  system  where  the  tank  is  periodically 
raised  up  the  hole  in  order  to  dispose  of  the  cuttings. 

5.5.3  Routing  of  Power,  Instrumentation  and  Other  Services 

A  major  consideration  in  the  design  of  the  field  equipment  will 
be  the  routing  of  the  various  cables  and  hydraulic  hoses  which  convey 
services  to  the  downhole  components.  A  convenient  means  of  routing  the 
instrumentation  cables  and  a  limited  amount  ot  hydraulic  power  is  afforded 
in  the  coupling  mechanism  by  the  space  at  the  bottom  of  the  grooves  in 
which  the  shear  coupling  keys  fit.  These  services  are  conveyed  by  means  of 
minature  cables  and  small  bore  hydraulic  tubes  which  easily  fit  into  the 
clearance  space  at  the  bottom  of  the  grooves  (see  Fig  36) .  In  the  case  of 
the  two  cylinder  mechanism,  the  actuator  power  hoses  will  require  considerably 
more  space.  The  only  solution  so  far  envisaged  here  is  that  they  encroach 
on  the  central  space  provided  for  the  helical  auger.  Tables  8,  9  and  10  are 
lists  of  the  services  coupled  between  the  uphole  and  downhole  systems.  The  list 
is  divided  between  the  two  alternative  systems,  i.e.:  the  double  and  single 
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cylinder  USDs  and  then  subdivided  into  individual  system  components. 
Certain  features  such  as  the  instrumentation  for  the  self-boring  equipment 
may  not  ultimately  be  deemed  necessary;  they  are  however,  included  in  this 
list  for  completeness. 

Table  8. _ Power  requirements  of  the  BSD  apparatus 


Torsional  Drive  Power 


Single  Cylinder  BSD  Two  Cylinder  BSD 

Rigid  torsion  tube  sectioned  High  pressure  hydraulic 

into  2  metre  lengths.  The  tube  supply  coupled  via  two 
is  thick  or  thin-walled  depending  30  mm  nominal  diameter 


upon  whether  oscillatory  or  direct  hoses, 
coupling  is  employed.  Outside  rating: 


Typical  supply 
2 . 5xl07  Pa, 


diameter  of  tubing:  100  to  15C  mm,  (250  bars)  at  0.5  litres/ 
max  torque  transmitted:  104  Nm.  second. 


Lateral  Jacking  and  Two  micro-bore  (3  diameter  maximum)  high  pressure  hydraulic 

and  Retracking  Power  lines  carrying  bi-directional  flow.  Maximum  jacking 

pressure:  150  bars,  maximum  compressible  displaced  volume 
of  fluid:  0.035  litres 


Hydraulic  fluid  conveyed  by  Hydraulic  fluid  conveyed 

tubing  attached  to  the  wall  of  by  flexible  high  pressure 
the  torsion  tube.  Quick  hose  fused  at  intervals 

release,  self  sealing  couplings  along  main  umbilical  line, 
installed  at  torsion  tube  joints. 


Auger  Drive  Power 


Sectioned  torsion  rod  fitted 
with  flush  screw  joints.  Web- 
type  centering  brushes  are  fitted 
at  torsion  tube  joints  to  main¬ 
tain  concentricity  with  torsion 
tube. 


Sectioned  torsion  rod 
running  in  bushed  radial- 
web  type  centering  fixtures 
to  maintain  concentricity 
with  borehole  casing. 


Cutting  Tip  Thrust  Applied  from  surface,  coupled  Applied  either  by  temporarily 

directly  through  torsion  tube.  installed  jacking  unit  acting 

against  borehole  casing  or 
coupled  from  surface  via  auger 
drive  rod. 

Table  P  Instrumentation  of  the  BSD  appa; ~atus 

Single  Cylinder  BSD  Two  cylinder  BSD 

Coupling  Cylinder  Drive  Multi -element  strain-gauge  torque  Multi-element  strain-gauge 

transducer  configured  on  torsion  torque  transducer  configured 
tube  close  to  coupling  auger.  on  the  actuator  shaft. 

Multi-element  strain  gauge  flexural  web  arrangement  coupled 
via  rigid  rod  to  surface  reference  frame.  Web  is  coupled 
to  inside  bore  of  coupling  cylinders  by  means  of  radially 
expending  grips. 


Coupling  Oy) inder 
Rotation 


Coupling  Cylinder  *• 


Alternatively  high  resolution  RVDT  or  LVDT  (Rotary  or  Linear 
Variable  Differential  Transformer)  coupled  as  above  to 
surface  reference  frame. 

(radius  arm  mechanism  would  be  required  with  the  linear 
device) . 


Shear  Coupling  Key 
Displacement. 


LVDT  fitted  inside  fixec  keys  at  both  ends  of  the  coupling 
cylinder/s.  24  LVDT's  per  coupling  cylinder. 


Lateral  Jacking 
Pressure 


Miniature  strain  gauge  diaphrara  pressure  tranducer  plumbed 
into  the  hydraulic  system  as  close  as  possible  to  the 
coupling  cylinder/s. 


Cutting  Tip  Thrust 


"Doughnut"  type  low  profile  multi-element  strain  gauge  force 
transducer  me  unted  behind  the  rim  blade  of  the  self-boring 
head. 


Auger  and  Breaker  Tip 
Drive  Torque 


Probably  most  conveniently  instrumented  at  the  surface 
(at  least  for  preliminary  trials) .  Conventional  slip-ring 
coupled  strain  gauge  torque  transducer. 


Axial  (Placement) 
Displacement 


Again,  probably  most  conveniently  instrumented  at  the  surfaces. 
Conventional  long-stroke  linear  potentiometer  displacement 
transducer . 


Pore  Water 
Pressure  and  Total 
Stress  Gauge 


Flush  contoured,  strain-gauge  diaphram  type  pressure 
transducers,  inset  int.  thick  conic  section  of  rim  blade 
due  to  axial  thrust. 


Plane  Strain 
Penetrometer.  Force/ 
displacement  measure¬ 
ments.  (see  Section 
5.6.2) 


One  LVDT  and  one  strain  gauge  force  transducer  channel  per 
plane  strain  penetrometer  probe. 


Table  10.  Miscellaneous  services  for  the  BSD  apparatus 


Single  Cylinder  BSD 


Two  Cylinder  BSD 


Bailer  Tank 
Lifting  Gear 


Bailer  tank  not  necessary 
for  Single  Cylinder  BSD 


Wire  line  used  to  hoist 
bailer  tank  uo  hole. 


Deadload  Relief 


Deadload  of  torsic- 
tube  can  be  relieved  by  pull 
applied  directly  to  tubing  at 
surface. 


Wire  line 


m 


Flushing/Lubricating 

Fluid 


If  required  temporarily  or  continuously  during  self 
boring,  flushing/lubricating  fluid  can  be  pumped  down 
the  centre  of  the  auger  drive  shaft. 


5.5.4  Practical  Aspects  Of  The  Single  And  Two  Cylinder  BSPs 

The  most  significant  difference  between  the  single  and  two 
cylinder  BSDs  in  terms  of  handling  and  placement  procedures,  is  the  means 
by  which  power  for  the  torsional  drive  is  transmitted  down  hole.  The 
single  cylinder  system  employs  a  torsion  tube  which  is  driven  from  a 
surface  mounted  actuator.  Reaction  torque  is  provided  by  a  system  of  ground 
anchors  and  the  torsion  tube  can  be  used  either  as  a  direct  coupling  or  as 
part  of  a  resonant  coupling  which  would  be  able  to  magnify  the  torque 
developed.  Drive  torque  for  the  coupling  mechanism  is  derived  from  a  motor 
or  set  of  motors  mounted  on  the  reaction  frame  and  coupled  directly  to  the 
torsion  tube.  The  system  is  shown  in  Fig  32  as  an  electrohydraulic  control 
system  which  would  typically  be  comprised  of  a  two  stage  flapper/spool 
servo  valve  driving  a  vane  type  torque  motor.  Hydraulic  power  would  be 
developed  at  the  system  pressure  by  a  screw  or  piston  type  pump.  Since  the 
dynamic  load  history  to  be  applied  to  the  soil  in  any  one  test  is  relatively 
short,  the  possibility  exists  for  accumulator  storage  of  the  total  energy 
for  each  test.  In  this  case,  hydraulic  power  could  be  developed  slowly 
from  a  high  pressure  pump  which  might  ultimately  be  powered  by,  say,  a 
vehicle  battery.  This  latter  feature  obviously  applies  to  both  the  single 
and  two  cylinder  BSD's. 


The  single  cylx..Jer  BSD  system  however  is  not  limited  to  the  use 
of  hydraulic  actuators  as  is  the  two  cylinder  device.  The  actuator  is 
surface  mounted  and  hence  is  not  constrained  by  problems  of  limited  space. 
It  may  be  more  appropriate  to  use  electric  torque  motors  with  switching 


drive  for  the  high  torque  *nd  linear  drive  for  smooth  control  at  low 


torque.  A  decision  on  the  type  of  drive  most  suited  to  this  system  is 


however,  not  critical  to  this  design  study  since  it  does  not  influence 
significantly  the  concept  of  the  single  cylinder  BSD.  It  is  better  relegated 
to  a  subsequent  phase  of  the  BSD  development  following  a  detailed  study  of 
these  devices. 

The  torsion  tube  in  the  single  cylinder  BSD,  fulfills  a  number  of 
roles.  Firstly,  and  most  obviously,  it  transmits  the  main-  drive  torque  for 
the  coupling  mechanism.  Secondly,  it  forms  a  catch  tank  to  retain  soil 
cuttings.  Finally  it  allows  direct  transmission  of  axial  thrust  to  the 
cutting  tip  during  the  self  boring  operation.  The  chief  penalty  inherent 
in  this  simple  configuration,  however,  is  the  necessity  for  couplings  in 
the  system.  The  torsion  tube  must  obviously  be  sectioned,  the  most 
appropriate  length  being  determined  by  factors  such  as  convenience  of 
handling  in  the  field  and  by  the  resolution  or  minimum  spacing  of  the 
measurement  locations.  These  joints  in  the  torsion  tube  must  be  able  to 
transmit  the  design  drive  torque  without  introducing  significant  backlash 
compared  with  the  rotation  of  the  coupling  mechanism.  This  may  not  seem  at 
first  sight  to  be  a  reasonable  requirement  since  the  drive  torque  and 
angular  displacement  of  the  coupling  cylinder  are  instrumented  directly  at 
the  measurment  location.  The  problem  arises,  however,  in  the  control 
system  which  under  torque  control,  experiences  a  loss  or  deterioration  in 
the  linearity  of  the  transfer  characteristics  when  the  torque  is  small  or 
at  the  point  of  reversing.  This  region  is  precisely  the  point  where  a 
well-behaved  system  is  required,  at  least  for  the  lower  strain  stage  of  the 
test.  The  high  loop  gains  required  to  compensate  for  such  non-linearities 
would  certainly  add  to  any  problems  with  stability  of  the  control  system. 


The  two  cylinder  system,  on  the  other  hand,  is  not  rigidly 
coi pled  to  the  surface.  It  is  suspended  on  an  umbilical  line  which  is 
comprised  of  hydraulic  power  and  instrumentation  cables  linked  to  a  strain 
relief  wire. 

Drive  torque  for  the  coupling  cylinders  is  derived  from  a 
centrally  located  torsional  actuator  coupled  directly  to  the  two  cylinders. 

Some  of  the  practical  difficulties  inherent  are  two  cylinder  BSD  configurations 
have  already  been  discussed,  notably,  those  associated  with  the  incorporation 
of  the  self-boring  equipment.  The  actuator  in  the  two  cylinder  BSD  would 
almost  certainly  have  to  be  a  hydraulic  motor.  The  extremely  limited 
radial  space  and  the  high  drive  torque  required  of  the  actuator  are 
constraints  that  only  a  high  pressure  hydraulic  system  could  meet  unless 
some  system  of  gearing  is  also  incorporated.  In  the  absence  of  radial 
space  the  length  of  the  motor  increases,  leading  to  further  drawbacks  with 
the  two  cylinder  BSD.  It  is  only  reasonable  to  expect  that  the  coupling 
mechanisms,  upon  expanding  in  the  locking-on  operation,  will  become 
slightly  misaligned.  This  will  result  in  bending  moments  applied  to  the 
motor  shaft  which,  of  course,  increase  with  the  length  of  the  motor.  Thus 
it  is  likely  that  a  system  of  coupling  would  be  required  between  the  motor 
and  the  coupling  cylinders  to  isolate  or  reduce  this  bending.  Those 
couplings,  however,  should  not  significantly  degrade  the  torsional 
stiffness  of  the  driving  system,  otherwise  reactive  power  consumption  of  the 
system  would  be  increased.  Reactive  power  is  here  defined  as  the  power 
observed  into  strain  energy  in  the  BSD  system.  It  is  wasted  power  since  the 
electrohydraulic  systems  Ao  not  recover  power  returned  from  the  load,  and 
the  low  stiffnesses  which  are  characteristic  of  their  systems  limit  the 
performance  obtainable  from  the  control  system. 


A  fundamental  problem  of  the  two  cylinder  BSD  concept  is  that  of 
limiting  the  dynamic  range  of  operation.  As  discussed  above,  there  is 
little  choice  in  the  type  of  actuator  that  can  be  used  and  little  space 
available  in  which  to  tandem  devices,  or  to  incorporate  alternative  drive 
systems  to  extend  the  dynamic  range  of  operation.  Furthermore,  the  actuator  is 
operating  under  relatively  adverse  conditions  which  may  include  bending  and 
lateral  shaft  loading.  The  dynamic  range  of  controlled  torque  available 
from  the  device  corresponds  directly  to  the  range  of  shear  stress  coupled 
to  the  soil  and  hence  to  the  range  of  soil  properties  that  the  BSD  can  cope 
with. 

Common  to  both  systems  is  the  requirement  to  relieve  dead  weight 
loading  in  order  to  avoid  inducing  unwanted  shear  stresses  in  the  soil. 

This  dead  weight  is  partly  contributed  uy  the  in-situ  apparatus,  and,  in 
the  case  of  the  single  cylinder  system,  by  the  torsion  tube  string.  If 
instrumentation  is  incorporated  in  the  downhole  equipment  to  measure  the 
cutting  tip  thrust,  it  may  well  serve  also  as  a  monitor  for  the  deadload 
applied  to  the  coupling  mechanisra/s. 

The  technique  proposed  for  the  measurement  of  angular  displacement 
is  again  the  same  for  both  the  single  and  two  cylinder  BSDs  (See  section  5.5.4) 
The  displacement  reference  for  the  system  is  derived  from  a  rigid  rod 
lowered  through  the  apparatus  and  attached  to  a  reference  frame  at  the 
surface.  This  rod  does  not  couple  directly  to  the  apparatus,  it  is  linked 
via  the  angular  displacement  transducer  and  an  expansion  mechanism  to  the 
central  access  bore  in  the  coupling  mechanism/s. 


The  following  points  discuss  in  more  detail  the  relative  merits 
of  the  two  systems  (SC=single  cylinder,  DC=double  cylinder). 

(a)  PLACEMENT:  The  SC  requires  a  substantially  shorter  self-bored  length 
This  expedites  and  reduces  the  cost  of  the  field  procedures.  Also,  in 
the  DC  it  is  almost  certain  that  the  upper  cylinder  will  be  required 
to  contact  soil  somewhat  more  disturbed  than  the  lower  one. 

(b)  MEASUREMENT  RESOHJTION  AND  EFFECTS  OF  SOIL  PROPERTIES:  The  SC  provides 
a  single  test  at  a  time,  performed  at  one  depth  and  repeatable  at 
approximately  1  metre  intervals.  The  DC  could  provide  two  simultaneous 
tests  at  2  metres  from  each  other,  repeatable  every  3  or  more  metres. 
When  each  shoe  is  in  a  different  layer,  or  at  shallow  depths  in  soils 
with  depth  dependent  properties,  the  interaction  between  the  two  shoes 
may  be  difficult  to  analyze#  particularly  when  the  constitutive  non- 
linearities  induce  deviations  from  the  sinusoidal  shape  of  the  actuator 
output.  In  other  words,  the  input  to  one  shoe  is  affected  not  only 

by  the  actuator  and  the  soil  surrounding  a  particular  shoe,  but  also 
by  the  soil  around  the  other  one.  Furthermore,  if  the  material  at 
one  shoe  is  significantly  weaker  than  that  at  the  other  one,  large 
strains  will  not  be  obtainable  at  that  shoe.  Hence,  the  DC  performance 
grows  worse  as  gradients  of  soil  strengths  increase.  This  problem 
is  not  shared  by  the  SC  configuration. 

(c)  DYNAMIC  RANGE  OF  OPERATION:  T^e  large  strain  operations  require 
torques  which  generate  stresses  close  to  the  producing  soil  failure. 
Since  shear  strengths  are  likely  to  increase  more  or  less  linearly 
with  depth,  the  torque  required  will  have  to  change  over  a  range  at 
least  similar  to  that  of  depths  of  operation  {2  to  50  metres) . 
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Since  it  would  be  ideal  to  be  able  to  use  the  BSD  for  a  large  range 
of  soils  (from  shallow  and  soft  to  deep  and  stiff) ,  and  also  since 
a  tain i cum  range  of  strains  has  to  be  examined  at  each  depth.-  it 
is  apparent  that  a  wide  dynamic  range  of  controlled  output  j.s  an 
important  requirement  of  the  BSD.  The  advantages  of  the  SC  with 
respect  to  the  DC  are  obvious  in  this  respect.  The  hydraulic  no*_v_r 
in  the  DC  system  would  have  a  probable  dynamic  range  of  torque  on 
the  order  of  or  little  above  two  decades.  This  system  is  very 
difficult  to  improve  or  change  since  it  is  specially  developed 
and  located  within  the  confines  of  a  borehole.  The  SC  on  the  other 
hand  is  powered  by  a  motor  on  the  surface.  This  actuator  is  not 
confined  by  the  borehole  and  allows  the  freedom  to  choose 
between  a  number  of  systems  (electro-hydrualic  or  electromagnetic) , 
or  indeed  a  combination  of  systems.  By  this  means  the  dynamic  range 
can  be  extended  and  much  less  special  systems  development  is  involved. 

(d)  FREQUENCY  OF  OPERATION:  The  upper  frequency  attainable  with  the  SC 
system  approaches  50Hz,  the  only  theoretical  limitation  being  given 

by  the  resonant  frequency  of  the  coupling  between  the  resonant  mechanism 
and  the  BSD.  However,  in  the  case  of  the  DC  the  torque-frequency  curve 
is  fixed;  the  maximum  frequency  is  determined  by  the  peak  angular 
amplitude,  the  torque  and  the  available  power.  This  frequency  is  not 
expected  to  approach  10Hz  at  the  maximum  design  torque. 

(e)  HANDLING  Ht  THE  FIELD:  The  SC  requires  a  special  torsional  tube 
coupling  with  backlash-free  connections  in  order  to  transmit  the  torque. 
This  tubing  has  to  be  transported  and  handled  carefully.  The  DC  on 

the  other  hand  is  -  -ore  compact  system  which  is  suspended  and  powered 
by  means  of  a  flexible  line.  However,  the  DC  is  a  far  heavier  single 
unit  tham  the  SC  increasing  problems  with  handling  and  transportation. 
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(f)  HARDWARE  DEVELOPMENT:  The  SC  can  be  more  readily  produced  and 
involves  less  special  component  development  that  the  DC.  In 
particular,  the  TC  requires  a  specially  developed  hydraulic  motor 
which  is  both  costly  and  complex  in  design. 
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5.6  PLACEMENT 

Following  discovery  of  the  fact  that  the  allowable  dimensionless 
frequencies  were  small  and  the  BSD  could  not  properly  be  considered  as  a 
wave  generating  device,  it  became  evident  that  placement  would  be  a  critical 
factor  affecting  the  test.  This  is  because  the  zone  of  high  strain  is  highly 
localized  at  the  coupling  cylinder.  Hence,  it  becomes  important  to  minimize 
the  disturbance  to  the  in-situ  stress  state  and  fabric  of  the  soil 
immediately  surrounding  the  apparatus  is  required. 

To  assist  with  these  problems,  the  following  features  were  incorporated 
into  the  design  of  the  BSD. : 

(a)  A  self-boring  mechanism  which,  as  an  integral  part  of  the 

downhole  equipment,  is  arranged  to  advance  into  the  soil  ahead  of  the 
coupling  cylinder. 

(b)  A  precision  engineered,  multi-element  concept  has  been  used  for 
the  design  of  the  coupling  cylinder,  in  which  the  device  is  capable 

of  making  some  sensitive  changes  in  surface  profile  to  achieve  uniform 
contact  with  the  soil. 


5.6.1  Self-Boring  Operation 

The  self  boring  attachment,  a  schematic  drawing  of  which  can  be 
seen  in  Figure  32  is  based  on  the  system  developed  at  Cambridge  University 
for  placement  of  in-situ  stress  gauges  and  pressuremeter  testing  apparatus. 

The  process  is  a  highly  cor  -oiled  and  steady  removal  of  soil  in  a  manner  similar 
to  a  tunnelling  machine.  By  careful  control  of  certain  parameters  in  the 
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apparatus  and  in  the  drilling  procedures,  it  is  possible  to  obtain  extremely 
small  zones  of  disturbance  in  the  soil  and  very  little  modification  from  the 
undisturbed  stress  state. 

A  thin-walled  circular  blade  or  cutting  tip  is  mounted  on  the  end 
of  the  downhole  assembly  the  outside  profile  of  which  is  a  smooth  and  continuous 
cylinder.  Axial  thrust  is  applied  to  advance  the  blade  into  the  soil  and  the 
isolated  central  core  is  broken  up  by  a  rotating  blade  located  behind  the  tip. 

The  total  thrust  required  to  advance  the  apparatus  below  the  borehole 
casing  will  be  comprised  of  the  following  components: 

(a)  A  net  cutting  force  developed  over  the  tip  of  the  rim  cutting  blade 

(b)  Adhesion  or  friction  forces  developed  over  the  surface  of  the  in-situ 

equipment  by  contact  with  the  soil  under  lateral  stress. 

The  maximum  friction  or  adhesion  forces  will  be  the  greatest  component 
of  the  two  above  (see  section  5.7.7);  these  forces  depend  heavily  upon  the  type 
of  soil,  the  depth  of  the  selected  measurement  location  and,  of  course,  on  the 
distance  of  placement  below  the  borehole  casing.  The  rim  blade  of  the  self 
boring  equipment  is  conic  within  the  cutting  head  leading  past  the  breaker  blade. 
The  blade  is  designed  such  that  the  stresses  in  the  soil  ahead  of  the  blade  are 
changed  as  little  as  possible.  Since  the  amount  of  load  applied  to  advance 
the  cutter  is  critical,  it  becomes  convenient  to  include  some  form  of  monitoring 
and  control  for  the  force  developed  at  the  cutting  head. 
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The  soil  cuttings  are  conveyed  by  the  advancing  core  pressure  to 
a  throated  section  at  which  point  they  are  picked  up  by  a  helical  auger.  The 
auger  is  housed  in  a  concentric  bore  which  runs  through  the  entire  length 
of  the  apparatus  and  opens  into  a  catch  tank  located  at  the  top  of  the  downhole 
assembly.  Hence  the  auger  drive  is  used  both  to  power  the  core  breaker  blade 
and  to  lift  the  soil  cuttings  through  the  apparatus.  It  is  conceivable  that 
it  could  also  be  used  to  transmit  the  axial  thrust  required  on  the  rim  blade 
during  the  placement  operation.  In  Figure  32  the  drive  rod  is  shown  socketed 
into  a  drive  coupling  within  the  cutting  head.  This  arrangement  allows  for 
the  auger  to  be  withdrawn  after  self-boring  leaving  an  uninterupted  central 
bore  through  the  downhole  assembly.  This  feature  assists  not  only  in  helping 
to  remedy  problems  with  jamming  or  fouling  of  the  placement  equipment,  but 
as  discussed  in  section  5.8.1,  it  permits  the  use  of  a  very  convenient  means 
of  instrumenting  the  rotation  of  the  coupling  cylinder. 


For  clayey  soils,  the  auger  may  not  be  the  most  appropriate  means 
for  lifting  the  cuttings.  Without  lubricating  and  flushing  with  water, 
it  is  highly  possible  that  these  soils  could  gather  in  coagulated  masses 
with  sufficient  adhesion  and  strength  to  foul  the  auger.  To  overcome  this 
problem  conventional  flush  boring  methods  can  easily  be  employed;  indeed, 
with  the  detachable  auger  system  already  described,  it  would  be  possible 
to  change  boring  techniques  even  during  a  particular  placement.  Flushing 
water  could  be  pumped  through  the  centre  of  the  auger  drive  rod  and  the 
slurry  returned  via  the  annular  space  between  the  rod  and  the  bore  of  the 
coupling  cylinder.  This  space  would  have  to  be  reduced  in  order  to  raise 
the  fluid  velocity  but,  again,  the  means  for  achieving  this  could  be 
incorporated  onto  the  hollow  drive  rod  assembly  to  allow  simple  conversion 
of  the  downhole  equipment  to  water  flush  boring. 
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The  self-boring  operation  is  likely  to  be  the  longest  and  most 
in  need  of  careful  control  of  the  complete  testing  procedure.  It  would  be 
worth  attempting  to  make  some  assessment  of  the  likelihood  of  an  acceptable 
placement  and  of  the  soil  quality  and  uniformity,  whilst  the  self  boring  process 
is  in  progress.  Hence  it  is  suggested  that  at  least  for  the  initial  trials 
of  the  field  equipment,  the  principal  parameters  associated  with  the  placement 
equipment  be  instrumented  and  provision  be  made  so  that  they  can  be  monitored 
directly  on  site.  If,  for  instance,  a  simple  stripchart  type  record  is  produced 
during  self-boring  operations  showing  penetration,  tip  thrust,  auger  drive  speed 
and  torque,  it  is  probable  that  an  immediate  on-site  impression  could  be  gained 
concerning  the  progress  of  the  placement.  Under  certain  conditions  these  data 
might  be  used  to  test  the  suitability  of  a  given  location  or  to  determine  whether 
or  not  to  abort  a  particular  placement. 

5.6.2  Coupling  to  the  Soil 

Having  reached  the  desired  measurement  location  in  the  virgin  soil, 
the  final  phase  of  the  placement  process  occurs:  that  of  jacking-out  or 
locking-on  to  the  borehole  wall.  Because  of  the  sensitivity  of  the  test  to 
the  state  of  soil  immediately  in  contact  with  the  coupling  cylinder,  the 
apparatus  must  replace  as  closely  as  possible  the  in-situ  or  undisturbed  lateral 
stress  condition.  In  this  respect,  two  factors  are  of  prime  importance: 

(a)  The  cylinder,  should  achieve  upon  jacking  out  a  contact  pressure 
with  the  soil  as  uniform  as  possible.  Thus,  the  shoe  should 
be  able  to  adhere  closely  to  the  uniform  cylindrical  cavity 
produced  by  the  self  boring  operation  and,  if  necessary,  should 
also  be  able  to  distort  slightly  in  order  to  accommodate  small 
anomalies  or  gradients  in  the  properties  of  the  surrounding  soil. 
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(b)  A  means  must  exist  to  determine  the  undisturbed  lateral  stress 

in  the  soil,  so  that  the  correct  jacking  pressures  can  be  applied. 


To  achieve  uniform  contact  stress  with  the  soil,  the  surface  of 
the  cylinder  is  composed  of  a  large  number  of  shoe  segments.  These 
segments  act  entirely  independently  and  are  free  to  both  tilt  and  displace 
radially.  The  large  number  of  shoes  employed  ensures  that  the  cylinder 
can,  if  necessary,  maintain  a  close  approximation  to  the  cylindrical 
profile,  whilst  the  independent  movement  of  the  segments  allows  the  system 
some  tolerance  to  variations  in  the  soil  properties.  For  the  undisturbed 
lateral  stress  in  soil,  the  following  is  a  list  of  viable  techniques: 


(a)  By  using  a  plane-strain  or  other  type  of  horizontal 
penetrometer  probe. 

(b)  By  monitoring  the  total  stress  and  the  pore  water  pressure, 
n,.d  subsequent  calculation  based  on  these  data  and  correction 
factors 

(c)  By  calculation  or  estimation  based  on  experience  and/or  ..ite 
data. 

(d)  By  observation  of  the  pressure-displacement  curve  during 
radial  expansion  of  the  shoe. 


The  penetrometer  test,  conducted  using  horizontally  action  probes 
located  behind  the  self  boring  head,  would  load  the  soil  entirely  through 
the  range  of  lateral  stress  from  the  initial  stress  up  to  yield.  The  resultant 
stress-strain  characterist’- c  would  provide  the  required  soil  data. 


The  stress  gauge  is  a  more  simple  expedient  than  the  penetrometer 
and  requires  no  intricate  mechanics  and  little  or  no  control  procedures. 

The  measurement  is  subject  to  errors,  however,  principally  arising  from 
disturbance  or  strpas: relief  due  to  placement  and  to  the  well-known  phenonomen 
of  soil  arching  eve*:  the  gauge.  With  experience,  the  output  of  the  stress 
gauges  could  no  doubt  be  used  to  achieve  quite  reasonable  results  for  a 
majority  of  cases. 

Without  any  measurements,  it  is  felt  that  little  better  than  20% 
accuracy  could  be  achieved  in  the  determination  of  the  undisturbed  lateral 
stress.  This,  it  should  be  emphasized,  assumes  a  uniform  uncomplicated 
site  with  well  known  geology.  It  may  occur  that  under  certain  conditions 
this  would  be  a  perfectly  adequate  method,  but  for  the  analysis  of  early 
data  from  the  BSD  test,  such  lack  of  certainty  can  only  hinder  the  rate 
at  which  familiarity  with  the  procedure  grows. 

As  a  further  indication  of  the  quality  or  state  of  coupling  to 
the  surrounding  soil,  it  is  proposed  that  displacements  of  the  expanding 
elements  of  the  coupling  cylinder  be  instrumented.  In  the  coupling  cylinder 
design  described  in  this  report,  expansion  is  achieved  by  means  of  rigid 
continuous  keys  moving  in  slots  in  the  surface  of  a  central  cylindrical 
core  (see  Fig  32) .  Displacement  of  these  keyways  is  instrumented  at 
both  ends  of  the  coupling  cylinder  giving  information  on  the  mean  displacement 
and  tilt  of  the  moving  elements.  Jacking  force  is  achieved  by  a  system  of 
hydraulic  pistons  and  cylinders  and  hence  the  mean  lateral  jacking  stress 
can  easily  be  determined  from  the  pressure  in  the  hydraulic  system. 
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As  described  in  Section  4.5,  it  is  also  possible  to  choose  an 
optimum  point  at  which  the  mean  lateral  stress  on  the  coupling  cylinder 
is  close  to  the  original  undisturbed  stress;  this  is  done  by  monitoring 
the  pressure/displacement  function  during  jacking  out.  Should  there  be 
an  inclusion  or  gross  discontinuity  in  the  soil  which  results  in  correspond¬ 
ingly  erratic  shoe  displacements,  the  proposed  monitoring  system  will  readily 
indicate  this  and  allow  a  quick  decision  to  be  made  concerning  the  success 
or  otherwise  of  a  particular  placement.  This  method  has  the  advantage 
of  providing  the  same  data  during  expansion  of  the  shoe  that  existing  self¬ 
boring  pressuremeters  produce.  The  procedure  hence  gives  more  information 
about  the  soil  and  creates  at  the  same  time  one  more  link  between  the  new 
data  and  other  results  which  already  count  with  years  of  expertise. 
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5.7  SYSTEMS  DESIGN  SPECIFICATION 


5.7.1  Size  Of  The  Coupling  Cylinder 


Probably  one  of  the  most  fundamental  questions  that  could  be 
asked  concerning  the  design  of  a  practical  BSD  is  what  are  the  factors 
which  govern  the  physical  size  of  the  coupling  mechanism.  Clearly  the 
dimensions  are  limited  at  the  lower  extremes  by  the  particle  size  of  the 
soils  in  which  the  device  is  expected  to  work.  For  the  upper  bound,  two 
significant  limitations  arise:  these  are  the  cost  of  the  preliminary  boring 
and  casing  operations,  and  the  increasing  power  requirements  of  the  system. 

The  dimensions  finally  chosen  for  the  coupling  cylinder  were  arrived  at  by  the 
reasoning  that  the  larger  the  device,  the  more  representative  the  test  and 
the  greater  the  range  of  materials  in  which  it  could  operate.  This, 
however,  was  tempered  by  the  correspondingly  higher  power  requirements  and 
the  fact  that  use  of  the  apparatus  might  be  limited  if  non-standard  or 
prohibitively  expensive  borings  were  required.  Thus  it  was  decided  that 
design  specifications  should  be  based  on  the  largest  readily  available 
cased  borings  that  could  be  produced  without  using  specialised  equipment. 

This  requirement  seems  to  be  satisfied  by  200  mm  (eight  inch)  boreholes  and 
hence  the  design  specifications  and  the  design  of  the  prototype  coupling 
mechanism  are  based  on  this  size  of  boring. 


5.7.2  Operational  Specifications 

The  specifications  for  the  power  supplies  to  the  equipment  and 
for  the  coupling  cylinder  design,  are  based  on  the  following  worst  case 
assumptions  for  the  operational  conditions: 


(a)  The  maximum  depth  of  operation  of  the  equipment  is  50  metres. 


(b)  The  water  table  is  not  lower  than  25  metres. 


(c)  The  maximum  coupling  cylinder  drive  torque  is  determined  by  that 


required  to  fail  the  soil  in  shear. 


(d)  The  coupling  system  is  to  be  operated  at  contact  pressures  close 


to  the  original  undisturbed  lateral  stress. 


For  the  parameters  of  the  different  soil  types  in  which  the  BSD 


has  potential  for  testing,  the  spread  of  values  is  large.  A  number  of 


practical  disadvantages  weigh  against  a  totally  worst  case  choice  for  the 


values  of  these  parameters  and  it  was  decided  that  calculations  should  be 


based  on  nominal  values  only.  This  is  justified  to  a  certain  extent  in 


that,  should  a  worse  case  or  extreme  value  situation  arise,  the  effect  on 


the  testing  program  is  not  catastrophic.  Testing  can  still  proceed  although 


some  limitation  will  occur  on  the  depth  to  which  the  apparatus  can  be 


operated  and  the  maximum  strain  at  which  the  test  can  be  run.  Furthermore, 


since  the  dynamic  range  is  finite,  design  for  an  extreme  worse  case  decreases 


the  resolution  of  the  system  thus  losing  some  of  its  small  strain  capabilities 


5.7.3  Calculation  of  Maximum  Drive  Toraue 


The  effective  lateral  stress  in  the  soil  before  placement  of  the 


coupling  cylinder  is: 


o '  =  gK  ( (p  -  p  }  d  +  n  d  ) 

h  o  s  w  w  w 


where 


7 '  =  effective  lateral  stress  (N/m^) 

n 

2 

g  =  gravitational  constant  (m/sec  ) 


K  =  lateral  earth  pressure  coefficient 

°  3  3 

Ps  =  mean  density  of  soil  overburden  (10  Kg/m  ) 


p  =  density  of  water 
w 


d  --  depth  of  measurement  from  the  soil  surface 


d^  =  depth  of  measurement  from  the  water  table 


Neglecting  soil  adhesion  to  the  shoe  and  assuming  that  the 
friction  coefficient  at  the  shoe-soil  interface  is  equal  to  the  angle 
of  internal  friction,  the  torque  required  to  yield  the  soil  approximately 
is  given  by: 


where  : 


T  =  2flr  ha’  tan  <t> 
m  oh 

T  =  peak  toraue  applied  to  the  coupling  cylinder  (Nm) 
m 

r  =  radius  of  the  coupling  cylinder 
o 

h  =  height  of  the  coupling  cylinder 
<{>  =  internal  friction  angle  of  the  soil 


Assuming  the  following  values  for  the  variables  in  the  above 

3  3 

equations:  tan  <J>  =  0.5,  K  =1.0,  p  =  1.8  x  10  Kg/m  ,  r  =  0.1  m  and  h  =  0.5  m, 

OS  o 

4 

a  maximum  design  torque  T  =  1x10  Nm  results. 

m 


5.7.4  Calculation  Of  Maximum  Jacking  Pressure 

For  the  lateral  jacking  pressure,  the  mean  lateral  stress  at  the 

surface  of  the  coupling  cylinder  must  equal  the  total  applied  stress.  Hence, 

3  3 

assuming  a  saturated  overburden  density  of  2.0  x  10  Kg/m  and  Kq  =  1.0, 
this  gives  a  maximum  jacking  pressure  of: 


a  =  lxlO6  N/m^ 


5.7.5  Calculation  of  Maximum  Range  of  Rotation 

Using  the  hyperbolic  soil  model  developed  in  Section  4.3  ,o  describe 
the  non-linear  behaviour  of  the  soil  at  high  strains,  an  approximate  calculation 
for  the  range  of  rotations  occuring  at  the  coupling  cylinder  can  be  derived. 

The  relationship  between  shear  stress  at  the  cylinder  surface  and  rotation  of 


the  cylinder  has  been  found  to  be: 
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where 


t  = 
o 


T  = 
m 


"o“ 

0  = 
o 


shear  stress  at  the  cylinder/soil  Interface 
yield  shear  stress  of  the  hyperbolic  soil 
elastic  shear  modulus  of  the  soil 
rotation  of  the  cylinder 


For  a  maximum  value  of  shear  stress  t  =  0.8  t  ,  the  rotation  is 

o  m 


0 


om  2)i 


In  5 


Substituting  for  t  and  p  gives 
m  o 


2  2 

T  3.2mr  hp  v 
m  oss 


0 


om 


In  s 


Hence,  the  apparent  rotational  stiffness  o?  the  in-situ  coupling 


cylinder  at  t  =  0.8  x  is: 

o  m 


A  _n  1  *\ 

- —  =  5.62  x  lO  v  Nm/rad 
”  __  s 

on 


Given  the  maximum  design  torque  derived  earlier  of  T  =  1  x  1(j  Nm  and 

m 

assuming  a  range  of  shear  wave  velocities  100  <  vg  <  5O0m/sec,  the  corresponding 
range  of  maximum  rotation  is: 

7.12  x  10~4  <0  <  1.78  x  lO-2  rad 

mo  — 


or 


0.04  <0  <1.0  degrees 

mo  ’ 
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5.7.6  Calculation  of  Drive  Power 

Assuming  a  sinusoidal  drive  function  and  a  purely  elastic  load,  the  peak 
mechanical  power  required  to  drive  the  coupling  cylinder  is: 


p  =  itt  e  f 

m  in  on* 
where  f  =  frequency  (hertz) 

For  the  range  of  maximum  rotations  given  above  the  peak  drive  power  is: 


22.5  <  P  <  560  Watts/Hertz 
m 


For  the  power  absorbed  in  strain  energy  in  the  BSD  apparatus, 
the  major  component  is  due  to  the  twist  in  the  torsion  coupling  tube  in 
the  single  cylinder  system.  The  torsional  stiffness  of  a  hollow  tube  is 
given  by: 


4  4 

"Vr2  'rl} 


where: 


Vr2 


=  torsional  stiffness  of  the  tube 


=  inside  and  outside  radi  of  the  tube 


=  shear  modulus  of  the  tube  material 


=  length  of  the  tube 


For  a  sinusoidal  torque  function  applied  to  the  torsion  tube,  the 
peak  power  exchanged  between  the  torsion  tube  and  the  actuator  (i.e. :  non 


contributory  or  reactive  power)  is: 


9 


where  f  =  frequency  of  operation  of  the  torsional  drive  (Hertz) 
o 

Assuming  a  torsion  tube  of  150aa  outside  diaaetex,  10m  vail 

lO  ° 

thickness  (34.5  kg/m  tubing)  with  a  nominal  of  SxlO  N/a  ,  the  aaxiaun 

4 

reactive  power  (i.e.  at  the  aaxlaua  drive  torque  T  =  lO  Mi)  due  to  the 

X£ 

torsion  tube  coupling  is: 

P  =  131  vatts/aetre  hertz  (0.24  h.p./oetre  hertz) 
n 

Hie  total  theoretical  energy  required  assuming  the  drive  systea 
is  incapable  of  recovering  the  reactive  energy  stored  in  the  tube  is: 

E  =  P  Joules/ cyclj  metre 

as 

4 

which,  at  the  systea  oaxiaua  torque  of  =  1x10  Mi  and  at  the  aaxiaua 

3 

depth  of  operation  of  50  aetres  is  9.0x10  Joules/cycle. 

It  oust  be  reacnbered  that  this  figure  is  an  upper  bound  and  that 
only  a  mall  number  of  cycles  involving  energies  approaching  this  figure 
are  required  in  the  BSD  test.  If,  as  suggested  in  an  earlier  section,  an 
acooulator  storage  systea  is  used  to  provide  hydraulic  power,  then  the 
power  supply  systea  will  not  have  to  be  generating  and  damping  the  naxicua 
required  power  Although  efficiencies  will  be  low,  the  total  energy  required 
for  the  test  will  not  be  inordinately  high  and  could,  it  is  felt  be  stored 
in  a  high  pressure  of  gas  bladder  type. 

~.i  the  torsional  coupling  is  operated  at  resonance,  the  power 
required  will  effectively  be  divided  by  Q,  the  resonant  signification 
of  the  systea.  Under  these  conditions  the  torsion  tribe  would  be  thin-walled 
and  *n  approximate  expression  for  the  resonant  frequency  of  the  systea  is: 


where 


f  =  resonant  frequency  of  the  coupling  in  Hertz 
r^,  r  =  mean  radius  and  wall  thickness  of  the  coupling  tube 

I  =  lumped  moment  of  inertia  of  the  coupling  tube  and  drive  system 

The  above  analysis  assumes  a  much  greater  stiffness  for  the  in-situ 

coupling  cylinder  than  for  the  torsion  tube,  he  ice  the  tube  is  effectively 

fixed  at  one  end  with  a  lumped  mass  attached  to  the  free  end.  To  serve  as  an 

example,  a  torsion  tube  of  2.5mra  wall  thickness,  and  a  lumped  moment  of  inertia 
2 

I  =  10  Kg  m  ,  the  resonant  frequency  of  operation  at  L  =  50  metres  is: 

S  t  * 

f  =  5.18  Hz 
o 

Hence,  at  resonance,  an  approximate  expression  for  the  peak  reactive 
power  due  to  the  coupling  is: 


T  2  I  L 

—  / - -T - 

Q  /\J  8tt  u  r  r  I 
*  T  m  t  s 


(18) 


where : 


P  =  peak  power  at  resonance  (watts) 
mo 

Q  =  resonant  torque  magnification  of  the  system. 

In  terms  of  the  linear  damping,  either  inherent  in  the  drive  system 


or  artificially  introduced  by  means  of  feedback,  a  simple  relationship  exists 
for  systems  with  low  damping: 


(19) 


where: 


£  =  coefficient  of  viscous  damping  in  the  overall  transfer 

function  of  the  torque  control  system  associated  with  the  peak 
due  to  resonance  of  the  torsional  coupling. 


By  controlling  certain  terms  of  the  feedback  signal  in  the  torque 
control  system  it  would  be  reasonable  to  expect  a  range  of  Q  values 
between: 


1  <  Q  <  25 


The  main  drawback  of  the  resonant  coupling  method  would  be  the 
fact  that  the  buildup  of  torque  at  the  output  of  the  system  would  depend 
upon  the  characteristic  damping.  If  a  constant  amplitude  sinusoidal 
command  signal  is  applied  to  the  system,  the  output  torque  will  build  up 
exponentially.  The  time  constant  of  this  response  in  terms  of  the  system 
damping  and  resonant  frequency  is: 


t 

c 


1 

2tt  f  t 

U 


(20) 


where : 

t  =  time  constant  of  the  exponential  envelope  of  the  response 

f  =  undamped  resonant  frequency  of  the  system 
u 
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5.7.7  Calculation  of  Placement  Forces 


As  mentioned  previously,  the  nett  thrust  applied  to  advance 
the  apparatus,  into  the  ground  during  the  self-boring  operation  is  comprised 
of  two  components: 

(a)  a  tip  thrust  due  to  the  action  of  the  cutting  head 

(b)  a  frictional  force  developed  on  the  walls  of  the  in-situ  apparatus 
due  to  contact  with  the  soil. 


In  order  to  disturb  the  material  ahead  of  the  cutter  as  little 
as  possible,  the  nett  tip  force  must  be  close  to  the  m-situ  vertical 


stress  times  the  sectional  area  of  the  apparatus,  i.e.: 


f  =  irr  g 
t  o^ 


(p  -  0  )  d  +  p 
s  w 


d 

w  w 


(21) 


where  f  =  tip  force 


thus  giving  a  maximum  tip  force  ffc=  2x10  N 


Taking  the  case  of  sand  as  providing  the  greatest  frictional 
forces  on  the  surface  of  the  in-situ  apparatus: 


f  =  2 'nr  h  c  ^an^ 
s  o  t  n 


(22) 


where  f  =  force  developed  by  frictional  contact  with  soil 


h  =  self-bored  depth 


6  =  friction  angle  at  the  soil-shoe  interface  governing  vertical  relative 
measurement 


Assuming  tan  6  =0.3  and  h^  =  i.O  metres,  the  component  of 


thrust  due  to  friction  between  the  in-situ  apparatus  and  the  surrounding  soil 


is: 


fs  =  1.4  x  10  N 


Hence,  the  approximate  total  thrust  required  for  a  -self-bored  depth  of 


one  metre  is  1.4  x  10  N. 
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at  this  point  that  it  is  necessary  or  indeed  relevant  to  begin  to  measure 
the  rotations  of  the  coupling  cylinder.  Unless  a  technique  such  as  this 
is  employed,  the  rotation  transducer  must  be  capable  of  handling  significantly 
larger  rotations  during  general  handling  and  placement  than  will  be 
realized  during  the  measurement  process.  Optical  transducers  in  general 
posseses  this  ability,  and  they  can  even  tolerate  the  unwanted  displacement 
as  an  arbitrary  offset  or  bias  to  the  measurement  range  in  which  they  are 
required  to  operate.  They  are,  however,  expensive  and,  for  the  very  high 
resolution  required  in  the  BSD  test,  they  are  also  bulky  and  fragile. 

An  appealing  solution,  which  is  in  fact  the  one  proposed  in  this  study 
is  to  couple  the  transducer  into  the  in-situ  apparatus  immediately  before  the 
dynamic  loading  history  is  to  begin.  By  this  method,  the  only  offsets 
introduced  will  be  those  inherent  in  the  transducer  itself  and  any  contribution 
that  may  arise  from  the  system  employed  to  couple  the  transducer  to  the 
apparatus.  The  transducer  would  include  a  self  locking  attachment  which 
would  radially  expand  to  grip  the  inside  of  the  central  bore  in  the  coupling 
cylinder.  Being  a  frictional  coupling,  the  device  would  require  no  initial 
orientation  or  location  and  it  would  simply  assume  a  zero  displacement  for 
the  cylinder  with  respect  to  the  reference  system,  at  the  point  of  locking 
on  to  the  coupling  cylinder.  The  displacement  reference  for  the  system  is 
derived  from  a  rigid  rod  coupled  to  an  anchored  frame  at  the  ground  surface. 

The  transducer  and  locking  system  are  mounted  directly  on  the  reference  rod 
and  are  lowered  into  the  in-situ  apparatus  following  the  removal  of  the  auger 
or  flush  boring  equipment. 


The  design  suggested  for  the  transducer  employs  a  series  of  strain 
gauged  cantilevers  as  displacement  sensors.  The  cantilevers  would  be 
arranged  as  radial  web  type  members  coupling  the  reference  rod  to  an  outer 
cylinder  (see  Fig.  37) .  This  outer  cylinder  is  attached  to  or  in  part  of 
the  locking  mechanism  and  hence  follows  the  rotations  of  the  coupling 
cylinder.  The  required  gain  of  the  transducer,  determined  by  the  strain 
in  the  cantilever  surfaces  for  a  given  rotation,  can  be  achieved  in  a  given 
design  by  varying  the  length  of  the  cantilevers  and  their  radial  position. 

With  the  method  proposed  here,  the  sensitivity  cf  the  system 
could  be  changed  very  easily  by  replacement  of  the  rotation  transducer 
assembly.  Alternatively  it  would  be  possible  to  employ  transducers  of 
different  sensitivites  in  parallel  to  extend  the  dynamic  range  of  measurement, 
although  it  should  be  remembered  that  the  higher  sensitivity  device  would 
have  to  be  decoupled  beyond  a  certain  amplitude  of  rotation. 


5.8.2  Measurement  of  the  Coupling  Cylinder  Drive  Torque 

Here  again,  a  strain  gauge  transducer  is  proposed.  The  most  suitable 
location  of  the  gauges  would  depend  on  the  BSD  configuration,  i.e.  the  single 
or  double  cylinder  device,  but  in  any  case  the  transducer  would  be  configured 
on  the  shaft  or  tube  driving  the  coupling  cylinders  on  a  uniform  or  gauge 
section  close  to  the  cylinder.  Since  the  dynamic  range  required  at  this 
measurement  is  more  a  function  of  the  depth  of  operation  of  the  equipment. 


it  is  likely  that  adequate  dynamic  ra..^e  will  be  provided  at  any  given 


measurement  location  by  a  single  transducer  system.  For  the  generally 
increasing  drive  torque  required  with  increased  depth  of  the  test,  it 
would  be  a  simple  matter  to  change  transducers  between  tests.  This  would 
not  significantly  inconvenience  the  testing  procedure  since  the  entire 


BSD  apparatus  is  extracted  fron  the  boring  between  measurements  in  order 
that  the  cased  hole  can  be  achieved. 


5.8.3  General 

Much  of  the  remaining  instrumentation  of  the  BSD  apparatus,  although 
operating  in  what  may  be  described  as  harsh  environmental  conditions,  is  not 
regarded  as  particularly  critical.  Many  of  the  measurements  can  be  achieved 
using  standard  process  control  or  geotechnical  instrumentation  equipment  and 
therefore  a  discussion  of  the  methods  is  not  attempted  here. 


6.0  COUPLING  CYLINDER  DESIGN 


6 . 1  INTRODUCTION 


The  design  of  the  coupling  mechanism  with  the  corresponding 
investigation  into  the  practical  problems  of  operating  and  handling  the 
proposed  device  in  the  field,  has  proved  an  invalaaf-j  exercise.  Much 
discussion  has  taken  place  on  matters  such  as  the  ini • uence  of  soil 
properties,  drilling  procedures  and  instrumentation  problems;  the  findings 
of  the  work  have  often  resulted  in  modifications  to  the  course  of  the 
project  development.  Equally,  the  findings  of  the  analytical  work  have 
had  direct  and  significant  influences  in  the  hardware  design. 

The  following  section  presents  a  discussion  of  these  developments 
and  their  influence  on  the  prototype  hardware  design.  Later,  the  overall 
design  considerations  for  the  coupling  mechanism  are  listed  and  discussed, 
and  the  final  design  of  the  coupling  mechanism  is  described.  Aspects  of 
the  design  affecting  performance  and  manufacture  are  then  dealt  with  and 
finally,  a  description  of  the  laboratory  testing  apparatus  is  given. 

6.2  DESIGN  CONSIDERATION 


As  originally  conceived,  the  BSD  operated  in  a  constrained  system 
where  two  coupling  cylinders  provided  complementary  reaction  torques  for  each 
other.  The  region  in  which  high  strains  were  to  be  induced  was  the  space  sur¬ 
rounding  the  borehole  be -ween  the  two  cylinders.  Hence  a  major  portion  of  the 
coupling  surface  was  distant  from  the  zone  of  high  strain  and  the  effects  of  both 


non-uniformity  of  contact  and  of  hysterysis  or  backlash  in  the  coupling  cylinders 
were  controllal «'  in  that  they  were  dependent  on  the  spacing  of  the  cylinders. 


However  it  was  found  that  for  most  realistic  configurations  of  the 
double  cylinder  BSD,  t  0  strains  predominate  over  the  f  strains.  Hence, 
as  discussed  earlier,  the  thrust  of  the  study  of  the  3SD  test  was  directed 
towards  the  r  strain  and  its  effect  on  the  design  specifications  for  the  apparatus. 


For  inducing  strain  the  dominant  plane,  the  characteristics  of 

the  coupling  mechanism  become  critical:  the  zone  of  high  strain  in  the  soil  is 
small  in  relation  to  the  coupling  cylinder  and  the  maximum  value  occurs  at  the 


interface  between  the  soil  and  the  cylinder.  The  state  of  stress  or  strain  in 


the  surrounding  soil  can  only  be  inferred  from  measurement  of  applied  torque  and 


the  resultant  rotation  of  the  cylinder.  Hence  the  limit  of  resolution  or  dynamic 


range  of  the  measurement  is  directly  linked  to  the  mechanical  characteristics  of  the 
coupling  cylinder.  Of  particular  importance  is  the  degree  of  backlash  or  hysteresis 


exhibited  by  the  system  under  the  reversing  drive  torque. 


Backlash  or  hysteretiu  effects  in  the  torque/displacement  characteristics 
of  the  cylinder  are  likely  to  arise  mainly  from  the  mechanism  which  allows 
radial  expansion  of  the  cylinder.  Some  clearance  must  exist  to  permit  this 
movement  and,  although  it  is  proposed  that  a  viscous  lubricant  be  used  between 
the  moving  surfaces,  it  must  be  accepted  that  some  non-linear  behavior  will 
occur.  In  the  section  on  design  specifications,  the  range  of  maximum  rotations 
of  the  coupling  cylinder  under  the  design  maximum  torque  were  found  to  be: 

-4  -2 

7.12  x  10  <  C  <  1.78  x  10  (rad) 

O 
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The  radial  movement  which  allows  expansion  of  the  coupling  cylinder 
occurs  near  the  surface  and  a  tolerable  lateral  movement  or  clearance  in  this 
system  can  be  determined  from  an  equivalent  rotation.  The  tangential  displace¬ 
ments  u  at  the  surface  of  the  coupling  cylinder  corresponding  to  the  above 
rotations  are 

0.07  <  u  <  1.8  (mm) 


The  required  resolution  of  measurement  must  be  applied  to  these  displace¬ 
ments  in  order  to  determine  the  allowable  lateral  movements  due  to  clearance  in  the 
coupling  cylinder  jacking  system.  The  upper  bound  to  the  maximum  displacements 
requires  a  total  movement  of  less  than  0.02  mm  for  1%  resolution  which  is  not 
a  particularly  difficult  accuracy  to  achieve  in  the  manufacture  of  the  coupling 
cylinder,  but  nonetheless  falls  within  the  realm  of  close- tolerance  engineering. 
Clearly,  for  the  lower  bound,  it  would  be  highly  impractical  to  expect  to  be 
able  to  achieve  so  small  a  backlash  from  purely  machined  and  fitted  components. 

Thus  the  coupling  cylinder  will  rely  heavily  on  the  ability  of  the  lubricant  film 
on  the  jacking  system  bearing  surfaces  to  take  up  the  manufacturing  clearances. 


The  critical  effects  of  disturbance  due  to  placement  have  been 
mentioned  earlier.  The  T/0  relation  derived  from  the  apparatus  will  depend 
heavily  on  the  state  of  the  soil  immediately  surrounding  the  coupling  cylinder, 
hence  the  magnitude  and  the  extent  of  the  disturbance  are  important  factors. 


For  the  lacking-on  part  of  the  placement  procedure,  the  major  contribution 
tc  disturbance  is  likely  arise  frcm  non-uniformities  on  the  contact  stresses 
at  the  soil/cylinder  interface.  The  cylinder  profile  must  be  as  regular  and  con¬ 
tinuous  as  possible  in  order  to  minimize  the  effects  of  local  yielding  and 
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weakening.  Such  effects  will  certainly  limit  the  degree  of  coupling  attainable 
with  the  soil  and,  being  in  the  critical  measurement  zone,  could  result  in 
anamolous  or  unrepresentative  data.  To  achieve  the  desired  highly  uniform 
contact  with  the  soil,  the  mechanism  which  allows  radial  expansion  of  the 
cylinder  must  be  comprised  of  a  relatively  large  number  of  elements.  Thus 
the  departure  from  the  ideal  circular  section  with  expansion  will  be  minimized 
at  least  to  the  extent  that  manufacturing  costs  allow. 

6.2.2  torsional  Stiffness 

Finally,  the  discussion  in  this  section  on  fundamental  design  considerations 
fo.  the  coupling  cylinder,  is  the  question  of  torsional  stiffness  of  the  device. 

In  the  process  of  coupling  shear  stress  into  the  soil,  the  coupling  cylinder 
will  itself  twist  under  the  applied  torsion.  Hence  a  rotation  along  the  cylinder 
will  occur  which  will  in  effect  reduce  the  nett  stress  applied  to  the  soil.  If 
rot  held  to  admissable  limits,  the  effect  would  certainly  degrade  the  measurements, 
particularly  since  it  would  be  difficult  bo  compensate  for  it  in  the  analysis  of 
the  field  data. 


In  order  to  assess  the  influence  of  twist  in  the  cylinder,  the  ratio 
of  twist  in  a  non-rigid  cylinder  with  one  end  fixed,  to  the  rotation  of  a 
perfectly  rigid  cylinder  in  an  elastic  space  (see  section  3  )  is  formulated. 
For  a  perfectly  rigid  cylinder  bonded  to  the  walls  of  a  continuous  hole  in  an 
Clastic  space,  the  ratio  of  rotation  to  applied  torque  is  given  by: 


0 

o 

T 


(1) 
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o 


133 


'll  ^  LtffhUt  2\l4^MWj4X& , 


=  1.44 


a 

e 

o 

which  is  felt  to  be  a  reasonable  starting  point  in  view  of  both  the  highly 
conservative  nature  of  the  calculations  and  of  the  worse  case  assumptions 
in  the  values  used  for  the  soil  properties. 

6.2.3  Summary  Of  Fundamental  Design  Considerations 

The  following  list  summarizes  the  fundamental  design  considerations 
discussed  above  together  with  additional  practical  considerations  relevant  to 
the  design  of  a  system  suitable  for  operation  in  the  field. 

(a)  Coupling  to  the  soil:  the  coupling  cylinder  must  be  able  to 
expand  radially  under  controlled  lateral  stress.  This  feature 
allows  the  system  to  accommodate  soil  compaction  and  deformation 
under  the  combined  static  lateral  and  dynamic  shear  stresses. 

(b)  Torsional  stiffness:  In  order  to  comply  with  the  analytical  model 
the  cylinder  should  be  torsionally  rigid. 

(c)  Load  characteristics :  The  design  should  yield  minimal  interaction 
between  the  dynamic  shear  stress  and  the  static  lateral  stress 
applied  to  the  soil. 

(d)  Coupling  characteristics:  non-linear  hysteretic  or  backlash  effects 
in  the  coupling  cylinder  load  at  least  over  the  range  of  operation 
of  the  test,  should  be  small  in  relation  to  the  rotations 
inducing  stra-'-  in  the  soil. 

(e)  Loading  considerations:  the  coupling  cylinder  must  withstand  the  load- 
applied  during  and  resulting  from  the  self  boring  operation.  Principally, 


6.3  DESIGN  GF  THE  PROTOTYPE  COUPLING  CYLINDER 


The  object  of  the  hardware  design  and  construction  aspect  of  the 
project  was  to  produce  a  prototype  version  of  the  coupling  cylinder,  upon 
which  laboratory  tests  could  be  performed  to  investigate  the  more  basic  charac¬ 
teristics  of  the  mechanism.  Apart  from  the  inherently  enlightening  process 
of  design  and  manufacture,  the  exercise  was  aimed  at  producing  some  hardware 
which  hopefully  can  be  used  in  subsequent  field  trials.  Certain  features 
essentia]  to  the  field  version  of  the  coupling  cylinder  are  not  incorporated 
in  the  first  prototype.  These  features,  it  was  felt,  did  not  present  a  significant 
challenge  in  their  design  or  manufacture  and  were  better  relegated  to  subsequent 
phases  of  the  BSD  development.  The  prototype  described  here,  however,  has  been 
designed  with  many  of  these  considerations  in  mind  and  should  require  no  major 
modifications  to  incorporate  any  additional  features  required  by  the  field  version. 

The  design  finally  chosen  for  the  coupling  cylinder  utilises  a 
series  of  beams  keyed  into  the  surface  of  a  hollow  cylinder  (see  Fig.  37).  The 
keys  are  rectangular  in  section  and  are  ground  to  a  precision  fit  with  mating 
grooves  in  the  surface  of  a  thick  walled  cylinder.  Lateral  jacking  forces  are 
derived  from  a  system  of  single-acting  hydraulic  pistons  and  cylinders  located 
within  the  body  of  the  keys  and  arranged  at  regular  intervals  along  their  length. 
Shear  stress  is  coupled  from  the  keys  to  the  cylinder/soil  interface  by  means 
of  shoes  which  are  attached  to  the  outer  faces  of  the  keys.  These  shoes  are 
curved  in  outside  profile  to  achieve  uniform  contact  with  a  circular  opening, 
and  are  provided  with  a  fine  splined  surface  in  order  to  enhance  frictional 


coupling  capabilities. 


The  concept  of  this  design  evolved  from  the  reasoning  that  the 
maximum  strength  and  torsional  stiffness  would  be  realized  from  a  mechanism 
that  approaches  as  closely  as  possible  the  optimum  configuration  of  a  thick 
walled  hollow  cylinder.  Apart  from  the  locations  of  the  jacking  pistons,  the 
section  of  the  design  described  above  is  a  continuous,  virtually  solid  thick- 
walled  tube.  Shear  and  compressive  forces  are  applied  to  the  components  of  the 
cylinder  uniformly  and  over  large  contact  areas.  Hence  the  design  attempts 
to  minimize  stress  concentrations  with  their  inherent  penalty  of  inefficiency 
of  load  rating  and  reduced  service  life. 

Torsion  is  applied  by  means  of  another  set  of  keys  located  at  one 
end  of  the  grooved  cylinder.  These  keys  are  similar  in  section  to  the  laterally 
moving  keys,  but  they  are  force-fitted  into  the  cylinder  and  are  utlised  only 
in  coupling  load  to  the  walls  of  the  cylinder  grooves.  The  keys  are  tapered 
in  section  leading  axially  out  of  the  cylinder,  the  taper  forming  part  of  an 
anti-backlash  high  torque  coupling  which  employs  wedges  to  eliminate  clearances. 
These  wedges  normally  bear  on  the  walls  of  the  keys  protruding  from  the  cylinder 
and  are  normally  preloaded  by  means  of  sprung  plates. 

For  retracing  the  laterally  moving  keys,  an  arrangement  similar 
to  the  integral  jacking  pistons  mentioned  earlier  is  proposed  (see  Fig  37.) 

This  would  again  be  a  single-acting  piston  and  cylinder  arrangement  machined 
into  the  fixed  keys  and  arranged  to  apply  load  laterally  inwards.  The  scheme 
suggested  for  the  final  design  utilizes  a  domed  headed  piston  acting  on 
a  heel  extension  of  the  laterally  moving  keys.  Provision  has  been  made  in  the 
design  of  tne  prototype  c<-"'pling  cylinder  to  incorporate  this  system  and  a  means 
of  measuring  the  positions  of  the  ends  of  the  moving  keys  into  the  apparatus. 


A  clearance  is  deliberately  provided  between  the  laterally  moving 
keys  and  the  bottom  of  the  grooves  in  the  central  cylinder.  This  space, 
being  relatively  distant  from  the  surface  of  the  cylinder  does  not 
significantly  degrade  the  torsional  stiffness;  it  does,  on  the  other  hand, 
serve  a  number  of  practical  uses: 

(a)  It  allows  a  loose  tolerance  on  the  depth  of  the  grooves  and  also 
allows  that  the  close  tolerance  on  the  groove  width  does  not  have 
to  be  maintained  to  the  bottom. 

(b)  The  bottom  corners  of  the  groove  can  be  rounded  to  minimize 
stress  concentrations . 

(c)  Transducer  cables  and  small  hydraulic  pipes  can  be  routed  through 
the  cylinder  and  afforded  good  protection. 

The  lateral  jacking  pistons  are  mushroom  shaped  structures 
protruding  from  the  bottom  of  the  cylinder  grooves.  The  necking  or  relief 
away  from  the  sealing  ring  on  the  piston  allows  the  key  to  tilt  whilst 
requiring  only  a  minimum  of  radial  clearance.  Five  pistons  are  arranged 
along  the  keys  to  distribute  the  jacking  load  and  each  unit  is  through-fed 
on  a  common  hydraulic  line.  For  the  laboratory  prototype  coupling 
cylinder,  no  provision  has  been  incorporated  to  isolate  these  pistons  from 
the  shear  stresses  that  would  result  from  any  loads  applied  to  the  keys  in 
the  direction  of  the  axis  of  the  coupling  cylinder.  These  loads,  as 
mentioned  in  the  section  on  design  considerations,  will  arise  from 
frictional  contact  of  the  shoe  segments  of  the  cylinder  surface  with  the 
surrounding  soil.  It  may  he  that  for  most  soil  conditions  encountered  in 
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testing,  these  shear  stresses  may  not  damage  the  piscons  and  that  no 
special  clamps  or  stops  are  required  to  absorb  the  load.  However,  if 
required,  it  is  suggested  that  for  the  field  design  the  fixed  keys  in  the 
top  of  the  cylinder  be  closely  butted  to  the  moving  keys  to  form  a  stop. 

Obviously  this  technique  cannot  also  be  used  to  absorb  loads  due  to 
extraction  of  the  apparatus,  however,  it  is  likely  that  upon  retraction  of 
the  coupling  keys,  the  frictional  contact  with  the  soil  will  be  considerably 
reduced.  Under  these  circumstances,  it  is  almost  certain  that  the  jacking 
pistons  will  easily  withstand  the  axial  loads  on  the  shoe  segments. 

6.4  ASSEMBLY  AMD  TESTING  OF  THE  PROTOTYPE  COUPLING  CYLINDER 

Assembly 

The  mechanical  functions  of  the  coupling  cylinder  are  threefold;  to  supply 
a  radial  pressure  to  the  soil  to  restore  the  lateral  in-situ  stresses,  to  key  into 
the  soil  so  that  torsional  loads  may  be  applied  and  then  to  transmit  those  loads. 

The  complete  set  of  design  drawings  are  included  in  Appendix  VI  together  with  a  set 
of  photographs  of  the  device  and  test  equipment.  The  reimposition  of  the  radial  soil 
stresses  and  the  keying  are  both  achieved  by  the  radial  expansion  of  the  twelve  shoes, 
whilst  the  torque  is  transmitted  by  imposing  a  small  rotation  about  the  vertical  axis. 
These  two  mechanical  actions  are  mutually  conflicting  in  the  configuration 
adopted.  Radial  movement  is  accommodated  by  the  provision  of  slides  which  must  have 
some  circumferential  clearance  to  avoid  binding  through  metal-to-metal  welding.  On 
the  other  hand,  the  anticipated  rotations  of  the  cylinder  when  applying  a  couple  to  the 
soil  would  be  of  the  ordej.  of  0.03  degrees,  corresponding  to  an  absolute  movement  at 
the  circumference  of  the  cylinder  of  about  0.05  mm,  comparable  with  the  clearance 
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necessary  to  ensure  satisfactory  slide  operation.  The  slides  therefore  had  to  be 
manufactured  with  minimal  clearance,  necessitating  careful  hand  fitting  and  scraping. 
Because  of  the  length  of  these  slides  the  accuracy  of  the  final  result  was  dependent 
upon  the  straightness  of  the  slots  and  of  the  slidi  jars.  Inevitably  each  bar  had 
to  be  selectively  fitted  into  a  specific  slot. 


The  workmanship  was  to  a  high  standard.  Only  two  bar  and  slot  pairs 
showed  any  circumferential  clearance  which  could  be  detected  by  hand  when  clean  and 
oil- free.  The  others  tended  to  bind  when  dry  but  moved  smoothly  and  progressively 
with  steady  resistance  when  lubricated  (there  qualitative  observations  were  made 
before  the  six  pistons  were  inserted) .  It  is  therefore  concluded  that  the  cylinder 
was  manufactured  to  a  standard  which  closely  matched  the  fine  tolerances  demanded 
by  the  design. 


The  fix  pistons  were  found  to  be  in  accurate  alignment  with  the  cylinders 
on  each  of  the  sliding  bars  (a  total  of  60)  and  no  problems  were  encountered  with 
their  operations.  Nevertheless,  it  is  now  realized  that  the  pistons  should  be  re¬ 
designed  because,  to  avoid  alignment  problems,  they  were  manufactured  0.08  mm  under¬ 
size  on  this  diameter  with  the  result  that  the  ’O'  ring  seals  would  be  intruded  at 
hydraulic  pressures  in  excess  of  35  MPa. 


Testinc 


The  objective  of  the  mechanical  testing  of  the  coupling  cylinder  was  to 
simulate  its  operation  under  laboratory  conditions  so  that  potential  problems  could 
be  identified.  The  tests  _re  conducted  using  a  dummy  borehole  fabricated  from  a 


length  of  203  mm  (8  inch)  bore  mild  steel  tube  with  support  flanges  welded  at  each  end. 
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In  the  initial  stages  of  testing,  the  tube  was  lined  with  wooden  strips 
to  reduce  its  internal  diameter  to  185  mm  and  to  provide  a  deformable  material  into 
which  the  shoes  could  key.  To  contact  the  wood  each  shoe  had  to  be  jacked  out 
through  a  mean  radial  distance  of  6  mm.  These  tests  proved  that  the  radial  loading 
system  worked  satisfactorily  and,  once  all  hydraulic  leaks  had  been  located,  it  was 
possible  to  hold  hydraulic  oil  pressures  of  14  MPa  in  the  shoe  rams  for  periods  of 
several  hours.  All  rams  functioned  smoothly  at  pressures  above  700  KPa  without  un¬ 
due  sticking. 

Torque  loading  was  applied  to  the  coupling  device  by  a  pair  of  hydraulic 
rams,  bolted  to  the  top  flange  of  the  dummy  borehole.  These  acted  tangentially  onto 
a  reaction  plate  bolted  to  the  cylinder.  To  ensure  that  the  torque  was  transmitted 
uniformly  to  each  shoe,  this  reaction  plate  was  positively  keyed  to  each  of  the  radial 
slide  slots  by  a  series  of  metal  wedges.  These  torque  rams  were  loaded  by  a  hand 
operated  hydraulic  pump  with  a  pressure  gauge  to  monitor  the  applied  load. 

When  attempting  to  apply  significant  torque  loads  to  the  device  it  was 
found  that  slipping  occurred  between  the  wood  and  the  shoes  because  of  their  failure 
to  key  by  mechanical  identation.  Methods  of  enhancing  the  coefficient  of  friction 
between  the  two,  such  as  interposing  a  thin  rubber  sheet,  proved  unsuccessful.  These 
difficulties  were  overcome  by  using  sand  as  the  locking  mechanism,  an  approach  which 
was  taken  with  reluctance  because  this  prototype  cylinder  incorporated  no  design 
features  for  excluding  sand  from  the  precision  sliding  surfaces. 

Prior  to  exposing  the  cylinder  to  the  sand,  the  gaps  between  the  shoes  were 
packed  with  mastic  sealin'-1  strip.  The  cylinder  was  supported  concentrically  within 
the  dummy  borehole  and  a  simple  annular  seal  was  made  with  modelling  clay  at  the  bot¬ 
tom  of  the  borehole  tube,  between  the  wood  lining  and  the  shoes.  Dry  sand  was  poured 
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down  the  annular  gap  and  manually  compacted.  The  shoes  were  then  radially  expanded 

to  key  into  the  sand  and  further  compact  it.  With  this  system  it  was  found  that 

slipping  occurred  at  a  torque  which  was  proportional  to  the  radial  jacking  pressure. 

At  the  maximum  hydraulic  pressure  which  could  be  applied  to  the  shoes  (18.5  MPa, 

an  effective  radial  pressure  on  the  sand  of  0.49  MPa)  slipping  occurred  at  an  applied 
3 

torque  of  5  x  10  Nm. 

Short  of  setting  the  cylinder  in  concrete  or  epoxy-resin,  the  testing 
has  demonstrated  that  the  cylinder  is  well  able  to  operate  over  the  range  of  applied 
loads  which  will  normally  be  encountered  in  service  and  the  sufficient  confidence 
has  been  gained  to  extrapolate  that,  where  a  positive  key  can  be  made  with  the  soil, 
the  device  will  operate  satisfactorily  at  its  maximum  torque. 
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7.0  CONCLUSIONS  AND  RECOMMENDATIONS 

A  concept  has  been  proposed  for  an  in-situ  measurement  device 
designed  to  determine  the  soil  secant  modulus  as  a  function  of  strain  for 
levels  up  to  failure.  This  device  first  bores  itself  into  the  soil  with 
minimal  disturbance;  it  then  lays  itself  into  the  soil  and  at  the  same 
time  re-imposes  the  lateral  stresses  by  radially  expanding  a  number  of 
axial  shoes  against  the  sides  of  the  bored  hole.  A  harmonic  rotation  of 
increasing  amplitude  is  then  applied  and  the  torque  and  corresponding 
angular  rotation  are  measured  to  determine  the  stiffness  characteristics 
of  the  soil.  The  apparatus  has  been  named  the  Borehole  Shear  Device  (BSD). 

Based  on  the  authors'  assessment  of  the  state-of-the-art,  such  a 
device  is  presently  required  for  producing  reliable  soil  properties  for  site 
characterisation  and  subsequent  analysis.  Existing  procedures  for  determining 
the  soil  modulus  over  the  required  range  of  strains  suffer  important  drawbacks 
past  attempts  to  produce  a  device  able  to  operate  in-situ  over  the  complete 
range  of  strains  have  not  been  totally  successful. 

This  report  summarises  the  research  efforts  conducted  during  the 
first  phase  of  work  orientated  towards  production  of  an  operational  3SD. 

The  activities  undertaken  were: 


A  literature  search  and  critical  assessment  of  the  state-of-the-art 
Analytical  and  numerical  studies  to  predict  the  characteristics 
of  the  BSD  operation  and  the  effects  on  the  soil;  also  to  aid  the 


configuration  of  the  test  and  the  instrumentation  required. 


■ 
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The  specification  of  a  procedure  to  derive  the  value  of  the 
soil  secant  shear  modulus  as  a  function  of  the  shear  strain  level, 
based  on  the  values  of  the  torque  and  the  rotation  measured  in  the 
field.  This  procedure  attempts  to  minimise  possible  errors  which 
might  result  from  the  presumed  constitutive  behaviour  of  the  soil. 
A  conceptual  design  of  the  BSD,  its  instrumentation  requirements 
and  mode  of  operation. 

Detailed  design,  construction  and  laboratory  testing  under  the 
design  torque  of  the  radially  expanding  cylinder. 


All  of  the  above  tests  have  been  duly  completed.  As  a  consequence 
of  this  work,  a  number  of  conclusions  have  been  reduced.  They  can  be  summarised 


thus: 


From  the  theoretical  studies  conducted  and  t-.»-.  analysis  of  the 
hardware  reduced,  the  construction  and  operation  of  the  BSD  seems 
to  be  feasible  and  desireable. 

A  single-point  measurement  technique  appears  more  adequate,  raxdier 
than  uphole,  crosshole  or  other  alternative  test  configurations.  This 
refers  only  to  normal  operation  of  the  device;  it  will  probably  be 
convenient  though  to  install  extra  instrumentation  during  initial 
field  trials  with  the  purpose  of  confirming  particular  details  of  the 
operat:  ?n  of  the  BSD. 

Two  concepts  have  been  provided:  a  double-cylinder  and  a  single-cylinder 
operation.  Altho  jh  no  decisions  on  this  matter  were  required  during 
the  present  phase  f  the  work,  the  latter  seems  to  be  more  advantageous 
both  technically  and  economically;  hence  the  full  development  of  the 
single-cylinder  BSD  is  that  proposed  for  future  work. 
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4.  At  the  design  frequencies  of  the  BSD  (about  10  Hz) ,  wave-lengths 
are  large  compared  to  all  the  dimensions  of  the  test.  Also,  contrary 
to  the  initial  concept  and  the  test  in  which  it  was  expected  to  impose 
major  e.  strains,  the  main  component  of  the  strain  was  found  to 

UZ 

be  e^0  and  the  shear  modulus  corresponds  to  this  type  of  deformations. 

5.  ‘  ll  is  important  to  minimise  the  disturbance  of  the  soil  around 

the  BSD.  As  a  consequence,  the  latter  has  been  provided  with  a 
self-boring  tool. 

6.  A  procedure  has  been  devised  for  deriving  other  shear  modulus  as 
a  function  of  strain  from  torque  and  rotation  measurements.  The 
procedure  is  based  in  progressive  hyperbolic  approximations  and 
has  oeen  successfully  tested  for  hyperbolic  and  non-hyperbolic 
soil  constitutive  laws. 

7.  It  is  also  concluded  from  the  analysis  that  stiff er  or  softer 
layers  close  to  the  BSD  will  affect  the  measurements  by  very  small 
amounts.  In  soils  of  low  sensitivity,  the  effect  of  disturbance  around 
the  borehole  will  also  be  small  due  to  the  self-boring  placement 
and  to  the  fact  that  i-he  lateral  normal  stresses  ar-.  theoretically 
recovered. 

8.  Close-form  sol  i cions  reported  include:  full  dynamic  long-cylinder 
solutions  for  several  constitutive  laws;  full  dynamic  long-cylinder 
solutions  for  elastic  soil  (with  or  without  damping) ;  and  torque- 
rotation  relationships  for  the  case  of  a  finite  cylinder  at  finite 
depth  in  elastic  soil. 

9.  Numerical  solutions  were  based  on  the  finite  element  method;  an 
equivalent  linear  procedure  was  implemented  to  produce  solutions 
in  the  frequency  domain.  Non-reflecting  boundaries  of  the  Lysmer- 
Waas  and  Lysmer-Kuhlemeyer  types  were  included  in  the  model.  Excellent 
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agreement  was  shown  with  existing  close-form  solutions-  Numerical 
solutions  were  used  for  all- non-linear  dynamic  cases  and  some 
elastic  ones.  At  the  frequencies  of  interest,  the  real  part  of 
the  dynamic  solutions  essentially  coincides  with  the  static  one; 
this  is  particularly  important  for  constructing  the  procedure 
for  derivation  of  the  shear  modulus  from  the  field  data. 

10.  Displacement  and  strain  distributions  have  been  obtained  for 

the  cases  of  interest.  Although  those  have  not  provided  quantitative 

data  in  usable  form,  they  have  produced  essential  qualitative  information, 

such  as  the  relative  importance  of  i  „  versus  e.  strains,  the 

r0  9z 

progressively  faster  decrease  of  the  strain  level  away  from  the  BSD 
as  the  soil  becomes  more  non-linear,  etc. 

11.  The  cylinder  and  radially  expanding  keys  were  designed  and  built. 
Laboratory  expansion  and  torsion  tests  showed  this  part  of  the  BSD 
to  perform  adequately  for  the  purposes  intended  up  to  the  design 
radial  pressures  and  torques. 


It  is  recognised  that  the  work  reported  here  constitutes  only 

a  first  step  towards  the  solution  of  a  difficult  problem.  In  order  to  continue 

♦ 

this  effort  in  the  most  appropriate  manner,  several  recommendations  are  proposed 
he.-e  based  on  the  results  and  experience  achieved  to  this  point: 


The  work  should  continue  towards  production  and  operation  of  the  BSD. 
No  more  theoretical  work  is  deemed  necessary  at  this  point,  at  least 
until  the  device  starts  undergoing  field  tests. 

It  is  recommended  that  the  single-cylinder  concept  be  preferred 
instead  of  the  double-cylinder  one  because  of  the  reasons  stated 
in  the  report. 
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3. 


Due  to  the  unknowns  present  in  all  development  work  it  is  proposed 
that  on  the  basis  of  experience  gained  to  date  the  remaining  work 
be  tackled  in  two  separate  phases:  to  design  and  manufacture 
the  mechanical  system  relating  to  the  torque  loading  so  that  the 
concept  can  be  tested  and  developed  under  more  realistic  conditions 
and,  the  final  phase,  to  manufacture  the  self -boring  device  and 
fit  all  instrumentations  so  that  full  field  trials  can  be  performed. 


14S 


8.  REFERENCES 


1  Schnabel,  P.B.,  Lyssner,  J.,  Seed,  H.B.,  (1972)  “SHAKE  a  Computer  Program 
for  Earthquake  Response  Analysis  of  Horizontally  Layered  Sites",  Report 
No.  EERC  72-12,  University  of  California,  Berkeley,  Dec.  1972 

2  Seed,  H.B.,  Idriss,  I.M,  (1970)  "Soil  Moduli  and  Damping  Factors  for 
Dynamic  Response  Analysis",  Report  No.  EERC  70-10,  University  of 
California,  Berkeley,  Dec.  1970. 

3  Stokoe,  K.H.,  Arnold,  E.J.,  Hoar,  R.J.,  (1978)  "Development  of  a 
Bottom-hole  Device  for  Offshore  Shear  Wave  Velocity  Measurement", 

Offshore  Tech.  Conf.  paper  3210,  May  1978 

4  Stokoe,  K.H.  Richart,  F.E.,  (1973)  "Shear  Moduli  of  Soils,  In-Situ 
and  from  Laboratory  Tests”,  Proc.  of  the  Fifth  World  Conf.  on 
Earthquake  Eng.,  Vol.  1,  Rose,  Italy,  (1973),  pp.  356-359 

5  Anderson,  D.G. ,  Stokoe,  K.H. ,  (197S)  "Shear  Modulus,  A  Tine-Dependent 
Material  Property",  "Dynamic  Soil  ar.d  Rock  Testing  in  Field  and 
Laboratory  for  Seisaic  Studies" ,ASTM  STP,  (1978) 

6  Anderson,  D.G. ,  Woods,  R.D. ,  (1976)  "Tine-Dependent  Increase  in  Shear 
Modulus  of  Clay",  Journal  of  the  Geotechnical  Eng.  Div. ,  ASCE,  May 
1976,  pp.  525-537 

7  Tinoshenko,  S.P.,  Goodier,  J.H.,  (1951)  "Theory  of  Elasticity", 

Third  Edition,  McGraw  Kill,  New  York. 

8  Ewing,  W.M. ,  Jardetzkv,  W.S.,  Press,  T. ,  (1957)  "Elastic  Waves  in 
Layered  Media",  McGraw-Hill,  New  York 

9  Duke,  H. ,  (1969),  "Techniques  for  Field  Measurement  of  Shear  Wave 
Velocity  in  Soils”,  Fourth  World  Conf.  in  Earthquake  Eng.,  Santiago, 

Chile,  1969 

10  Richart,  F.E. ,  Hall,  J.R.,  Woods,  R.D. ,  (1970)  "Vibrations  of 
Soils  and  Foundations",  Prentice  Hall,  1970 

11  Mooney,  H.H.,  (1974)  "Seismic  Shear  Waves  in  Engineering",  Jour,  of  Soil  Hech 
and  Found.  Div.,  ASCE,  Aug  1974,  pp.  905-923 

12  Schwarz,  S.D.,  Conweil,  F.R.,  (1974)  "A  Technique  for  In-Situ  Measurement 
of  Shear  Wave  Velocities  (Vs)  for  Deep  Marine  Foundations” ,  Offshore 
Tech.  Conf.,  paper  2014,  May  1974 

13  Lance,  G. ,  (1974)  "Seisaologisk  Metod  for  In-Si tubes tamaning  av  Jord- 
och  Bergdynaaiska  "  as-.-ter" ,  3jerking  Ingenjorsbyra,  Uppsala,  Sweden 


«ae@s«@ss^ 


Kurzeme,  M. ,  (1971)  "In-Situ  Investigations  Using  SH-Waves" ,  J.  of  the 
Soil  Mech.  and  Found.  Div. ,  ASCE,  Feb.  1971,  pp.  341-356 

Jones,  R. ,  (1960)  "Measurement  and  Interpretation  of  Surface  Vibrations 
on  Soil  and  Roads",  Highway  Research  Board  Bulletin,  1960 

Jones,  T.R. ,  (1962)  "Surface  Wave  Technique  for  Measuring  the 
Elastic  Properties  and  Thickness  of  Roads:  Theoretical  Development", 

Brit.  J.  Appl.  Physics,  1962,  Vol.  13,  pp.  21-29 

Fowler,  J. ,  (1972)  "An  Evaluation  of  an  Existing  Procedure  for  Determining 
Shear  Moduli  at  Depths  by  In  Situ  Vibration  Technique",  Technical  Report 
S-72-5,  WES,  Vicksburg,  Missisippi,  April  1972 

Ballard,  R.F.,  McLean,  F.G.,  (1975)  "Seismic  Field  Methods  for  In  Situ 
Moduli",  Proc.  of  the  Conf.  on  In  Situ  Measurement  of  Soil  Properities, 
ASCE,  June  1975,  Vol'.  1,  pp.  121-150 

Stokoe,  K.H.,  Woods,  R.D. ,  (1972)  "In-Situ  Shear  Wave  Velocity  by  the 
Cross-Hole  Method",  J.  of  Soil  Mech.  and  Found.  Div.  ASCE,  May  1972, 
pp.  443-460 

Swain,  R.J.,  (1962)  "Recent  Techniques  for  Determination  of  In-situ  Elastic 
Properties  and  Measurement  of  Motion  Amplification  in  Layered  Media", 
Geophysics,  Vol.  27,  No.  2,  April  1962,  pp.  237-241 

Tanimoto,  K. ,  Kurzeme,  M. ,  (1973)  Discussion  of  Stokoe  et  al,  J.  of 
Soil  Mech.  and  Found.  Div.,  ASCE,  April  1973,  pp.  351-353 

Auld,  B. ,  (1977)  "Cross-Hole  and  Down-Hole  vs  by  Mechanical  Impulse", 
Journal  of  the  Geotechnical  Eng.  Div.,  ASCE,  Dec.  1977,  pp.  1381-1398 

Ballard,  R.F.,  Leach,  R.E.,  (1969)  "Development  of  a  Vibropacker 
System  for  Inducing  Polarize  Shear  Waves  and  Compression  Waves  at 
Depth",  S-69-30,  WES,  Vicksburg,  Missisippi,  July  1969 

Miller,  R.P. ,  Brown,  F.R.,  (1972)  "Shea--  Modulus  Determination  of  Soils 
by  In-Situ  Methods  for  Earthquake  Engineering",  Proc.  of  the  Int.  Conf. 
of  Micromization  for  Safer  Construction,  Research  and  Application,  Seattle, 
Washington,  Nov.  1972,  pp.  319-335 

Miller,  R.P. ,  Troncoso,  J.H.,  Brown,  F.R.,  (1975)  "In-Situ  Impulse  Test 
for  nynamic  Shear  Modulus  of  Soils",  Proc.  of  the  Conf.  on  In  Situ 
Measurement  of  Soil  Properties,  ASCE,  June  1975,  Vol  1,  pp.  319-335 


Troncoso,  J..-  (1975)  "In-Situ  Impulse  Test  for  Determination  of  Soil 
Shear  Modulus  as  a  Function  of  Strain",  Ph.D.  Thesis,  Dept,  of 
Civil  Eng.,  University  of  Illinois  at  Urbana-Champaign,  1975 

Boyer,  G.R.,  Oien.  M.A. ,  (1972)  "Steady  Motion  of  a  Torsional  Oscillator 
Clamped  in  a  Bore! ole",  Journ.  of  Sound  and  Vibration,  (1972)  23  (2),  pp.  17u-186 


151 


28  Lysmer,  J.,  Waas,  G. ,  (1972)  "Shear  Waves  in  Plane  Infinite  Structures", 

Journ.  of  Eng.  Mech.  Div. ,  ASCE,  Feb.  1972,  pp.  85-105 

29  Lysmer,  J.,  Kuhlemeyer,  R.L.,  (1969)  "Finite  Dynamic  Model  for  Infinite 
Media",  Journal  of  Eng.  Mech.  Div.,  ASCE,  Aug.  1972,  pp.  859-877 

30  Hughes,  J.M.O. ,  Wroth,  C.P.,  Windle,  D. ,  (1977)  "Pressuremeter  Tests 
in  Sands",  Geotechnique  27,  No.  4,  po  455-477 

31  Windle,  D. ,  Wroth,  C.P.,  (1977)  "The  Use  of  Self-Boring  Pressuremeter 

to  Determine  the  Undrained  Properties  of  Clays",  Ground  Engineer. ng.  Sept.  1977 


i 


I 


4 


APPENDIX  A 


vertical  force  developed  by  friction 
tip  force 

undamped  resonant  frequency  of  the  system 
influence  function 

influence  function  for  the  elastic  soil 
influence  function  for  the  hyperbolic  soil 
gravitational  constant 
height  of  the  BSD's  snoe 
height  of  the  self  boring 

Hankel  function  of  the  first  kind  and  order  one 
Hankel  function  of  the  second  kind  and  order  one 
Hankel  function  of  the  second  kind  and  order  two 
complex  constant 
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lumped  moment  of  inertia  of  the  coupling 
tube  and  drive  system 

Bessel  function  of  the  first  kind  and  order  zero 

Bessel  function  of  the  first  kind  and  order  one 

coefficient  of  lateral  earth  pressure 

torsional  stiffness  of  the  hollow  tube 

length  of  the  hollow  tube 

lumped  mass  in  a  finite  element  node 

maximum  reactive  forces  power 

peak  power  at  resonance 

resonance  maqnification  of  the  system 

radial  coordinan.e 

radius  of  the  borehole 

internal  radius  of  a  hollow  tube 

external  radius  of  a  hollow  tube 

mean  radius  of  the  coupling  tube 

thickness  of  the  coupling  tube 

normalized  radius  of  the  BSD 

time 

time  constant  of  the  exponential  envelope 
of  the  response 

torque  per  unit  length 

torque  per  unit  length  in  a  finite  cylinder 
with  elastic  soil 

torque  per  unit  length  in  an  infinite 
cylinder  with  elastic  soil 

torque  per  unit  length  in  a  finite  cylinder 
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torque  per  unit  length  in  a  finite  |j 

cylinder  with  hyperbolic  soil  M 


torque  per  unit  length  in  an  infinite 
cylinder  with  hyperbolic  soil 

torque  per  unit  length 

torque  for  an  infinite  cylinder 
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maximum  torque  applied  to  the  coupling 
cylinder 

tangential  displacement 
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tangential  displacement  at  a  FE  node 

tangential  displacement  at  the  borehole 
wall 

radial  displacement 
tangential  displacement 
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velocity  of  compressional  waves 

velocity  of  Rayleigh  waves 

velocity  of  shear  waves 

Bessel  function  of  the  second  kind  and 
order  zero 

Bessel  function  of  the  second  kind 
and  order  one 

vertical  coordinate 

angular  twist  over  the  height  of  the 
cylinder 

friction  angle  in  the  soil-shoe  interface 
shear  strain 
equivalent  strain 

shear  strain  in  the  borehole  wails 
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shear  strain  components 

coefficient  of  linear  viscous  damping 
rotation  angle 

rotation  angle  at  the  borehole  wall 

maximum  rotation  angle  in  the  borehole 

internal  friction  angle  of  the  soil 

Lame  constant 

secant  shear  modulus 

secant  shear  modulus  of  the  tube 

shear  modulus  of  the  soil  in  the  finite 
element  zone 

shear  modulus  of  the  soil  at  low  strain 
(elastic) 

reference  shear  modulus 

Poisson's  ratio 

density  of  ''he  material 

density  of  the  soil 

density  of  water 

maximum  lateral  jacking  pressure 

radial  normal  stress 

effective  confining  pressure 

effective  lateral  stress 

reference  effective  confining  pressure 

shear  stress 

limiting  shear  stress 

shear  stress  at  the  cylinder-soil  interface 
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NON-REFLECTING  BOUNDARIES 


S3 


As  mentioned  in  Chapter  4,  two  different  non-reflecting  boundaries 


have  been  used: 


A  permanent  non-reflecting  boundary  of  the  Lysmer-Waas  type  has 
been  provided  at  the  cylindrical  boundary. 

An  optional  non-reflecting  boundary  of  the  Lysmer-Kuhlemeyer  type 
was  implemented  at  the  upper  boundary. 


In  either  case,  the  purpose  of  such  boundaries  is  to  guarantee 
that  the  phenomena  predicted  inside  the  finite  domain  modelled  are  those 
which  would  be  observed  if  the  domain  modelled  extended  to  infinity.  That 
is  to  say,  that  the  relatively  arbitrary  location  of  the  boundaries  does  not 
affect  in  any  way  the  calculations. 


6.1  THE  LYSMER-WAAS  AXISYMMETRIC  BOUNDARY 


The  basic  assumptions  embodied  in  the  axisymmetric  version  of  this 
boundary  are: 

The  soil  properties  outside  the  boundaries  are  linear  and 
change  only  with  depth. 

The  particle  velocity  takes  place  only  in  the  0  direction. 

Waves  arriving  at  the  boundary  are  locally  cylindrical. 


The  implementation  of  the  Lysmer-Waas  boundary  requires  the  determination 
of  the  following: 

The  stiffness  matrices  at  the  boundary  elements,  that  is,  the 
relation  between  forces  and  displacements  at  the  boundary  nodes 
accounting  for  the  fact  that  the  soil  extends  all  the  way  to  infinity. 


The  forces  introduced  at  the  boundary  of  the  finite  domain 
due  to  the  presence  of  the  external  1  which  is  not 
modelled  by  the  mesh. 

The  following  derivations  will  be  pursured  in  parallel  with  those 

presented  in  the  original  plane-strain  version  of  this  boundary  (Lysmer 
28 

and  Waas,  1972) .  Figure  B  1  shows  the  finite-element  mesh. 


Due  to  one  of  the  basic  assumptions  of  the  FE  method,,  the  displace¬ 
ment  at  any  point  of  a  rectangular  axisymmetric  element  are  (since  o...’’'  cut 
of  plane  displacements  occur) : 


6  (r,z)  =  || d ||  |  6  } 

where  ||d||  =  —  ||  (b-z)  (h-r)  z(h-r)  zr  r(b-z)|| 


(1) 


j  6  }  is  the  matrix  of  the  displacements  at  the  four  corners 
(Fig.  b  2) 

The  strains  can  be  be  derived  from  the  displacements: 
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All  other  strain  conponents  vanish  everywhere. 


The  stress  components  are: 


(3) 


where  p.  is  the  corresponding  shear  modulus. 


The  forces  per  unit  radian  acting  on  the  corners  of  each  element 
can  then  be  related  to  the  nodal  displacements: 
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is  called  the  stiffness  matrix  of 


the  element. 


When  the  above  integration  is  carried  out,  the  stiffness  matrix  can 


be  derived. 
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where  m  =  b  r  ln(l+a) 
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B  .1.1  Stiffness  Matrix  for  External  Elements 

Let  us  now  concentrate  in  the  region  outside  the  boundary.  In  any 
two  elements  (A  and  B  in  Fig.  B . 3) . 
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The  coefficients  D  and  D'  are  easily  obtained  from  the  solution  for 
cylindrical  Love-wave  propogation  in  the  positive  r-direction: 


u=u  (z)  H  (2)(kr)e1Ut  (7) 

o  1 

where  uq  (z)  is  an  amplitude  depending  on  depth 

k  is  the  wave  number  (w/v  ) 

s 


is  the  Hankel  function  of  second  kind  and  order  1. 


Then: 


.  .  all  values  of  h: 

D  =  1  +  khF  +  (kh)2S 
d'=  1  -  khF  +  (kh)2S 
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In  order  to  derive  the  forces  in  the  two  nodes  common  to  A  and 
B  (nodes  1  and  2) : 
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where  | J J j  stands  for  the  rth  row  of  the  stiffness  matrix  of  element  A 


By  performing  the  operations  indicated  above  and  taking  h->0- 
the  following  expression  can  be  obtained  relating  the  force  and  the 
displacements  at  nodes  1  and  2. 
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It  can  be  readily  seen  from  the  recurrence  formulae  between 
Hankel  functions  that: 
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which  substituted  above,  yields: 
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B  .1.2  External  Boundary  Forces 

Consider  again  Fig.  B.3  where  the  boundary  now  splits  the  two 
zones  A  and  B  (Fig.  :i  .4)  . 

From  equation  (3) ,  for  z  =  0  and  r  =  0. 
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where  D  =  1  +  Kh  F  since  second  order  terms  can  now  be  neglected. 


Equation  (15)  gives  the  stress  at  the  top  of  layer  j.  At  the 


bottom: 
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The  forces  per  unit  radian  at  both  nodes  can  then  be  obtained: 
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B.1.3  Concluding  Remarks 


Equations  (12)  and  (17)  provide  '.ne  stiffness  matrix  of  the  boundary 
elements  and  the  equivalent  boundary  forces,  respectively.  The  problem  can 


then  be  completely  solved  with  the  aid  of  those  expressions.  It  is  interesting 

2 


to  notice  that,  although  the  eigenvalues  (k  )  of  the  resulting  problem  are  all 
real,  some  of  them  are  positive  and  some  are  negative.  The  positive  eigenvalues 
yield  real  wave  numbers  which  correspond  to  travelling  waves;  the  sign  of.  the 
square  root  must  be  selected  so  that  the  wave  propagates  away  from  the  axis. 
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Negative  eigenvalues  generate  imaginary  wave  numbers,  which  correspond  to 
standing  waves;  the  sign  of  the  square  root  should  be  taken  such  that  the 
standing  wave  decays  towards  infinity. 


B  .2  THE  LYSMER-KUHLEMEYER  PLANE  BOUNDARY 


The  full  discussion  of  this  boundary  can  be  found  in  the  original 

29 

paper  by  Lysmer  and  Kuhlemeyer  (1969) .  The  reader  is  referred  there  for 
details  on  the  characteristics  and  behaviour  of  the  boundary. 


For  the  purpose  of  the  present  report  it  is  sufficient  to  state 
that  the  boundary  consists  of  a  series  of  dashpots  connected  to  the  boundary 
nodes.  The  viscous  constants  of  these  dashpots  are  taken: 


C  =  p  V  A 


where 


p  is  the  density. 


V  is  the  velocity  of  propagation  of  the  radiated  pulse  (shear  waves 
in  this  case) . 

A  is  the  area  lumped  in  the  boundary  node  to  which  the  dashpot  of 
constant  C  is  connected. 


Da.-hpots  chosen  in  such  way  dissipate  exactly  all  the  energy 
carried  by  waves  in  the  case  of  one-dimensional  propagation.  For  this 
case  it  constitutes  a  perfect  non-reflecting  boundary.  In  the  case  of  more 
than  one  dimension,  some  of  the  incoming  energy  is  reflected  back  into  the 
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The  following  figure,  taken  from  the  referenced  paper,  shows  the 
ratio  of  reflected  to  incident  energy  ratio  as  a  function  of  the  angle  of 
incidence. 
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FIG.  B.5.  EFFECTIVE  ENERGY  RATIO  FOR  INCIDENT  P-WAVE 
(AFTER  LYSMER  AND  KUHLEMEYER,  1972) 


Given  the  facts  that  in  the  finite  element  model  proposed  the 
angle  of  incidence  is  never  small  and  that  the  effects  of  the  BSD  vanish 
rapidly  above  and  below  the  depth  of  operation,  the  proposed  form  of  viscous 
boundary  is  thought  to  be  sufficiently  accurate  to  simulate  the  radiation  of 
energy  across  the  norizontal  boundaries  of  the  model. 


APPENDIX  C 


MISCELLANEOUS  CALCULATIONS 


C  . 1  ELASTIC  DISPLACEMENT  DISTRIBUTION 


From  Chapter  3,  equation  (9)  is: 
2 

.  t  r  “ 

du  _  u  o  o 
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where  the  partial  has  been  replaced  by  an  ordinary  derivative.  This  equation 
is  subjected  to  the  boundary  condition: 
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The  general  solution  to  ti  .  above  linear  equation  can  be  obtained 

as: 
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The  boundary  condition  requires: 
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Hence  the  final  solution  is: 


APPENDIX  D 


GENERAL  DESCRIPTION  OF  FINEL 

I  INTRODUCTION 

FINEL  is  a  finite-element  system  which  has  been  designed  to  meet 
the  requirements  of  teaching  and  research  establishments.  It  has  been 
developed  over  a  period  of  years  by  the  Department  of  Aeronautics,  Imperial 
College,  London.  The  system  has  most  of  the  features  found  in  large  general 
purpose  finite-element  programs,  but  also  has  additional  facilities  allowing 
the  user  to  add  features  pertinent  to  his  own  requirements.  In  adding  these, 
the  interfacing  between  the  user  and  the  program  is  a  minimum  and  occurs  in 
such  a  way  that  it  is  not  possible  for  any  user  to  alter  permanently  the  system. 

The  FINEL  input  has  been  designed  to  be  of  a  general  form  such  that 
it  is  capable  of  being  adapted  to  suit  any  problem  which  can  be  solved  by 
the  finite-element  method.  This  ensures  that  it  is  completely  compatible 
with  the  open-endedness  built  into  the  system.  The  input  is  also  of  such 
a  form  that  it  is  readily  understood,  making  the  system  easj  to  use. 

A  further  feature  of  the  system  is  that  it  has  been  specifically 
designed  to  fit  into  small  core,  shorttime  partitions  within  a  computer, 
but  with  the  user  able  to  expand  core  allociation  at  run  time  for  large 
analyses.  This  feature  is  essential  for  a  teaching  and  research  orientated 
program  since  a  large  amount  of  its  usage  will  involve  debugging,  requiring 


a  rapid  turnaround  time. 


D.2  SPECIAL  FEATURES  OF  FINEL 


FINEL  can  be  used  in  the  seine  way  as  other  finite  element  systems 
for  analysing  structures  and  other  problems.  In  addition  to  this  it  has 
particular  features  which  enable  it  to  be  used  in  a  rather  more  general 
fashion.  The  most  important  of  these  are: 

a)  It  has  a  highly  modular  structure,  including  sub-modules 
(termed  libraries  within  FINEL) . 

b)  A  general  form  of  free  fc_-mat  data  input  using  the  command 
and  label  concept  to  define  a  block  of  data. 

c)  A  central  data  filing  system,  forming  a  problem  data  base. 
This  has  a  two  level  structure  of  files  and  pages  within 
files.  Only  one  subroutine  is  used  to  control  the  data  base 
and  this  of  such  a  form  that  it  is  easy  for  the  user  to 
employ. 

a)  There  are  no  assumptions  built  into  the  system  that  fix  *t 
to  a  particular  type  of  analysis.  This  applies  tc  both 
the  input  formats  and  the  execution  sequence. 

During  the  analysis  carried  out  for  the  BSD  project,  several 
features  vere  added  to  the  program: 

An  optional  axisymmetric,  plane,  non-reflecting  boundary  of  the 

Lysroer-Kuhlemeyer  type. 

An  axisymmetriu ,  cylindrical,  non-reflecting  boundary  of  the 

Lysaer-Waas  type. 


A  procedure  to  incorporate  “equivalent- linear"  calculations  for 

frequency-domain  solution  of  non-linear  problems. 

D.3  EXECUTION  SEQUENCE  OF  A  FINEL  RUM 

The  basic  execution  sequence  of  a  FINEL  is  shown  in  Figure  D.l. 

If  the  user  is  adding  facilities  to  the  basic  program  then  the  routines 
describing  these  are  compiled  first.  These  are  written  in  FOETPAN. 

A  large  amount  of  data  checking  is  carried  out  at  the  preprocessor 
stage.  An  interactive  form  of  this  is  available  for  on  line  data  input  and 
correction. 

The  finite-element  analysis  is  carried  out  in  the  FINEL  main 
program.  This  is  structured  in  a  highly  modular  form.  Each  module  is 
only  entered  through  the  control  program.  The  modules  which  are  called, 
and  their  calling  sequence,  are  basically  defined  at  runtime  but  there 
are  facilities  for  altering  this  depending  upon  intermediate  results. 

This  allows  iterative  and  non-linear  analysis  to  be  performed. 

The  modular  structure  has  various  advantages.  It  enables  simple 
overlay  sequences  to  be  set  up,  allowing  FINEL  to  be  run  on  a  relatively 
small  machine.  More  importantly,  for  a  research  orientated  system,  modules 
can  be  added  at  runtime.  This  allows  the  user  to  modify  the  system  to  suit 
his  own  requirements  without  actually  altering  any  part  of  the  system.  The 
modules  are  written  as  essentially  independent  programs  (Fig.  D.2). 
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The  modular  concept  is  repeated  at  a  lower  level  in  the  form  of 
FINEL  libraries.  These  have  the  access  structure  snown  in  Figure  d  -3. 
The  libraries  contain  a  series  of  standard  entries  which  can  always  be 
used.  There  are  also  facilities  for  the  user  to  add  to  his  own  entries 
are  runtime. 


FINEL  contains  six  libraries,  these  being: 


a)  Element  library. 


b)  Load  library. 

c)  Mesh  generation  library. 

d)  Plotting  library. 

e)  Program  library. 


f)  Module  library. 


All  libraries  have  the  same  form  and  are  used  in  the  same  way.  For 
teaching  purposes  students  can,  typically,  add  their  own  elements  to  gain 
practical  experience  in  F.E.  theory.  The  coding  required  for  this  is  minimal. 


D.4  FINEL  DATA  INPUT 


Data  is  inputted  to  FINEL  in  free  format.  Each  blcck  of  data 
consists  of  a  command  word,  a  label  and  numerical  data.  The  commands  and 
labels  are  usually  English  words,  making  the  input  itself  explanatory  in 
most  cases.  This  form  of  input  was  chosen  since  it  is  adaptable  to  any 
type  of  problem.  The  command  word  is  usually  defined  by  F.E.  theory, 
whereas  the  label  can  be  chosen  to  suit  the  problem  nomenclature. 
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A  simple  stress  analysis  problem,  together  with  the  FINEL  input 
is  shown  in  Figure  d.4. 

D. 5  THE  F1NEL  FILING  SYSTEM 

Data  is  transferred  within  FINEL  either  in  labelled  common  areas 
for  small  (single  words)  amounts  of  data  or  through  a  direct  access  filing 
system  (data  base)  for  large  blocks  of  data.  The  data  base  comprises  a  two 
level  system  of  files  and  pages  with  files.  This  is  illustrated  in  Figure 
0.5. 

A  set  of  files  are  dedicated  to  storing  blocks  of  data  that  are 
always  generated,  by  the  nature  of  the  F.E.  method.  Other  files  are 
available  to  the  user  for  intermediate  scratch  use. 

Any  operation  on  the  filing  system  is  carried  out  through  a  single 
routine  called  ZFILE.  This  is  called  as: 

CALL  ZFILE  (VAR , NVAK , NFILE , NPAGE , IV7TRD ) 
where  VAR  is  the  first  word  address  for  the  data  transfer 

NVAR  is  the  number  of  words  transferred  on  page  NPAGE  of  file  NFILE. 

IWTRD  is  set  to  1  for  writing  and  2  for  reading. 

The  filing  system  has  proved  to  be  easy  to  use  and  efficient  in  operation. 
Since  all  data  base  operations  are  controlled  through  a  single  routine  it  becomes 
easy  to  transfer  the  system  to  other  computers.  The  filing  routine  can  then  be 
written  separately  for  each  machine  to  give  optimum  efficiency.  The  rest  of  FINEL 
is  in  standard  ASA  FORTRAN  so  that  no  other  changes  have  to  be  made  when  transferring 


to  other  machines. 
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FiG  D.4.  EXAMPLE 
OF  FINEL  INPUT 


D . 6  ASSUMPTIONS  WITHIN  FINEL 


Every  attempt  has  been  made  to  make  FINEL  free  from  any  restrictive 
asumptions.  The  data  input  allows  problems  to  be  specified  in  their 
natural  form.  Elements  can  be  used  to  define  up  to  ten  types  of  coefficient 
matrices,  although  only  up  to  four  can  be  used  in  any  one  run.  It  is 
assumed  that  every  node  on  an  element,  for  a  given  matrix  type,  has  the 
same  number  of  freedoms  per  node  and  this  is  not  always  so,  however  in  such 
cases  dummy  freedoms  can  do  specified  without  loss  of  efficiency.  There 
is  no  restriction  upon  the  number  of  freedoms  per  node,  or  the  types  of 
freedoms.  The  pre-processor  automatically  allocates  core,  optimizing  array 
sizes  to  suit  the  problem.  The  largest  single  problem  is  dictated  by  the 
requirement  to  hold  a  completed  load  vector  in  core,  but  the  use  of  sub¬ 
structuring  allows  any  size  problem  to  be  solved. 

P.7  SOLUTION  PROCEDURE  IN  THE  PRESENT  FINEL  VERSION 

The  calculation  procedure  employed  in  FINEL  is  based  on  the 
technique  of  modal  decomposition  which  is  the  general  scheme  normaly  used 
in  other  structural  dynamics  frequency  domain  computer  routines.  The 
principles  of  the  technique  are  illustrated  in  Figure  D.6  in  which  the  logical 
sequence  of  the  main  operations  are  shown. 

To  expand  upon  this  illustration ,  it  was  shown  in  Appendix  B 
that  the  values  of  displacement,  strains  and  stresses  at  any  point  in  an 
element  can  be  derived  as  a  function  of  the  displacements  of  the  nodes  of 
that  element: 

6  (r,z)  =  ||  d  ]{  {  6  }  (1) 


185 


start; 


'READ  INPUT  DATA: 

GEOMETRY  AND  PROPERTIES/ 


SET  UP  INITIAL  PROBLEM 


READ  INPUT  FREQUENCY 


PRINT  RESULTS 

3Sjcwvergence  re/ 

HTyes 

[  END  1 


FIG.  D.6.  GENERAL  FLOW-CHART  FOR  FINEL  VERSION  IN  THE 
BSD  PROJECT. 


{c|  =  [D]  (6|  (2) 

|  T  |  -  Uj[D]  1  S  |  (3) 

where  INI  and  are  the  matrices  introduced  in  the  equations 

(1)  and  (2)  of  Appendix  B 

| o |  is  a  column  vector  representing  the  displacement  in  the  element 
nodes . 


Likewise,  the  stiffness  and  mass  matrices  for  the  elements  of  the  soil 


can  be  defined: 


H-jC  ”5 

[*]  •  fv  Oj  IUII  idl  dv 


The  soil  will  in  general  display  some  internal  damping  characteristics. 

In  FE  calculations,  the  stiffness  matrix  of  the  soil  is  often  expressed  as  a 
linear  combination  of  the  mass  and  stiffness  matrices  (Rayleigh  damping) .  Such  is 
the  procedure  used  in  FINEL-  However,  preliminary  estimates  showed  internal  damping 
to  be  unnecessary  in  the  BSD  runs.  This  is  due  to  the  small  size  of  the  domain 
modelled  compared  to  the  wave  lengths  induced  by  the  device. 


Radiation  damping  is  however  important.  It  appears  in  the  formulation 
through  the  complex  stiffnesses  associated  with  the  non-reflecting  boundaries,  both 
the  Lysraer  dashpots  and  the  Lysmer-Waas  boundary. 


The  equations  can  then  be  formulated  as: 


+  iuC  -  uj2M  J  j  r  j  =  |  R  | 
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Where  K  is  composed  of  the  soil  stiffness  matrix  plus  the  real 
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a  global  displacement  vector  r  has  been  obtained  it  is  then  possible  to  deduce 
the  strains  in  every  element. 


The  constitutive  non-linearity  is  introduced  approximately  via  the 

equivalent  linear  method.  This  method  has  been  generalized  to  two  dimensions 

by  assuming  that  the  shear  modulus  is  only  a  function  of  position  and  time 

and  is  related  to  the  strain  level  e  via 

e 


u  =  u(ae  ) 
e 


where  a  has  been  taken  as  1 


c  is  the  maximum  value  of  the  shear  strain  at  the  ooint  (in  anv 
e 

direction  at  any  time) . 
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APPENDIX  F 

SET  OF  DRAWINGS  OF  THE  PROTOTYPE  COUPLING 
MECHANISM  AND  TESTING  EQUIPMENT 
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